[comp.lang.fortran] IMSL

muttiah@fescue.ecn.purdue.edu (Ranjan S Muttiah) (02/13/91)

main()
{
for(j=0;j < N;j++)
	hammertime();
}

hammertime()
{
float X[3][25]; ...

ido = 0; ncol=3;
nef=1; nvef[0] = 2; indef[0] = 2; indef[1] = 3;
inddep[0] = 1;

rglm_(.....);	/* imsl call */

}

Note: I am calling IMSL from within C.

------------
I get an error about indef[2] = 3 being greater than ncol=2 (at runtime).  
Unbelievable! I just assigned ncol=3.  When ido = 1, it works but wrong
regression coefficients.