  a=eps/4;
  y=(2-(1/a))/(1-(1/a));
  x=(1-y)/a;
  fprintf('  x = %.52e \n',x);
  fprintf('  y = %.52e \n',y);
  fprintf('1-u = %.52e \n',1-eps/2);