program Console2
include 'link_fnl_shared.h'
use imsl_libraries
integer,parameter::m=5
integer,parameter::n=2
integer::i
double precision::a(m,n),b(m),x(n)
double precision::c(n,m),d(n,n),e(n)
data a/2,8,2,7,4,-1,4,1,-1,0/
data b/1,0,1,8,3/
if(m.gt.n)then
c=transpose(a)
d=matmul(c,a)
e=matmul(c,b)
x=d.ix.e
else if(m.eq.n)then
x=a.ix.b
else
write(*,*)'no roots!'
end if
do i=1,n
write(*,*)x(i)
end do
end program Console2