y0 = [ 5 -4.4 ]';
tspan = 0:0.1:10;
[ t y ] = ode45(@diffeqs,tspan,y0);
plot(t,y(:,1),'k-'); grid
xlabel('t')
ylabel('y')

