| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 682 人关注过本帖
标题:同样是求积分,为什么被积函数换了就不能运算了呢?
只看楼主 加入收藏
xiaolaopo
Rank: 1
等 级:新手上路
帖 子:20
专家分:0
注 册:2008-7-31
收藏
 问题点数:0 回复次数:0 
同样是求积分,为什么被积函数换了就不能运算了呢?
#include "math.h"
  double ffts(a,b,eps,f)
  double a,b,eps,(*f)();
  { int n,k;
    double fa,fb,h,t1,p,s,x,t;
    fa=(*f)(a); fb=(*f)(b);
    n=1; h=b-a;
    t1=h*(fa+fb)/2.0;
    p=eps+1.0;
    while (p>=eps)
      { s=0.0;
        for (k=0;k<=n-1;k++)
          { x=a+(k+0.5)*h;
            s=s+(*f)(x);
          }
        t=(t1+h*s)/2.0;
        p=fabs(t1-t);
        t1=t; n=n+n; h=h/2.0;
      }
    return(t);
  }


#include "stdio.h"
  #include "9ffts.c"(w为以上函数文件名"9ffts.c")
  main()
  { double a,b,eps,h,t1,t2,sum,fftsf(double);
    a=0.0; b=1.0; eps=0.000001;
    t1=ffts(a,b,eps,fftsf);
    a=1.0; b=2.0; eps=0.000001;
    t2=ffts(a,b,eps,fftsf);
   sum=t1+t2;
   printf("\n");
   printf("sum=%e\n",sum);
    printf("\n");
    return 0;
}

  #include "math.h"
  double fftsf(x)
  double x;
  { double y;
    y=(x*x);            (这里)
    return(y);
  }
当被积函数为y=x*x,结果为2.555556e+00
但被积函数为y=16.7*((2.0*(0.0003-0.5*x))/0.002-((0.0003-0.5*x)/0.002)*((0.0003-0.5*x)/0.002))怎么就不能算了呢?
搜索更多相关主题的帖子: 函数 eps 积分 double 运算 
2008-08-01 10:14
快速回复:同样是求积分,为什么被积函数换了就不能运算了呢?
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.018416 second(s), 9 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved