
[t, y] = ode45(@solutions_Final_p1a, [0 1], 0);
fprintf('Solution at t = 1: %.5f\n', y(end));
fprintf('   Error at t = 1: %.5e\n', abs(y(end) - exp(-1)));