program main
use imsl_libraries
implicit none
integer,parameter::n=3
real::a(n,n)=(/1.0,2.0,3.0,4.0,5.0,6.0,7.0,8.0,9.0/)
real::b(n,n)=(/5.0,6.0,4.0,8.0,9.0,7.0,2.0,3.0,1.0/)
real::c(n,n)
c=.i.a
write(*,*)c
end program
program main
include 'link_fnl_shared.h'
include 'link_fnl_static.h'
use operation_x !//这句
implicit none
integer,parameter::m=3,n=3
integer::mi,ni,iaq,jaq
real,parameter::hoa=7.0,h=0.1
real::a(m,n)
real::t(m,n)
real::e(m,n)
!
open(24,file='test12311339.txt')
do mi=1,m
do ni=1,n
t(mi,ni)=hoa*ni+hoa*mi
a(mi,ni)=hoa*(ni-1)+hoa*mi
end do
end do
e=a.x.t
write(24,*)e
end program