lab9_2

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('Comparsion of the decisiov at various errors')
legend('10^{-3}','10^{-4}','10^{-6}','The exact decision',4)
grid on
hold off