%% Řešení ODR pomocí ode45

f=@(x,y)-y.^2;
I=[1,5];
y0=1;
[xx,yy]=ode45(f,I,y0);
plot(xx,yy,xx,1./xx,'--')
%xx'         % vypis uzlu
diff(xx')   % vzdalenosti uzlu

%
% Zobrazení chyby

plot(xx,yy-1./xx)

%
% Výpočet pro danou síť uzlů

x1=linspace(1,5,9);
[xx,yy]=ode45(f,x1,y0);
plot(xx,yy,xx,1./xx,'--')
plot(xx,yy-1./xx)

%
% Více počátečních podmínek

y0=[1 2 3];
[xx,yy]=ode45(f,I,y0);
plot(xx,yy)