function [T_off T_error] = find_T_off(Tecg, maxT, maxT2, maxT2_rel,zaporne_T)

T_error=0;

%Computations of gradients
gradient2 = zeros(1,length(Tecg)-10); gradient4 = zeros(1,length(Tecg)-10); gradient6 = zeros(1,length(Tecg)-10);
gradient8 = zeros(1,length(Tecg)-10); gradient10 = zeros(1,length(Tecg)-10);
for i = 1:(length(Tecg)-21),
    gradient2(i) = (Tecg(i+2) - Tecg(i));
    gradient4(i) = (Tecg(i+4) - Tecg(i));
    gradient6(i) = (Tecg(i+6) - Tecg(i));
    gradient8(i) = (Tecg(i+8) - Tecg(i));
    gradient10(i) = (Tecg(i+10) - Tecg(i));
end

[unused1 posNejvetsi_gradient2] = max(gradient2);
[unused2 posNejvetsi_gradient4] = max(gradient4);
[unused3 posNejvetsi_gradient6] = max(gradient6);
[unused4 posNejvetsi_gradient8] = max(gradient8);
[unused5 posNejvetsi_gradient10] = max(gradient10);

pocet_minim_Ton2 = 1; pocet_minim_Ton4 = 1; pocet_minim_Ton6 = 1; pocet_minim_Ton8 = 1; pocet_minim_Ton10 = 1;
for i = 2:length(gradient2),
    if zaporne_T==0,
        if (((gradient2(i-1)<0)&&(gradient2(i+1)>0))||((gradient2(i-1)<0)&&(abs(gradient2(i+1))<abs(mean(gradient2(1:11))/5))))
            posT_on_minima2(pocet_minim_Ton2) = i;
            pocet_minim_Ton2 = pocet_minim_Ton2 + 1;
        end
        if ((gradient4(i-1)<0)&&(gradient4(i+1)>0))||((gradient2(i-1)<0)&&(abs(gradient4(i+1))<abs(mean(gradient4(1:11))/5)))
            posT_on_minima4(pocet_minim_Ton4)= i;
            pocet_minim_Ton4 = pocet_minim_Ton4 + 1;
        end
        if ((gradient6(i-1)<0)&&(gradient6(i)>0))||((gradient2(i-1)<0)&&(abs(gradient6(i+1))<abs(mean(gradient6(1:11))/5)))
            posT_on_minima6(pocet_minim_Ton6) = i;
            pocet_minim_Ton6 = pocet_minim_Ton6 + 1;
        end
        if ((gradient8(i-1)<0)&&(gradient8(i)>0))||((gradient2(i-1)<0)&&(abs(gradient8(i+1))<abs(mean(gradient8(1:11))/5)))
            posT_on_minima8(pocet_minim_Ton8) = i;
            pocet_minim_Ton8 = pocet_minim_Ton8 + 1;
        end
        if ((gradient10(i-1)<0)&&(gradient10(i)>0))||((gradient2(i-1)<0)&&(abs(gradient10(i+1))<abs(mean(gradient10(1:11))/5)))
            posT_on_minima10(pocet_minim_Ton10) = i;
            pocet_minim_Ton10 = pocet_minim_Ton10 + 1;
        end
    else
        if (i<maxT2_rel)||(maxT==maxT2),
            if ((gradient2(i-1)<0)&&(gradient2(i+1)>0))||((gradient2(i-1)>0)&&(abs(gradient2(i+1))<abs(mean(gradient2(1:11))/4)))
                posT_on_minima2(pocet_minim_Ton2) = i;
                pocet_minim_Ton2 = pocet_minim_Ton2 + 1;
            end
            if ((gradient4(i-1)<0)&&(gradient4(i+1)>0))||((gradient2(i-1)>0)&&(abs(gradient4(i+1))<abs(mean(gradient4(1:11))/4)))
                posT_on_minima4(pocet_minim_Ton4)= i;
                pocet_minim_Ton4 = pocet_minim_Ton4 + 1;
            end
            if ((gradient6(i-1)<0)&&(gradient6(i)>0))||((gradient2(i-1)>0)&&(abs(gradient6(i+1))<abs(mean(gradient6(1:11))/4)))
                posT_on_minima6(pocet_minim_Ton6) = i;
                pocet_minim_Ton6 = pocet_minim_Ton6 + 1;
            end
            if ((gradient8(i-1)<0)&&(gradient8(i)>0))||((gradient2(i-1)>0)&&(abs(gradient8(i+1))<abs(mean(gradient8(1:11))/4)))
                posT_on_minima8(pocet_minim_Ton8) = i;
                pocet_minim_Ton8 = pocet_minim_Ton8 + 1;
            end
            if ((gradient10(i-1)<0)&&(gradient10(i)>0))||((gradient2(i-1)>0)&&(abs(gradient10(i+1))<abs(mean(gradient10(1:11))/4)))
                posT_on_minima10(pocet_minim_Ton10) = i;
                pocet_minim_Ton10 = pocet_minim_Ton10 + 1;
            end
        else
            if (((gradient2(i-1)<0)&&(gradient2(i+1)>0))||((gradient2(i-1)<0)&&(abs(gradient2(i+1))<abs(mean(gradient2(1:11))/5))))
                posT_on_minima2(pocet_minim_Ton2) = i;
                pocet_minim_Ton2 = pocet_minim_Ton2 + 1;
            end
            if ((gradient4(i-1)<0)&&(gradient4(i+1)>0))||((gradient2(i-1)<0)&&(abs(gradient4(i+1))<abs(mean(gradient4(1:11))/5)))
                posT_on_minima4(pocet_minim_Ton4)= i;
                pocet_minim_Ton4 = pocet_minim_Ton4 + 1;
            end
            if ((gradient6(i-1)<0)&&(gradient6(i)>0))||((gradient2(i-1)<0)&&(abs(gradient6(i+1))<abs(mean(gradient6(1:11))/5)))
                posT_on_minima6(pocet_minim_Ton6) = i;
                pocet_minim_Ton6 = pocet_minim_Ton6 + 1;
            end
            if ((gradient8(i-1)<0)&&(gradient8(i)>0))||((gradient2(i-1)<0)&&(abs(gradient8(i+1))<abs(mean(gradient8(1:11))/5)))
                posT_on_minima8(pocet_minim_Ton8) = i;
                pocet_minim_Ton8 = pocet_minim_Ton8 + 1;
            end
            if ((gradient10(i-1)<0)&&(gradient10(i)>0))||((gradient2(i-1)<0)&&(abs(gradient10(i+1))<abs(mean(gradient10(1:11))/5)))
                posT_on_minima10(pocet_minim_Ton10) = i;
                pocet_minim_Ton10 = pocet_minim_Ton10 + 1;
            end
        end
    end
