We build our best fitting model with parameters a and b.

Input := 

u[t_] = f[t]/.sol[[2]]

Output =

       1
----------------
      53.9126
1 + ------------
     0.0401633 t
    E