Using simult_ with third order approximation

Please take a look at [Choosing initial values for simult_)