%Class 27 Warm Up
% Analytical method using laplace transform
yfunc = @(t) t^2/2-exp(-2*t)/4-t/2+1/4+exp(-2*t);
fplot(yfunc,[0 5])
hold on
% Use ode45 for the RK method in this case:
dydt = @(t,y) -2*y+t^2;
[dy,dt] = ode45(dydt,[0 5],1);
plot(dy,dt,'o')
%
xlabel('t')
ylabel('y')