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.