/*
** Example 1.1
*/
one=ones(1,4);
e=1.0;
do until e < 1.0e-20;
    x=one|(e.*eye(4));
    print "e = " e;
    r=eigrs(x'x); print r';
    invx=invpd(x'x); print invx;
    e=e/10;
endo;

end;