Лаба 9 / Lab9

a = 0.001;
Y0 = [log(a); 1/a];
options = odeset('RelTol', 1.0e-3);
[T, Y] = ode45('diff2', [a 100], Y0, options);
plot(T, Y(:, 1), 'k:')
options = odeset('RelTol', 1.0e-4);
[T, Y] = ode45('diff2', [a 100], Y0, options);
hold on
plot(T, Y(:, 1), 'k--')
options = odeset('RelTol', 1.0e-6);
[T, Y] = ode45('diff2', [a 100], Y0, options);
hold on
plot(T, Y(:, 1), 'k-')
t = [a:4:100];
y = log(t) ;
hold on
plot(t, y, 'ko')
xlabel('\itt')
ylabel('\ity')
title('Comparison of the decision at various errors')
legend('10^{-3}', '10^{-4}', '10^{-6}', 'The exact decision')
grid on
hold off