标题:用矩形法求定积分输出结果不对?
只看楼主
zhengzhihui
Rank: 2
来 自:汕头大学
等 级:论坛游民
帖 子:36
专家分:57
注 册:2009-7-4
结帖率:90%
已结贴  问题点数:10 回复次数:3 
用矩形法求定积分输出结果不对?
写一个用矩形法求定积分的通用函数,所谓矩形法就是根据极限法则,把欲求的面积分成n份矩形,然后累加起来。
以下程序为求sin,cos,exp的定积分程序,代码本身没错,可输出结果却不对??


#include<stdio.h>
#include<math.h>
#define n 1000
float sum1=0,sum2=0,sum3=0;
float collect_fun(float up,float down)
{
float i,temp;
temp=up-down;
for(i=1;i<=n;i++)
  {
    sum1=sum1+sin(i*temp/n)/n*temp;   
    sum2=sum2+cos(i*temp/n)/n*temp;
    sum3=sum3+exp(i*temp/n)/n*temp;
   }
}
main()
{
float up,down;
scanf("%f%f",&up,&down);
collect_fun(up,down);
printf("%f,%f,%f",sum1,sum2,sum3);
getch();
}

[ 本帖最后由 zhengzhihui 于 2009-8-28 14:45 编辑 ]

[ 本帖最后由 zhengzhihui 于 2009-8-28 14:49 编辑 ]
搜索更多相关主题的帖子: 矩形 积分 结果 输出 
2009-08-28 14:42
arthaszu
Rank: 3Rank: 3
等 级:论坛游侠
帖 子:90
专家分:163
注 册:2009-6-29
得分:7 
楼主,你代码本身就错了。
sin和cos是三角函数,所以sin(x+PI/2)=cos(x)
你写的程序应该是单独求sin函数,cos函数和指数函数的定积分程序。
因为这三种函数都是连续曲线,所以LZ想写成∫ _a^b(f(x) dx ) = F(a)-F(b) 的形势吧?
PS:由于本人还没高中毕业,所以我也不清楚应该怎么写这个代码,然,以你之所学应该可以完成的这个任务的。


To  four  years  in  each  other's
2009-08-28 23:45
zhengzhihui
Rank: 2
来 自:汕头大学
等 级:论坛游民
帖 子:36
专家分:57
注 册:2009-7-4
得分:0 
当然不是想写“∫ _a^b(f(x) dx ) = F(a)-F(b) ”了,都说了是矩形法,就是分成n个小区间再累加!!!
2009-08-30 14:28
zhengzhihui
Rank: 2
来 自:汕头大学
等 级:论坛游民
帖 子:36
专家分:57
注 册:2009-7-4
得分:0 
???
2009-08-30 14:29



参与讨论请移步原网站贴子:https://bbs.bccn.net/thread-283968-1-1.html




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

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