function prvdvoj
% Tato funkce hleda prvociselna dvojcata(prvocilsa, jejichz rozdil je 2) po
% maximalni zadanou mez
hm = input('Zadej horni mez prvociselnych dvojcat: ');
celkem=0;
if (hm<0)           %kontrola, zda je dobre zadana horni mez
    disp('Byla zadana spatne horni mez');
else
if (hm==1)
    vysl=0;
    vysl1=0;
else vysl=[];
     vysl1=[];
end

for i=1:hm
    d=i-1;
    j=i-2;
    e=j-1;
    while (d>1)             %kontoroluje zda je cislo i prvocislo
        if (mod(i,d)~=0)
            d = d - 1;
        else d = 0;
        end
    end
    if (d==1) 
        while (e>1)         %kontroluje zda je cislo j prvocislo
            if (mod(j,e)~=0)
                e = e - 1;
            else e = 0;
            end
        end
        if (e==1)
            vysl = [vysl,i];    %ulozi vetsi z prvociselnych dvojcat
            vysl1 = [vysl1,j];  %ulozi mensi z prvociselnych dvojcat
            celkem = celkem + 1;
        end
    end
end
disp(sprintf('V intervalu od 1 do %d bylo nalezeno celkem %d prvociselnych dvojcat.',hm,celkem))
if (celkem~=0)
    vysl                %vypise vetsi cislo z prvociselnych dvojcat
    vysl1               %vypise mensi cislo z prvociselnych dvojcat
end
end
end