def rotation(axis_x,axis_y,axis_z,angle): s = sin(angle/2) return Quarternion(cos(angle/2),s*axis_x,s*axis_y,s*axis_z)