A0[0][0] = 1.0/(q0*q0-q1*q1-q2*q2+q3*q3);