dN = 150e-3; dX = 30e-3; A = 30;
Kov  =1; Tx = 42; Tr = 36;
Snom = 20;

Eps = 0.6;
nTakt = 10;
sigma = []; miu = []; ksi = []; vreme = []; lambda = [];
tau = zeros(1, nTakt);
tick = 0;
H = -Kov*A;
L=0;
G1 = Eps - dN/2+dX;
G2 = Eps - dN/2;
level = 0;
for i = 1:nTakt
    if (level == 0)
        level = 1;
            for t = 0:0.001:5
            epsilona = H*exp(-t/Tz)+Kov*A;
                vreme = [vreme t+sum(tau(1:i-1))]
                ksi = [ksi epsilona];
                sigma = [sigma Eps-epsilona];
                miut = Snom*(tick+t);
                miu = [miu miut];
                if (epsilona > G1)
                    tau(i) = t;
                    L = epsilona;
                    break;
                end
            end
            tick = tick + tau(i);
    else
        level = 0;
            for t = 0:0.001:5
            epsilona = L*exp(-t/Tr);
                vreme = [vreme t+sum(tau(1:i-1))];
                ksi = [ksi epsilona];
                sigma = [sigma Eps-epsilona];
                lambda = [lambda 0];
                miu = [miu miut];
                if(epsilona < G2)
                    tau(i) = t;
                    H = epsilona - Kov*A;
                    break;
                end
            end
    end
end