求助Gauss消元法求解线性方程组的数值不稳定问题

Transpose

[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