I can do the plot in MATLAB. Just solving the equation is enough.
Design a Cartesian space trajectory planning system for a planar 2-link RR (link lengths are 1 meter each) robot that moves the tool from an initial position (-0.1, 0.8) to a final position (-0.1, -0.8) with respect to the base frame B, which is assumed to be at (0,0) along a straight line in 5 seconds. Set up all the equations properly. Determine and plot (using MATLAB) the joint space trajectory and velocity profile (using J^-1 for joint rates and INVKIN for joint values at an appropriate sampling rate) as a function of time. What happens to joint rates near a singularity? Explain why it happens. Suggest a possible mechanism to deal with it. Please refer to Equation 5.74 in your text for J^-i and to Equations 4.16 and 4.27 for INVKIN.