In this paper, we derive an algorithm and develope a computer program which analyze rapidly and precisely the inverse kinematics of robotic mechanism with spatial complex chain structure based on the relative coordinates. We represent the inverse kinematic problem as an optimization problem with the kinematic constraint equations. The inverse kinematic analysis algorithm, therefore, consists of two algorithms, the main, an optimization algorithm finding the motion of independent joints from that of an end-effector and the sub, a forward kinematic analysis algorithm computing the motion of dependent joints. We accomplish simulations for the investigation upon the accuracy and efficiency of the algorithm.