th 1[m] "Thickness of circular sheet"
Rin 5[m] "Inner radius of circular sheet"
Rout 12.5[m] "Outer radius of circular sheet"
theta 0[deg] Rotation
rxx 0 "Intermediate help variable"
ryy 0 "Intermediate help variable"
rzz 1 "Intermediate help variable"
rx rxx/sqrt(rxx^2+ryy^2+rzz^2) "Intermediate help variable"
ry ryy/sqrt(rxx^2+ryy^2+rzz^2) "Intermediate help variable"
rz rzz/sqrt(rxx^2+ryy^2+rzz^2) "Intermediate help variable"
ct cos(theta) "Intermediate help variable"
st sin(theta) "Intermediate help variable"
R11 ct+rx^2*(1-ct) "Rotation matrix, 11 component"
R12 rx*ry*(1-ct)-rz*st "Rotation matrix, 12 component"
R13 rx*rz*(1-ct)+ry*st "Rotation matrix, 13 component"
R21 rx*ry*(1-ct)+rz*st "Rotation matrix, 21 component"
R22 ct+ry^2*(1-ct) "Rotation matrix, 22 component"
R23 rz*ry*(1-ct)-rx*st "Rotation matrix, 23 component"
R31 rx*rz*(1-ct)-ry*st "Rotation matrix, 31 component"
R32 rz*ry*(1-ct)+rx*st "Rotation matrix, 32 component"
R33 ct+rz^2*(1-ct) "Rotation matrix, 33 component"