share/contrib/gentran/test/array.mac

   1 mat : genmatrix(mat, 3,3, 1,1)$
   2
   3 mat[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
   8 mat[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
  13 mat[3,1] :
  14      mat[1,3] :
  15      -9*sin(q3)*sin(q2)*m30*p^2$
  16 mat[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
  20 mat[3,2] : mat[2,3] : 0$
  21
  22 mat[3,3] : 9*m30*p^2 + j30x$
  23
  24 gentranout("array.f")$
  25 for i:1 thru 3 do
  26           for j:i thru 3 do
  27               gentran( lrsetq(mat[i,j], mat[i,j]) )$
  28
  29 gentranshut("array.f")$