parameter( MAX = 100)
double precision A(MAX,MAX),B(MAX)
integer IPIV(MAX)
integer N
integer INFO
integer LWORK
integer NSINFO
...
call DGESV(N,1,A,MAX,IPIV,B,MAX,INFO)
...
...
call NETSL('DGESV',NSINFO,
$ N,1,A,MAX,IPIV,B,MAX,INFO)
...
...