Sorry,

I actually had a mistake in the definition of the matrix. I spent all day long... well, it's a good way to learn

The code is:

Then, I compute the Jacobian with the package VectorCalculus:

What I optain is a very complex matrix. Is there any way to have the matrix more "readable" ?