/*#include <stdio.h>*/
/*#include <stdlib.h>*/
/*#include <math.h>*/
/*#define fun(x) x+log(x)*/
double x;
main()
{
 int i=0;
 double v=1.0e-6,a=1.0e-7;
 double w=1.0e-9,c,e;
 double g,f;
 int z=0;
 double gosub830();

 g=v;
 f=v;
 for (i=0;i<10;i++) {
  if (z==0) {
   c=g;
   z=1;
  }
 g=g-w;
 c=g;

 e=gosub830(c,a);
 printf("%f\n",e);

 f=f+w;
 c=f;

 e=gosub830(c,a);
 printf("%f\n",e);
 }
}

double gosub830(h,m)
double h,m;
{
 double y,d;
 double b;

 do {
  x=h;

  b=x+log(x);
  y=b;

  x=m+h;
  b=x+log(x);
  d=h;
  h=d-m*y/(b-y);
 } while ((d-h)>0.00000001 || (h-d)>0.00000001);
 return (h);
}