end

if ~exist('posT_on_minima2','var'),
    posT_on_minima2 = 1;
end
if ~exist('posT_on_minima4','var'),
    posT_on_minima4 = 1;
end
if ~exist('posT_on_minima6','var'),
    posT_on_minima6 = 1;
end
if ~exist('posT_on_minima8','var'),
    posT_on_minima8 = 1;
end
if ~exist('posT_on_minima10','var'),
    posT_on_minima10 = 1;
end

[pocet_minim_Ton nejvice_minim_z_gradientu] = max([length(posT_on_minima2),length(posT_on_minima4),length(posT_on_minima6),length(posT_on_minima8),length(posT_on_minima10)]);
offset_min = 4; byl_jsem_tu = 0; global_min = 0;
if pocet_minim_Ton>1,
    switch nejvice_minim_z_gradientu,
        case 1,
            for i = 1:length(posT_on_minima2),
                if (posT_on_minima2(i) > maxT2_rel+10),
                    fired_at_least = 0;
                    if pocet_minim_Ton4>1,
                        if (find((posT_on_minima4<posT_on_minima2(i)+offset_min)&(posT_on_minima4>posT_on_minima2(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton6>1,
                        if (find((posT_on_minima6<posT_on_minima2(i)+offset_min)&(posT_on_minima6>posT_on_minima2(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton8>1,
                        if (find((posT_on_minima8<posT_on_minima2(i)+offset_min)&(posT_on_minima8>posT_on_minima2(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton10>1,
                        if (find((posT_on_minima10<posT_on_minima2(i)+offset_min)&(posT_on_minima10>posT_on_minima2(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end

                    if ((fired_at_least > 3)&&(byl_jsem_tu == 0)),
                        T_off = posT_on_minima2(i);
                        byl_jsem_tu = 1;
                    end
                end
            end
        case 2,
            for i = 1:length(posT_on_minima4),
                if (posT_on_minima4(i) > maxT2_rel+10),
                    fired_at_least = 0;
                    if pocet_minim_Ton2>1,
                        if (find((posT_on_minima2<posT_on_minima4(i)+offset_min)&(posT_on_minima2>posT_on_minima4(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton6>1,
                        if (find((posT_on_minima6<posT_on_minima4(i)+offset_min)&(posT_on_minima6>posT_on_minima4(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton8>1,
                        if (find((posT_on_minima8<posT_on_minima4(i)+offset_min)&(posT_on_minima8>posT_on_minima4(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton10>1,
                        if (find((posT_on_minima10<posT_on_minima4(i)+offset_min)&(posT_on_minima10>posT_on_minima4(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if ((fired_at_least > 3)&&(byl_jsem_tu == 0)),
                        T_off = posT_on_minima4(i);
                        byl_jsem_tu = 1;
                    end
                end
            end
        case 3,
            for i = 1:length(posT_on_minima6),
                if (posT_on_minima6(i) > maxT2_rel+10),
                    fired_at_least = 0;
                    if pocet_minim_Ton4>1,
                        if (find((posT_on_minima4<posT_on_minima6(i)+offset_min)&(posT_on_minima4>posT_on_minima6(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton2>1,
                        if (find((posT_on_minima2<posT_on_minima6(i)+offset_min)&(posT_on_minima2>posT_on_minima6(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton8>1,
                        if (find((posT_on_minima8<posT_on_minima6(i)+offset_min)&(posT_on_minima8>posT_on_minima6(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton10>1,
                        if (find((posT_on_minima10<posT_on_minima6(i)+offset_min)&(posT_on_minima10>posT_on_minima6(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if ((fired_at_least >= 3)&&(byl_jsem_tu == 0)),
                        T_off = posT_on_minima6(i);
                        byl_jsem_tu = 1;
                    end
                end
            end
        case 4,
            for i = 1:length(posT_on_minima8),
                if (posT_on_minima8(i) > maxT2_rel+10),
                    fired_at_least = 0;
                    if pocet_minim_Ton4>1,
                        if (find((posT_on_minima4<posT_on_minima8(i)+offset_min)&(posT_on_minima4>posT_on_minima8(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton6>1,
                        if (find((posT_on_minima6<posT_on_minima8(i)+offset_min)&(posT_on_minima6>posT_on_minima8(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton2>1,
                        if (find((posT_on_minima2<posT_on_minima8(i)+offset_min)&(posT_on_minima2>posT_on_minima8(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton10>1,
                        if (find((posT_on_minima10<posT_on_minima8(i)+offset_min)&(posT_on_minima10>posT_on_minima8(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if ((fired_at_least >= 3)&&(byl_jsem_tu == 0)),
                        T_off = posT_on_minima8(i);
                        byl_jsem_tu = 1;
                    end
                end
            end
        case 5,
            for i = 1:length(posT_on_minima10),
                if (posT_on_minima10(i) > maxT2_rel+10),
                    fired_at_least = 0;
                    if pocet_minim_Ton4>1,
                        if (find((posT_on_minima4<posT_on_minima10(i)+offset_min)&(posT_on_minima4>posT_on_minima10(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton6>1,
                        if (find((posT_on_minima6<posT_on_minima10(i)+offset_min)&(posT_on_minima6>posT_on_minima10(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton8>1,
                        if (find((posT_on_minima8<posT_on_minima10(i)+offset_min)&(posT_on_minima8>posT_on_minima10(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if pocet_minim_Ton2>1,
                        if (find((posT_on_minima2<posT_on_minima10(i)+offset_min)&(posT_on_minima2>posT_on_minima10(i)-offset_min))),
                            fired_at_least = fired_at_least + 1;
                        end
                    end
                    if ((fired_at_least >= 3)&&(byl_jsem_tu == 0)),
                        T_off = posT_on_minima10(i);
                        byl_jsem_tu = 1;
                    end
                end
            end
    end
end

if ~exist('T_off','var');
    disp('T-off was not found in find_T_off...')
    T_off = length(Tecg);
    T_error = 1;
end