title: Scipy - solving ODE system of the form: A*dy/dt = f(t, y)
body: <p>What are my options to solve in scipy the ODE system:</p>
  
  <p><img src="https://latex.codecogs.com/gif.latex?A%5Cfrac%7Bdy%7D%7Bdt%7D&space;=&space;f(t,y)" title="A\frac{dy}{dt} = f(y,t)" /></p>
  
  <p>Obviously if A is invertibile then i can solve equivalently:</p>
  
  <p><img src="https://latex.codecogs.com/gif.latex?%5Cfrac%7Bdy%7D%7Bdt%7D&space;=&space;A%5E%7B-1%7D&space;f(t,%20y)" title="\frac{dy}{dt} = A^{-1} f(t, y)" /></p>
  
  <p>using standard methods from scipy.integrate. But still remain the following questions: what about efficiency and what if A is singular and cannot be inverted?</p>