Use solve_ivp without the py-pde wrapper, but retain fields and grids from py-pde, for more readable code

Use solve_ivp without the py-pde wrapper, but retain fields and grids from py-pde, for more readable code #76