% initial values
t0 = 0;
y0 = 1;
tend = 1;

% compute ivp solution
h = 0.01;
nsteps = tend/h;
[timepoints, yapprox] = blackbox_method(@rhs, t0, y0, nsteps, h);

exact_sol = y0*exp(tend);

fprintf('y(%f) = %f\n', tend, yapprox(end))

blackbox_error = abs(yapprox(end)-exact_sol)