A quick idea.
Since you make assumptions when creating a differential equations by hand. I assume this can be generalised by using machine learning models.
So I guess its possible to simulate a slope field. The slope field is the equation. I picture a grid with randomly rotating vectors that the model then can iterate each a better and better value for.
This way your not bound by symbol math. Only by the data?