Geum.ru

Вычисление интеграла с помощью метода трапеций на компьютере

Курсовой проект - Компьютеры, программирование

Другие курсовые по предмету Компьютеры, программирование

 

 

 

 

 

 

 

i=1

i=n/2

 

 

 

 

 

 

 

 

 

 

 

 

 

  1. :

 

 

 

#include

#include

#include

main()

{

double a,b,er,eps,f(double),s,trap(double,double,double,double(*)(double));

clrscr();

printf("\n : ");

scanf ("%lf%lf%lf",&a,&b,&eps);

s=trap(a,b,eps,f);

printf("\n a=%3.2lf b=%3.2lf %lf",a,b,s);

getch();

}

double f(double x)

{

return x*x*x+2*(x*x)-3*x-8;

}

double trap(double a,double b,double eps,double(*f)(double))

{

double h,s,s0,s1,sn;

int i,n;

s=1; sn=101;

n=4;

s0=(f(a)+f(b))/2;

s1=f((a+b)/2);

while(fabs(s-sn)>eps){

sn=s;

h=(b-a)/n;

for(i=0; i<n/2; i++)

s1+=f(a+(2*i+1)*h);

s=h*(s0+s1);

n*=2;

}

return s;

}

 

 

  1. :

 

nint [a, b]iint adouble bdouble hdouble epsdouble fdouble(*) - xdouble - fsdouble s0double a bs1double sndouble

 

  1. :

 

XiYi0-80,75-8,7031251,5-4,6252,256,7656253283,7561,6093754,5110,1255,25176,0781256262

 

 

 

 

 

 

  1. :

eps = 0.1 eps = 0.001

a, b, eps: a, b, eps:

0 0

6 6

.1 .001

= 366.024170 = 366.000094

 

. 1/10000.