数值积分与微分上机报告.pdfVIP

数值积分

6.1目的与要求

（1）通过实际计算体会各种方法的精确度；

（2）会编写用算法求定积分的程序。

6.2复化

6.2.1算法

n−1

hh

复化为S=fx+4fx++fx，计算过程如下：

n(k)k(k+1)

6k=02

令=−=+=

1.hbans,fah2,s0;

()1()2

2.对=−计算

k1,2,,n1

s=s+fa+kh+hs=s+fa+kh

2,

11()22()

h

3.s=fa+4s+2s+fb。

(()12())

6

6.2.2程序与实例

1dx

例用复化计算积分I=03。

1+x

C语言程序如下：

#includestdio.h

#includemath.h

floatSimpson(float(*f)(float),floata,floatb,intn)

{

intk;

floats,s1,s2=0;

floath=(b-a)/n;

s1=(*f)(a+h/2);

for(k=1;k=n-1;k++)

{

s1+=(*f)(a+k*h+h/2);

s2+=(*f)(a+k*h);

}

s=h/6*((*f)(a)+4*s1+2*s2+(*f)(b));

returns;

}

floatf(floatx)

{

return1/sqrt(1+x*x*x);

}

voidmain()

{

inti,n=2;

floats;

for(i=0;i=4;i++)

{

s=Simpson(f,0,1,n);

printf(s(%d)=%f\n,n,s);

n*=2;

}

}

6.3算法

6.3.1算法

1n−1