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

[复制链接]

字体大小: 正常放大

maleesky

3 主题	2 听众	19 积分

升级 14.74%

该用户从未签到

电梯直达

1^#

发表于 2005-4-1 16:34 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

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

//////////////////////////////////////////////////////////////////////
// 龙贝格求积法
//
// 调用时，须覆盖计算函数f(x)值的虚函数double Func(double x)
//
// 参数：
// 1. a - Double型变量，积分下限
// 2. b - Double型变量，积分上限，要求b>a
// 3. eps - Double型变量，积分精度要求
//
// 返回值：double 型，积分值
//////////////////////////////////////////////////////////////////////
double 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>
<

> // 迭代初值
h=b-a;
y[0]=h*(Func(a)+Func(b))/2.0;
m=1;
n=1;
ep=eps+1.0;

// 迭代计算
while ((ep>=eps)&&(m<=9))
{
  p=0.0;
      for (i=0;i<=n-1;i++)
      {
x=a+(i+0.5)*h;
         p=p+Func(x);
      }

  p=(y[0]+h*p)/2.0;
      s=1.0;
      for (k=1;k<=m;k++)
      {
s=4.0*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;
}

return(q);
}</P>

zan