function ferr = solveODEs(y10)
y0 = [ 5 y10 ]';
tspan = [0 10];
[ t y ] = ode45(@diffeqs,tspan,y0);
ferr = y(end,1) - 8;

