function eps = macheps()
    
    one = 1;
    for i = 1:100
    
        if ( (1+2^(-i)) == 1 )
            eps = 2^(-i+1);
            break
        end
    end
end
