HOW TO CALL THE IMSL FORTRAN LIBRARIES FROM C by ...

3.4 Example
/*
This example solves a linear system of 100 by 100
using the IMSL routine LSLRG.
*/
#include
#include
#include
#include
#define MAXN 100
main()
{
int
long int
long int
double
double
i, nerr;
n=MAXN, nn=MAXN*MAXN, lda=MAXN;
ipath=2, seed=123457;
ermax=l.e-7;
a[MAXN][MAXN], b[MAXN], x[MAXN], ans[MAXN];
/* generate matrix a with uniform random numbers */
RNSET (&seed);
RNUN (&nn, a);
/* generate a pre-determined answer, ans */
for(i=0;i