function prvocisla
od = input('Od kolika se budou prvocisla pocitat: ');
do = input('Do kolika se budou prvocisla pocitat: ');

celkem = 0;

if (od==1) 
    Vysl = 1;
else
    Vysl = [];
end

for x = od:do
    d = x-1;    
    while (d>1)
        if (mod(x,d)~=0)
            d = d - 1;
        else d = 0;
        end
    end
    if (d==1) 
        Vysl = [Vysl,x];
        celkem = celkem + 1;
    end
end
disp(sprintf('V intervalu od %d do %d bylo nalezeno celkem %d prvocisel.',od,do,celkem))
plot(Vysl,'r+ ')
end