let u(t) be continuously differentiable numerical solution to the problem y'(t)=f(t,y) , y(t0)=y0 on the interval [0,100] for time. The defect of u(t) is defined as :

s(t)=u'(t)-f(t,u(t)).

Let s1 be the defect of u1 and s2 be the defect of u2 . write the matlab function to calculate u'. using ODE45.