function strojeps

cislo   = 1;
epsilon = 0.1 * cislo;
while cislo+epsilon > cislo
    epsilon = epsilon / 2;
end
epsilon = epsilon * 2
eps(1)
