/* timing gausswin vs 3.1.5*/
tim = hsec;
seed1 = 7654321;
i = 1;
do until i > 10;

x = rndns(100,100,seed1);
xi = inv(x);
errx = maxc(abs(vecr(x*xi - eye(100) )));
"error" errx;

i = i + 1;
endo;
tim2 = hsec - tim;

"time" 0.01*tim2;