title: Compute the Jacobian matrix in Python
body: <pre><code>import numpy as np
  
  
  a = np.array([[1,2,3],
                [4,5,6],
                [7,8,9]])
  
  
  b = np.array([[1,2,3]]).T
  
  c = a.dot(b) #function
  
  jacobian = a # as partial derivative of c w.r.t to b is a.
  </code></pre>
  
  <p>I am reading about jacobian Matrix, trying to build one and from what I have read so far, this python code should be considered as jacobian. Am I understanding this right?</p>