| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 683 人关注过本帖
标题:定积分的通用函数
取消只看楼主 加入收藏
303770957
Rank: 11Rank: 11Rank: 11Rank: 11
等 级:小飞侠
威 望:6
帖 子:838
专家分:2125
注 册:2005-9-10
收藏
 问题点数:0 回复次数:0 
定积分的通用函数

/* 题目:编写一个求定积分的通用函数,用它分别求以下5个函数的定积分:
∫(1+x)dx,∫(2x+3)xdx,∫(exp(x)+1)dx(exp表示e为底x为指数),∫(1+x)(1+x)dx,∫x*x*xdx .上-下限为a和b.
*/
程序清单如下:
#include "stdio.h"
#include "conio.h"
#include "math.h"

float f1(float x)
{return (1+x);}

float f2(float x)
{return (3+2*x);}

float f3(float x)
{return (1+exp(x));}

float f4(float x)
{return ((1+x)*(1+x));}

float f5(float x)
{return (x*x*x);}

float Integral(float a,float b,float (*fun)(float))
{
float fsum=0.0;
long N=400000L,i;
for(i=0;i<=N;i++)
fsum+=(*fun)(a+i*(b-a)/N);
return (fsum*(b-a)/N);
}

void main()
{
printf("f1=1+x 时∫f1dx=%.3lf\n",Integral(2,3.5,f1));
printf("f2=3+x*x 时∫f2dx=%.3lf\n",Integral(-6,1,f3));
printf("f4=(1+x)*(1+x) 时∫f4dx=%.3lf\n",Integral(3,4,f4));
printf("f5=x*x*x 时∫f5dx=%.3lf\n",Integral(1,2,f5));
getch();
}

运行结果如下:
*************************************
f1=1+x 时∫f1dx=5.625
f2=3+x*x 时∫f2dx=9.716
f4=(1+x)*(1+x) 时∫f4dx=20.333
f5=x*x*x 时∫f5dx=3.750
*************************************

搜索更多相关主题的帖子: 函数 积分 
2007-04-22 15:07
快速回复:定积分的通用函数
数据加载中...
 
   



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

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