If you are using the symbolic toolbox, then MATLAB has jacobian()
Otherwise, you need to estimate the jacobian, which you can do by calculating the gradient() with respect to each variable. You can do that by iterating through the variables, leaving all the parameters except that one fixed.