求一个fortran解决线性方程的代码.

来源：学生作业帮编辑：作业帮分类：综合作业时间：2024/05/23 22:13:02

求一个fortran解决线性方程的代码.
补充：用高斯消元法,要求有例子的.

我运行通过了~~
program main
implicit none
integer :: N = 3

real*8 :: A(3,3) = (/ 2, 3, 1, 1, 1, 2, 1, 2, 2 /)
real*8 :: B(3) = (/ 4, 6, 5/), X(3)

call LinearGauss (A, B, X, N)

write(*,*) X

end
!-------------------------------------------------
! 行尺度高斯消元法
!-------------------------------------------------
subroutine LinearGauss (A, B, X, N)
implicit none
integer :: N
real*8 :: A(N, N), B(N), X(N)

real*8 :: z
integer,allocatable :: p(:), s(:)
integer :: i, j, k, sl

allocate (p(N))
allocate (s(N))

do i=1, N
p(i) = i
s(i) = MaxVal ( Abs(A(i, :)) )
end do

do k=1, N-1

j = k
do i=k+1, N
if ( Abs(A(p(j), k)) / s(p(j)) < Abs(A(p(i), k)) / s(p(i)) ) then
j = i
end if
end do
sl = p(k)
p(k) = p(j)
p(j) = sl

do i=k+1, N
z = A(p(i), k) / A(p(k), k)
A(p(i), k) = z
do j=k+1, N
A(p(i), j) = A(p(i), j) - z * A(p(k), j)
end do
end do

end do

do k=1, N-1
do i=k+1, N
B(p(i)) = B(p(i)) - A(p(i), k) * b(p(k))
end do
end do

do i=N, 1, -1
z = 0
do j=i+1, N
z = z + A(p(i), j) * X(j)
end do
X(i) = (B(p(i)) - z) / A(p(i), i)
end do

end subroutine

急需用matlab解线性方程的代码! 一个公式写成Fortran语言代码 LINGO 模型的代码求高手版解决代码错误 fortran的一个问题求大神解答求一个Fortran程序求fortran高手帮我看一下这段用高斯乔丹消元法解方程组的代码为什么总是得不到想要的结果一个简单的FORTRAN程序, 一个fortran语句的意思求fortran大神帮我编一个Fortran程序计算无理数π的1-100的小数位, 矩阵,线性方程的一个简单题目线性代数.求下列齐次线性方程的通解和一个基础解析急求lingo代码解决极值问题