Default value for t_eval for solve_ivp

[KochdumperNiklas]

I think it would be nice that if t_eval is not explicitly specified by the user, it is automatically computed by the time-step size of the model and the size of the input vector. Currently an error message is outputted if t_eval is not provided, which is a little bit annoying in my opinion.