3 void eraRz(
double psi,
double r[3][3])
35 double s, c, a00, a01, a02, a10, a11, a12;
40 a00 = c*r[0][0] + s*r[1][0];
41 a01 = c*r[0][1] + s*r[1][1];
42 a02 = c*r[0][2] + s*r[1][2];
43 a10 = - s*r[0][0] + c*r[1][0];
44 a11 = - s*r[0][1] + c*r[1][1];
45 a12 = - s*r[0][2] + c*r[1][2];
void eraRz(double psi, double r[3][3])