Информатика_1 / C / lecture7 / l7_2

/*1*/ int simpson(double a,double b,double e,double (*g)(double),double *s)
/*2*/ {int n=20,nmax=16384,i; double h,s0;
/*3*/ h=(b-a)/n;
/*4*/ for(i=1,*s=0;i<=n-1;i++)if(i%2==0)*s+=2*g(a+i*h); else *s+=4*g(a+i*h);
/*5*/ *s+=g(a)+g(b); *s*=h/3;
/*6*/ do { s0=*s; n*=2; h/=2;
/*7*/ for(i=1,*s=0;i<=n-1;i++)if(i%2==0)*s+=2*g(a+i*h); else *s+=4*g(a+i*h);
/*8*/ *s+=g(a)+g(b); *s*=h/3;
/*9*/ if(fabs(*s-s0)<e)return 0;
/*10*/ } while(n<nmax);
/*11*/ // ў аЁ©­л© ўл室
/*12*/ return 1;
/*13*/ }
/*14*/