Instead of solving the inverse kinematics analytically, we will use the differential kinematics approach to solve the inverse problem numerically. Here, the operational space is set as the cartesian tip position = [, ]and the joint angles are = [1 2 3]T . In the simulation, assume all the link lengths are set to 1