function y=machepsdouble
    eps = 1.d+0;
    lasteps = 2.*eps;
    while ((1.d0+eps) > 1.d0) do
      lasteps = eps;
      eps = eps/2;
  end;
  y = lasteps;
  
endfunction

function y=machepszero
    eps = 1.d+0;
    lasteps = 2.*eps;
    while ((eps) > 0) do
      lasteps = eps;
      eps = eps/2;
  end;
  y = lasteps;
  
endfunction