How can I include the orientation in the Jacobian matrix for differential kinematics?

Question

I\'m writing a position controller for a 7 DOF manipulator (Kuka LWR) using the pseudo-inverse Jacobian. I already have a Matlab function that generates a 3x7 matrix that ac