title: Solve ax+by+cz = k using numpy
body: <p>I have scenario where I need to solve following equation and have some data.
  We should have only integer value. </p>
  
  <p>ax+by = c
      The value of which we have a,b and c and need to determine x and y.</p>
  
  <pre><code>396x+264y = 1980
  Ans :The value of x=5 and y=0
  
  396x+264y = 2508
  Ans : The vaule of x=5 and y=2
  
  In future, we may add dz.
  396+264y+dz = c(some value)
  </code></pre>
  
  <p>I understand, it can contain multiple solution but I would like to how can I determine those solution. In my scenario is I should be get x, and y integer positive value </p>
  
  <p>is it possible to solve this equation using numpy? I thought, it's like gcd
  factor.</p>