一个定积分计算的类的实现!经典!
2005-01-29 18:16
211 查看
可以用于近似计算一元函数的定积分,类的实现上采用了函数指针类型,使得应用更为灵活!
typedef double (*pfn)(double);
class integral{
double lower;//积分下限
double upper;//积分上限
pfn integrand;//积分函数
public:
integral(double a,double b,pfn f){
lower=a;
upper=b;
integrand=f;
}
double lowbd()const {return lower;}
double upbd()const {return upper;}
void changebd(double,double);//改变积分上下限
double trapezoidal(int) const;//积分求解,参数限定了步长的大小
friend double simpson(integral,int);//友元函数声明
};
inline void integral::changebd(double a,double b)
{
lower=a;
upper=b;
}
double integral::trapezoidal(int n) const
{
double h=(upper-lower)/n;
double sum=integrand(lower)*0.5;
for(int i=1;i<n;i++)sum+=integrand(lower+i*h);
sum+=integrand(upper)*0.5;
return sum*h;
}
double simpson(integral ig,int n)
{
double h=(ig.upper-ig.lower)/n;
double sum=ig.integrand(ig.lower)*0.5;
for(int i=1;i<n;i++)
sum+=ig.integrand(ig.lower+i*h);
sum+=ig.integrand(ig.upper)*0.5;
double summid=0.0;
for(i=1;i<=n;i++)
summid+=ig.integrand(ig.lower+(i-0.5)*h);
return (sum+2*summid)*h/3.0;
}
typedef double (*pfn)(double);
class integral{
double lower;//积分下限
double upper;//积分上限
pfn integrand;//积分函数
public:
integral(double a,double b,pfn f){
lower=a;
upper=b;
integrand=f;
}
double lowbd()const {return lower;}
double upbd()const {return upper;}
void changebd(double,double);//改变积分上下限
double trapezoidal(int) const;//积分求解,参数限定了步长的大小
friend double simpson(integral,int);//友元函数声明
};
inline void integral::changebd(double a,double b)
{
lower=a;
upper=b;
}
double integral::trapezoidal(int n) const
{
double h=(upper-lower)/n;
double sum=integrand(lower)*0.5;
for(int i=1;i<n;i++)sum+=integrand(lower+i*h);
sum+=integrand(upper)*0.5;
return sum*h;
}
double simpson(integral ig,int n)
{
double h=(ig.upper-ig.lower)/n;
double sum=ig.integrand(ig.lower)*0.5;
for(int i=1;i<n;i++)
sum+=ig.integrand(ig.lower+i*h);
sum+=ig.integrand(ig.upper)*0.5;
double summid=0.0;
for(i=1;i<=n;i++)
summid+=ig.integrand(ig.lower+(i-0.5)*h);
return (sum+2*summid)*h/3.0;
}
相关文章推荐
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。 假设每度电的价格为1.2元,计算并显示本月电费
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数, 假设每度电的价格为1.2元,计算并显示本月电费。
- 习题 9.2 写一个函数days, 实现第1题的计算。由主函数将年、月、日传递给days函数,计算后将日子数传回主函数输出。
- javascript 网页中实现一个计算当年还剩多少时间的倒数计时程序(转)
- 一个简单的大数实现方案,计算斐波纳契数列
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。 假设每度电的价格为1.2元,计算并显示本月电费
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。
- 一个简单的Java计算器实现
- Java经典实例:实现一个简单堆栈
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。 假设每度电的价格为1.2元,计算并显示本月电费
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。 假设每度电的价格为1.2元,计算并显示本月电费
- [家里蹲大学数学杂志]第425期一个定积分的计算
- C#实现计算一个点围绕另一个点旋转指定弧度后坐标值的方法
- 赛码网,简单计算,用C语言自己实现的一个程序,仅供参考
- 编写一个函数days,实现上面的计算。。由主函数将年月日传递给days,计算出该日在本年中是第几天并将结果传回主函数输出
- 【转载】如何搭建一个快速的分布式计算平台?Parallel python提供了简易的方式来实现此目的。
- 【java学习记录】2.定义一个计算矩形面积、立方体和球体体积的类,该类完成计算的方法用静态方法实现
- 编写一个程序,实现设置上月、本月电表读数,显示上月、本月电表读数,计算并显示本月用电数。 假设每度电的价格为1.2元,计算并显示本月电费
- 白话经典算法系列之中的一个 冒泡排序的三种实现
- 在linux中用C程序实现读取一个文件中的字符,计算它的个数,并把个数显示在另外一个文件夹内