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 ...

感谢老师指正，我搞错了ans赋值语句的位置，现在程序运行正常了。
还想麻烦老师帮忙看看这个程序在数值算法和语法规范等方面有没有可以优化的地方，非常感谢！

[Fortran] 查看源码 复制源码

   !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