title: Is there a way to evaluate the Jacobian in Matlab?
body: <p>I have this code:</p>
  
  <pre><code>syms x y z
  f =[x^2+y^2+z^2-100 ,x*y*z-1 ,x-y-sin(z)];
  v = [x, y, z];
  fp = jacobian(f,v)
  </code></pre>
  
  <p>This gives:</p>
  
  <pre><code>fp =
  [ 2*x, 2*y,     2*z]
  [ y*z, x*z,     x*y]
  [   1,  -1, -cos(z)]
  </code></pre>
  
  <p>Now I want to evaluate this at, say, x=y=z=1. But</p>
  
  <pre><code>fp([1,1,1])
  </code></pre>
  
  <p>gives me</p>
  
  <pre><code>[ 2*x, 2*x, 2*x]
  </code></pre>
  
  <p>Is there a way to do this?</p>