Generated program
c
c *** compute values for matrix m ***
c
m(1,1)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+
. sin(real(q3))**2*j30z+18.0*cos(real(q3))*cos(real(q2))*p**2*m30+
. 18.0*p**2*m30+p**2*m10+j30y+j10y
m(1,2)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+
. sin(real(q3))**2*j30z+9.0*cos(real(q3))*cos(real(q2))*p**2*m30+
. 9.0*p**2*m30+j30y
m(1,3)=-(9.0*sin(real(q3))*sin(real(q2))*p**2*m30)
m(2,2)=-(9.0*sin(real(q3))**2*p**2*m30)-(sin(real(q3))**2*j30y)+
. sin(real(q3))**2*j30z+9.0*p**2*m30+j30y
m(2,3)=0.0
m(3,3)=9.0*p**2*m30+j30x
c
c *** compute values for inverse matrix ***
c
t0=m(1,1)
t1=m(1,2)
t2=m(1,3)
t3=m(2,2)
t4=m(3,3)
mxinv(1,1)=-(t3*t4)/(t1**2*t4+t2**2*t3-(t0*t3*t4))
mxinv(1,2)=(t1*t4)/(t1**2*t4+t2**2*t3-(t0*t3*t4))
mxinv(1,3)=(t2*t3)/(t1**2*t4+t2**2*t3-(t0*t3*t4))
mxinv(2,2)=(t2**2-(t0*t4))/(t1**2*t4+t2**2*t3-(t0*t3*t4))
mxinv(2,3)=-(t1*t2)/(t1**2*t4+t2**2*t3-(t0*t3*t4))
mxinv(3,3)=(t1**2-(t0*t3))/(t1**2*t4+t2**2*t3-(t0*t3*t4))
do 25001 j=1,3
do 25002 k=j+1,3
m(k,j)=m(j,k)
mxinv(k,j)=mxinv(j,k)
25002 continue
25001 continue