clear;
clf;

t = 0;
y = 1;
tstop = 2;

tspan = [t tstop];
y0 = y;

[tt,yy] = ode23('func',tspan,y0);


% plot numerical solution
plot(tt,yy, 'ro');
axis([min(tt)-1 max(tt)+1 min(yy)-1 max(yy)+1]);
xlabel('tt');
ylabel('yy');

% plot exact solution
hold on;
ttt = tt(1) : 0.01 : tt(length(tt));
yyy = exct(ttt);
plot(ttt,yyy);