龙贝格求积法源代码,效率很好的

maleesky

<>积分的源代码,效率很好的

2 E, v. F- E, K7 ~' g6 x; X
! i2 m' |. W7 G* \& L5 I//////////////////////////////////////////////////////////////////////
2 L+ O) z% L# O// 龙贝格求积法
//
- J  ^4 I3 h+ X1 z. F: m& T// 调用时，须覆盖计算函数f(x)值的虚函数double Func(double x)
2 ]* }8 H: l; h- \4 [- h//
) T! k& @! X& r' ?/ X" d9 E2 W/ W// 参数：
- o) h& @, S. V// 1. a - Double型变量，积分下限
8 ?- y8 @! S# ]( g/ \( ]// 2. b - Double型变量，积分上限，要求b&gt;a
: O- c' X2 @2 X2 z+ H; ]  K// 3. eps - Double型变量，积分精度要求
1 a' \: D: W4 t/ f) C. M//
" n& Q* D$ P& P5 ?* `9 o// 返回值：double 型，积分值
5 ?1 s% h: X5 ^. w5 O//////////////////////////////////////////////////////////////////////
1 n0 U* s5 U$ |0 h5 idouble Integral_Romberg(double a, double b, double eps /*= 0.000001*/)
{
    int m,n,i,k;
    double y[10],h,ep,p,x,s,q;</P>
# C: f; w: X' d5 r( D6 ^) k" u<> // 迭代初值
    h=b-a;
    y[0]=h*(Func(a)+Func(b))/2.0;
* A& f  L6 k6 j+ V    m=1;
n=1;
ep=eps+1.0;
2 f& ]3 ]7 N3 Y9 i   
  e$ x5 D( Y2 O& y9 w4 r; g7 X // 迭代计算
" f8 [, h: w0 B6 B0 ` while ((ep&gt;=eps)&amp;&amp;(m&lt;=9))
    {
) W1 F6 ?7 G' @2 H1 R3 _( \  p=0.0;
, u) ~% q+ v4 F* O7 a        for (i=0;i&lt;=n-1;i++)
        {
: k9 N# E# j: S/ }! S  G   x=a+(i+0.5)*h;
            p=p+Func(x);
        }
6 |% M$ R( K. g7 E* P( e# M3 W1 i        
  p=(y[0]+h*p)/2.0;
        s=1.0;
, a% q  @$ d5 j6 j: c: N' I. w        for (k=1;k&lt;=m;k++)
        {
5 h  ?) h- ^& p% l   s=4.0*s;
0 x8 g) D3 `" e$ V% j& s            q=(s*p-y[k-1])/(s-1.0);
            y[k-1]=p; p=q;
        }</P>
<>        ep=fabs(q-y[m-1]);
        m=m+1;
  y[m-1]=q;
  n=n+n;
  h=h/2.0;
    }
0 Y# T/ t# e# L  l   
return(q);
% M5 N) i' G2 Q" d; a: N, D}</P>