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

