龙格
龙格—库塔法分析lorenz方程课程设计(编辑修改稿)
(1,i)。 end for r=2:n L(r,1)=A(r,1)/U(1,1)。 end for i=2:n for j=i:n z=0。 for r=1:i1 z=z+L(i,r)*U(r,j)。 end U(i,j)=A(i,j)z。 end if abs(U(i,i))eps flag=39。 failure39。 return。 end for k=i+1:n m=0。 for