Very simple question - can I use a non symbolic version of the rhs function and if so, how?
Specifically, a case of ode(t ,y, p) where t, y and p are lists or numpy arrays, and the function returns dYdt (a list or numpy array where each element is a different state equation)