Hello Massimiliano,
Please refer to example given in the following documentation:
a = jacobian([x*y*z, y^2, x + z], [x, y, z]);
subs(a,{x,y,z},{1,2,3});
Above code represents how to use 'subs' function to get the value at given point.
For more details please refer following documentation: