10 double x;
20 main()
30 {
40  int i=0;
50  double v=1.0e-6,a=1.0e-7;
60  double w=1.0e-9,c,e;
70  double g,f;
80  int z=0;
90  double gosub830();
100  g=v;
110  f=v;
120  for (i=0;i<10;i++) {
130   if (z==0) {
140    c=g;
141    z=1;
142  }
143  g=g-w;
144  c=g;
150  e=gosub830(c,a);
160  printf("%f\n",e);
170  f=f+w;
180  c=f;
190  e=gosub830(c,a);
200  printf("%f\n",e);
210  }
220 }
230 double gosub830(h,m)
231 double h,m;
240 {
250  double y,d;
260  double b;
270  while ((d-h)>0.00000001||(h-d)>0.00000001) {
280   x=h;
290   b=x+log(x);
300   y=b;
310   x=m+h;
320   b=x+log(x);
330   d=h;
340   h=d-m*y/(b-y);
350  }
360  return (h);
370 }