function [t,y] = ukol2(t0,b,n)
% [t,y] = ukol2(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;
z0 = 1;
g = 10;
l = 1;

% doplnte kod
%
%
%
%
%
% doplnte kod 

%figure
%plot(t,y)

% konec