Курсовая работа: Розвязок інтеграла методом Чебиша Гауса Сімпсона

h=(b9-b0)/n;

s=F(m,b0)+F(m,b9);

k=1;

for(j=1;j<=n-1;j++){

a=a+h;

if(k!=1) s=s+2*F(m,a); else s=s+4*F(m,a);

k=-k;

}

s=s*h/3; f[i-1]=s;

cprintf("¦%2dп/6¦ %f ¦",i,s); gotoxy(1,6+i);

/*if(m==1)*/

if(i==1){

gotoxy(20,5+i);

ds=zr-s;

e=fabs(ds/zr)*100.0;

cprintf("табл. знач.=%f абс. пох.=%f вiдн. пох.=%f%%",zr,ds,e);

gotoxy(1,6+i);

}

}

cprintf("+-----------------+");

}

else{

b0=aS; b9=bS;

a=b0;

h=(b9-b0)/n;

s=F(m,b0)+F(m,b9);

k=1;

for(j=1;j<=n-1;j++){

a=a+h;