Switch from internal ODE solver to scipy

Created by: gcarleo

At this point we should move from our home-made implementation of Runge Kutta, and just rely on scipy.integrate.ode. Using to_sparse, to_linear_operator etc, it should be really straightforward to rewrite the exact real/imaginary time evolution in a pure python function!