!----------------------------------------------------------------------- !Main program !----------------------------------------------------------------------- program block_alg ![specification part] implicit none integer, parameter:: M=3 integer ::k,n real(kind=kind(1.0d0)), dimension(1:M,1:M) :: A,B ![execution part] A=reshape ((/12,32,43,4,5,6,7,8,9/), (/3,3/)) ! open (7, file='initiall.dat', status='old', action='read') ! read (7,*) ((a(i ,j), j=1,n), i=1,n) B=myfunc(A,M) do n = 1, M write(*,*) (B(k,n), k=1,M) enddo contains ![your functions] function myfunc(A,N) implicit none integer, intent(in):: N real(kind=kind(1.0d0)), intent(in), dimension(1:N,1:N) :: A integer :: i,j,k real(kind=kind(1.0d0)), dimension(1:N,1:N) :: G, myfunc G=0 !cholesky here извинте, но я тоже ленив myfunc=G end function myfunc end program block_alg |