linear.f90

来自「Fortran 95程序设计【彭国伦】 Source_Code.rarFortr」· F90 代码 · 共 16 行

program main
  use IMSL
  implicit none
  real :: A(3,3) = (/ 1,3,2,&
                      1,2,1,&
					  2,1,3 /)
  real :: B(3,1) = (/ 4,6,6 /)
  real :: X(3,1)

  call lin_sol_gen(A,B,X) ! A*X=B,解X
  write(*,"(3F5.2)") X

  stop
end program