When i try to get all roots for an ODE, i try to use fsolve:
fsolve((rhs(ODE_4))(t)=0,t=10..20);
where ODE_4 is an ODE (dsolve, numeric, output=listprocedure), where i want to find the root of the 4th variable in the ODE (4 should be an index).
Then i only get the fist root after t=10 and the computation stops. What can i do to get ALL roots in the given interval?
