function [t,y,cas] = priklad2()
% [t,y,cas] = priklad2()
%
% Funkce vyresi rovnici y(t)' = k(Tv - y(t)) s pocatecni podminkou y(t0)=y0 na <a,b>
% Eulerovou metodou
%
% a ... dolni mez intervalu
% b ... horni mez intervalu
% n ... pocet deleni
%
% t ... diskretizace <a,b>
% y ... hodnoty hledane funkce

y0=100;
k=0.04;

h=0.05;
t(1) = 0;

y(1) = y0;
i=1;

while y(i)>50
  y(i+1) = y(i) + h*(k*(20 - y(i)));
  t(i+1)=t(i)+h;
  i=i+1;
end

plot(t,y)
cas=t(i);
fprintf('Hrnek se ochladi na 50 stupnu Celsia za %f minut.\n',cas)

%konec
