I try to solve rate equation using simulink in Matlab, more precisely,
the simulink functions ode23 and ode45 have been used to solve rate equation to verify chaotic dynamics of VCSEL? each model gives different behaviours and results,
in the first model, some results seem to be OK (which is not in the second model) but the rest not sure! (but seem to be OK in the second model!)