Calculating Jacobian without DH Parameters

This is not exactly a development contribution in terms of code, it's more like an idea. Nevertheless, I wanted to note that I was helpful to someone doing some development by leading them to the right direction on the ode-users mailing list.