3 void eraRxr(
double a[3][3],
double b[3][3],
double atb[3][3])
32 for (i = 0; i < 3; i++) {
33 for (j = 0; j < 3; j++) {
35 for (k = 0; k < 3; k++) {
36 w += a[
i][k] * b[k][j];
void eraRxr(double a[3][3], double b[3][3], double atb[3][3])
void eraCr(double r[3][3], double c[3][3])