Use solve in lti_disc
Now the model for each dt is calculated only once