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