function diffeq_unidirectional

close all

tic

y0=0

options = odeset('RelTol',1e-4,'AbsTol',[1e-4],'stats','on','outputfcn','odeplot');

[t,y]=ode15s(@unidirectional_solve,[0 100],y0)

plot(t,y,'-o')

toc