matrix.
get w  /file = * /variables = q1 q2 q3 q4 q5 q6.
get z /file = * /variables = p1 p2 p3 p4 p5 p6.
 
* print w.
* print z.
 
compute a = (t(w)*w).
compute b = (t(z)*z).
 
compute a_g =ginv(a).
compute b_g = ginv(b).
 
compute delta = a_g*t(w)*z.
 
compute gamma = b_g*t(z)*w.
 
* this the new betas in term of u1 u2 u3 b1 b2 b3.
 
print delta.
 
 
* this is u1 u2 u3 b1 b2 b3 in terms of u a1 a2 b d1 d2.
 
print gamma.
 
compute gamma_i = ginv(gamma).
 
print gamma_i.
 
compute  l =gamma_i*gamma.
print l.
end matrix.