#include #include #include #define fun(X) X+log(X) /* bousou simasu. */ double X; main() { int i=0; double V=1.0e-6,A=1.0e-2,W=1.0e-1,C,E; /* char B[30]="(X-1)*(X-2)"; */ double G,F; int Z=0; double GOSUB830(double, double); printf("Content-type: text/html\n\n"); printf("\n"); G=V; F=V; for (i=0;i<10;i++) { if (Z==0) { C=G; Z=1; } G=G-W; C=G; /* X=C; printf("X=%f\n",X); */ E=GOSUB830(C,A); printf("%f\n",E); F=F+W; C=F; E=GOSUB830(C,A); printf("%f\n",E); } printf("\n"); } double GOSUB830(double H,double M) { double Y,D; double B; do { X=H; /* printf("X2=%f\n",X); */ B=fun(X); Y=B; /* printf("Y=%f\n",Y); */ X=M+H; B=fun(X); D=H; H=D-M*Y/(B-Y); } while (fabs(D-H)>0.00000001) ; return (H); }