24小时热门版块排行榜    

CyRhmU.jpeg
查看: 325  |  回复: 0

springer_

木虫 (著名写手)

[求助] 上三角法解线性方程组,虽然没报错,但结果不正确。大神们帮忙看看

上三角法解线性方程组,虽然没报错,但结果不正确,大神们帮忙看看
program main
     implicit none
         integer,parameter :: n=4
         real :: A(n,n),b(n),x(n)
         A = reshape([2,1,3,4,0,3,2,5,0,0,7,3,0,0,0,2],[4,4])
         b = [50,49,53,12]
         
         call uptri(A,b,x,n)
         write(*,101) x
     101 format (T5,'计算结果',/,'x=',4(/F12.6))
end program main



subroutine uptri(A,b,x,n)
     implicit none
     integer :: i,j,n
     real :: A(n,n),b(n),x(n)
         
     x(n) = b(n) / A(n,n)
         do i = n-1,1,-1
             x(i) = b(i)
                 do j =i+1,n
                     x(i) = x(i) - A(i,j) * x(j)
                 end do       
                 x(i) = x(i) / A(i,i)
         end do         
end subroutine uptri
回复此楼
坐在墙头等红杏
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 springer_ 的主题更新
信息提示
请填处理意见