function [t,N] = atomy(tN,n,N0,r)
% [t,N] = atomy(tN,n,N0,r)
%
% Funkce vyresi rovnici N'(t) = -lambda*N(t) + r popisujici pokles poctu
% atomu Pb-210 s pocatecni podminkou N(0)=N0 na <0,tN> v 1 gramu oxidu
% olovnateho Eulerovou metodou.
%
% tN  ... horni mez intervalu, na kterem resime nasi rovnici (v letech)
% n   ... pocet deleni
% N0  ... pocet atomu Pb-210 v 1 gramu oxidu olovnateho v case 0
% r ... pocet atomu Pb-210, ktere vzniknou v 1 gramu oxidu olovnateho za 1 minutu rozpadem Ra-226 
%
% t ... diskretizace <0,tN>
% N ... pocty atomu Pb-210 v 1 gramu oxidu olovnateho v casech t

p_r = 22; % polocas rozpadu (v letech)
lambda = log(2)/(p_r*525600);

tN = tN * 525600;

%dodelat
%
%
%
%
%
%dodelat

plot(t,N)

% konec
