Added function to calculate quaternion axes