create.jacob function

Internal function to compute Jacobian