How can I find the partial derivative w.r.t angles

Good evening!

I have this:

And I want to know how can I find this with Matlab, w.r.t q1 and q2

in order to verify my results.

Thanks in advance!