add a homotopy mode for simul() and for steady with endval
Currently homotopy only exist for steady with initval.
These features would be useful for GIMF at the IMF.
Currently homotopy only exist for steady with initval.
These features would be useful for GIMF at the IMF.