function [t,C] = koncentrace_reseni(tN,n,C0,Cl,v_d)
% [t,C] = koncentrace_reseni(tN,n,C0,Cl,v_d)
%
% Funkce vyresi rovnici C'(t) = -k*C(t) popisujici pokles koncentrace leku 
% v krvi s pocatecni podminkou C(0) = C0 na <0,tN> Eulerovou metodou
%
% tN  ... horni mez intervalu, na kterem resime nasi rovnici (v hodinach)
% n   ... pocet deleni
% C0  ... koncentrace leku v case 0
% Cl  ... clearance
% v_d ... distribucni objem
%
% t ... diskretizace <0,tN>
% C ... hodnoty hledane koncentrace

k = Cl*60/1000/v_d;

h = tN / n;
t = 0 : h : tN;

C(1) = C0;
for i = 1 : n
  C(i+1) = C(i) - h*k*C(i);
end

plot(t,C)

% konec
