Numerical method for inverse kinematics using an extended angle-axis vector to avoid deadlock caused by joint limits*. (3rd August 2021)