I have been working on a small Python package to solve a class of PDEs using scipy.integrate.solve_ivp.