FORTRAN90写的雅克比和高斯赛德尔迭代程序.docVIP

!JACOBI迭代求解方程组 program jacobi implicit none integer :: i,j integer :: k !迭代次数旗标 save k real,parameter :: e=0.001 integer,parameter :: n=4 real :: x(n),y(n),b(n) data b /1.38,-0.34,0.67,1.52/ real :: D real :: a(n,n) !!!千万千万要注意Fortran的数据存放是按照列来存放的，如果变成还按照习惯的数组存放方式输入数据，会导致错误!!! data a /2.52,0.39,0.55,0.23,0.95,1.69,-1.25,-1.15,1.25,-0.45,1.96,-0.45,-0.85,0.49,-0.98,2.31/ write(*,*)*********矩阵A的形式为********* write(*,(1x,4f6.2,/))a forall(i=1:n) x(i)=0 end forall k=1 100 D=0 do i=1