function [t] = stari(N_tN,r)
% [t] = stari(N_tN,r)
%
% Funkce zjisti, zda je obraz se zadanymi parametry stary 20 ci 300 let.
% Pokles poctu atomu Pb-210 v 1 gramu oxidu olovnateho popisujeme
% rovnici N'(t) = -lambda*N(t) + r s pocatecni podminkou N(0) = N0 na <0,tN> 
%
% N_tN  ... pocet atomu Pb-210 v 1 gramu oxidu olovnateho v case tN
% r ... pocet atomu Pb-210, ktere vzniknou v 1 gramu oxidu olovnateho za 1 minutu rozpadem Ra-226 
%
% t ... zjistene stari obrazu (20 nebo 300 let)

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

tN = 300;
n = 10000; % pocet deleni <0,tN>

tN = tN * 525600;

%dodelat
%
%
%
%
%
%dodelat


disp(['Obraz je pravdepodobne stary 300 let a N(0) = ', num2str(N0,'%1.2g'),'.'])
return

%dodelat
%
%
%dodelat

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55

tN = 20;
n = 10000;

tN = tN*525600;

%dodelat
%
%
%
%
%
%dodelat

disp(['Obraz je pravdepodobne stary 20 let a N(0) = ', num2str(N0,'%1.2g'),'.'])
return
    
%dodelat
%
%
%dodelat

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%55

error(['Nejsme schopni rozhodnout, zda je obraz stary 300 nebo 20 let.'])

% konec
