1module solv_cap 2 integer, parameter, public :: dp = selected_real_kind(5) 3contains 4 subroutine prod0( G, X ) 5 real(kind=dp), intent(in out), dimension(:,:) :: X 6 real(kind=dp), dimension(size(X,1),size(X,2)) :: Y 7 X = Y 8 end subroutine prod0 9 function Ginteg(xq1,yq1, xq2,yq2, xp,yp) result(G) 10 end function Ginteg 11 subroutine fourir(A,ntot,kconjg, E,useold) 12 end subroutine fourir 13end module solv_cap 14