1mat : genmatrix(mat, 3,3, 1,1)$
2
3mat[1,1] :
4     18*cos(q3)*cos(q2)*m30*p^2 - 9*sin(q3)^2*p^2*m30
5     - sin(q3)^2*j30y + sin(q3)^2*j30z + p^2*m10 + 18*p^2*m30
6     + j10y + j30y$
7
8mat[2,1] :
9     mat[1,2] :
10     9*cos(q3)*cos(q2)*m30*p^2 - sin(q3)^2*j30y + sin(q3)^2*j30z
11     - 9*sin(q3)^2*m30*p^2 + j30y + 9*m30*p^2$
12
13mat[3,1] :
14     mat[1,3] :
15     -9*sin(q3)*sin(q2)*m30*p^2$
16mat[2,2] :
17     -sin(q3)^2*j30y + sin(q3)^2*j30z - 9*sin(q3)^2 *m30*p^2
18     + j30y + 9*m30*p^2$
19
20mat[3,2] : mat[2,3] : 0$
21
22mat[3,3] : 9*m30*p^2 + j30x$
23
24gentranout("array.f")$
25for i:1 thru 3 do
26	  for j:i thru 3 do
27	      gentran( lrsetq(mat[i,j], mat[i,j]) )$
28
29gentranshut("array.f")$
30