doublec1=sin(phi)*cos(kappa)+cos(phi)*sin(omega)*sin(kappa);doublec2=-sin(phi)*sin(kappa)+cos(phi)*sin(omega)*cos(kappa);doublec3=cos(phi)*cos(omega);