I need to do numerical simulations such as phase plots, bifurcation diagrams, parametric basins of attraction, and Lyapunov exponent diagrams. I already have the fixed points though.
With regard to Numerical Simulation, I have something to offer. But, what type of Differential Equations are included in the model ? is it ODE, stiff ODE or PDE and of what order?, what is the form of the DE? is it delay , fussy, Wave/Heat etc?. Let see the model [email protected]
Since this is a quasilinear system of ODEs, you could simply use an implicit Euler scheme so you end up with an algebraic system of two equations and two unknowns that can be easily solved with Newton's method.