%
% odetest.m
%
% Test fuer ode23, ode45;
% verwendet Funktion ftest
t0 = 0;     % Startwert f"ur die Zeit 
y0 = [0  1]; % Startwert f"ur y0
tfinal = 40.*pi;
tspan=[t0 tfinal]; 
%[t,y] = ode23(@ftest,tspan,y0);
%
opts = odeset('RelTol',1e-5);
[t,y] = ode23(@ftest,tspan,y0,opts);
%
%[t,y] = ode45(@ftest,tspan,y0);

dis = [sin(t) cos(t)]-y; % Diskrepanz zur exakten Loesung
figure(1);
plot(t,dis);
title('Diskrepanz zur Loesung')
figure(2)
plot(y(:,1),y(:,2));
axis square   
title('Loesung, Komponeten gegeneinander'); 
