function [t,y] = ukol1_reseni(t0,b,n)
% [t,y] = ukol1_reseni(t0,b,n)
%
% funkce vyresi priblizne soustavu rovnic y'(t) = z(t), z'(t) = -(g/l)y(t) 
% s pocatecnimi podminkami y(t0) = y0, z(t0) = z0 na <t0,b> Eulerovou metodou
% (tj. aproximuje vychylku daneho kyvadla pomoci Eulerovy metody)
%
% t0 ... dolni mez intervalu
% b ... horni mez intervalu
% n ... pocet deleni intervalu <t0,b>
%
% t ... diskretizace <t0,b>
% y ... hodnoty hledane funkce

% parametry kyvadla
y0 = 0.1;
z0 = 0;
g = 10;
l = 1;

h = (b - t0) / n; % krok
t = t0 : h : b;

y(1) = y0;
z(1) = z0;
for i = 1 : n
  y(i+1) = y(i) + h * z(i);
  z(i+1) = z(i) + h * (-g/l)*y(i);
end

figure
plot(t,y)

% konec
