!ans=0 这个ans的位置错误
X(4) = B(4)/A(4,4)
X(3) = (B(3)-A(3,4)*X(4))/A(3,3)
X(2) = (B(2)-A(2,4)*X(4)-A(2,3)*X(3))/a(2,2)
do j=3,1,-1
ans = 0 !应该在这里
do i=j+1,4,1
ans = ans + A(j,i)*X(i)
end do
X(j) = (B(j)-ans)/A(j,j)
end do
Transpose 发表于 2024-3-28 13:55
[mw_shl_code=fortran,true] !ans=0 这个ans的位置错误
X(4) = B(4)/A(4,4)
X(3) = (B(3)-A(3,4 ...
欢迎光临 Fortran Coder (http://bbs.fcode.cn/) | Powered by Discuz! X3.2 |