matrix.
get x/variables=all/file=*.
compute a = x*t(x).
print a.
compute lambda =0.
compute nr = nrow(a).
call eigen(a, eigenvec,eigenval).
print eigenval.
*print eigenvec.
compute id=make(nr,nr,0).
*print id.
loop i = 1 to nr.
compute id(i,i)=1.
end loop.
*print id.
loop j = 1 to 4.
compute lambda(1,1)=eigenval(j,1).
print lambda.
compute k =(id-lambda*a).
compute f = ginv(k).
compute w=k*f.
compute tr = trace(w).
print tr.
print w.
end loop.
end matrix.