%% t = 0.4

[a, b, c] = euler_explicit(@(x, t) 3*exp(-201*t) - 5*exp(-t), @(x, t) -201*x - 100*exp(-t), -2, 0, 1, 0.1);
% [a, b, c] = euler_explicit(@(x, y) -0.5*x^4 + 4*x^3 - 10*x^2 + 8.5*x + 1, @(x, y) -2*x^3 + 12*x^2 - 20*x + 8.5, 1, 0, 4, 0.005);

disp(['The Explicit Euler estimation at ', num2str(10)]);
disp(['is: ', num2str(a)]);
disp(['The actual value is ', num2str(b)]);
disp(['The relative percent of error is ', num2str(c)]);

