I am trying to make a function that can solve general ODEs in python. My problem has nothing to do with ODEs themselves so to make it easier to understand the problem i have simplified it. I want the user of the function to define ODEs as functions which takes in time(scalar value), value(numpy array) and some arguments specific to the ODE. Consider the examples below: