急！！！利用函数指针变量编写一个求定积分的通用函数，
答案:4信息版本：手机版
解决时间 -05-05 09:17
已解决
-05-05 02:15
用它分别求5个函数的定积分:每次需要求定积分的函数是不一样的。可以编写一个求定积分的通用函数integral，它有3个形参： 下限a、上限b，以及指向函数的指针变量fun。函数原型可写为
double integral (double a，double b，double (*fun)(double));
分别编写5个函数f1，f2，f3，f4，f5， 用来求上面5个函数的值。然后先后调用integral函数5次，每次调用时把a，b以及f1，f2，f3，f4，f5之一作为实参，即把上限、下限以及有关函数的入口地址传送给形参fun。在执行integral函数过程中求出各函数定积分的值。 五个函数分别。为: 1+x; 2x+3; e的x次方+1；(1+x)*(1+x); x*x*x。 谢谢大神了， 我急用
每一步最好有注释。 谢谢了。
最佳答案
-05-05 02:41
//利用Simpson公式来求定积分
#include
#include
doubleT(doublex,doubley,intz,double(*fun)(double));
doubleintegral(doublea,doubleb,double(*fun)(double));
doublef1(doublet);
doublef2(doublet);
doublef3(doublet);
doublef4(doublet);
doublef5(doublet);
intn=0;//用来记录积分区间划分的间隔数，数量越大，越精确
voidmain()
{
doublea,b,s;
printf("积分下限a:\n");
scanf("%lf",&a);
printf("积分上限b:\n");
scanf("%lf",&b);
printf("区间等分个数n:\n");
scanf("%d",&n);
s=integral(a,b,f1);//用函数f1来验证
printf("函数f(x)在区间%f到%f的积分值为s=%f\n",a,b,s);
}
doublef1(doublet)
{
return1+t;
}
doublef2(doublet)
{
return3+2*t;
}
doublef3(doublet)
{
returnpow(2.71828,t)+1;//自然常数e，取了一个近似值2.71828
}
doublef4(doublet)
{
return(1+t)*(1+t);
}
doublef5(doublet)
{
returnt*t*t;
}
//辛普森公式：
doubleT(doublex,doubley,intz,double(*fun)(double))
{
doubleh,Tn;
inti;
h=(y-x)/z;
Tn=(fun(x)+fun(y))/2;
for(i=1;i
Tn=Tn+fun(x+i*h);
Tn=Tn*h;
return(Tn);
}
doubleintegral(doublex,doubley,double(*fun)(double))
{
return(4*T(x,y,2*n,fun)-T(x,y,n,fun))/3;
}
验证结果：
通过数学知识，我们可以知道
f1(t)=1+t;(积分应该是t+0.5*t*t+任意常数)在区间0到1之间定积分的确是1.5
全部回答
1楼
-05-05 05:17
//简单的写了一下按照最原始的定积分思想来的
#include
//定积分通用函数
double integral (double a,double b,double (*fun)(double)); //函数声明
double func1(double x); //函数声明
// 主函数
void main()
{
double dintegral;
double a = 1,b = 2; //上下限
dintegral = integral(a,b,func1); //调用定积分函数
printf("积分上限：%.2f,下限 %.2f所求函数f(x) = x + 1 积分为：%.2f \n",a,b,dintegral);
}
//定积分函数 a下限 b 上限 (*fun)(double) 函数指针 梯形法求积分
double integral (double a,double b,double (*fun)(double))
{
//设置积分精度
int iAccuracy = 10000;// 及把给定区域10000等分
double dPart = (b - a) / 10000;
double dTotalArea = 0;
//循环求面积
for(int iPartNum = 1; iPartNum <= iAccuracy; iPartNum ++ )
{
//一个梯形面积 (上底加下底)* 高 /2
double dArea = ( fun(a) + fun(a+dPart) ) * dPart /2;
a+=dPart;
dTotalArea += dArea;
}
//返回所求积分
return dTotalArea;
}
//函数取一个做示范 其余类似
double func1(double x)
{
return (x + 1);
}
2楼
-05-05 04:43
//利用Simpson公式来求定积分
#include
#include
double f1(double);
double f2(double);
double f3(double);
double f4(double);
double f5(double);
double T(double x,double y,int z,double (*fun)(double));
double integral(double a,double b,double (*fun)(double));
int n=0;//用来记录积分区间划分的间隔数，数量越大，越精确
int main()
{
double a,b,s1,s2,s3,s4,s5;
printf("积分上限 a：");
scanf("%lf",&a);
printf("积分下限 b：");
scanf("%lf",&b);
printf("区间等分个数 n:");
scanf("%d",&n);
s1=integral(a,b,f1);
s2=integral(a,b,f2);
s3=integral(a,b,f3);
s4=integral(a,b,f4);
s5=integral(a,b,f5);
printf("函数f1在区间%f到%f 的积分值为 s=%f\n",a,b,s1);
printf("函数f2在区间%f到%f 的积分值为 s=%f\n",a,b,s2);
printf("函数f3在区间%f到%f 的积分值为 s=%f\n",a,b,s3);
printf("函数f4在区间%f到%f 的积分值为 s=%f\n",a,b,s4);
printf("函数f5在区间%f到%f 的积分值为 s=%f\n",a,b,s5);
return 0;
}
double f1(double x)
{
return x+1;
}
double f2(double x)
{
return 2*x+3;
}
double f3(double x)
{
return pow(2.71828,x)+1;
}
double f4(double x)
{
return pow(x+1,2);
}
double f5(double x)
{
return pow(x,3);
}
double T(double x,double y,int z,double (*fun)(double))
{
double h,Tn;
int i;
h=(x+y)/z;
Tn=((*fun)(x)+(*fun)(y))/2;
for(i=1;i
Tn+=(*fun)(x+i*h);
Tn=Tn*h;
return Tn;
}
double integral(double x,double y,double (*fun)(double))
{
return (4*T(x,y,2*n,fun)-T(x,y,n,fun))/3;
}
//利用Simpson公式来求定积分
#include
#include
double f1(double);
double f2(double);
double f3(double);
double f4(double);
double f5(double);
double T(double x,double y,int z,double (*fun)(double));
double integral(double a,double b,double (*fun)(double));
int n=0;//用来记录积分区间划分的间隔数，数量越大，越精确
int main()
{
double a,b,s1,s2,s3,s4,s5;
printf("积分上限 a：");
scanf("%lf",&a);
printf("积分下限 b：");
scanf("%lf",&b);
printf("区间等分个数 n:");
scanf("%d",&n);
s1=integral(a,b,f1);
s2=integral(a,b,f2);
s3=integral(a,b,f3);
s4=integral(a,b,f4);
s5=integral(a,b,f5);
printf("函数f1在区间%f到%f 的积分值为 s=%f\n",a,b,s1);
printf("函数f2在区间%f到%f 的积分值为 s=%f\n",a,b,s2);
printf("函数f3在区间%f到%f 的积分值为 s=%f\n",a,b,s3);
printf("函数f4在区间%f到%f 的积分值为 s=%f\n",a,b,s4);
printf("函数f5在区间%f到%f 的积分值为 s=%f\n",a,b,s5);
return 0;
}
double f1(double x)
{
return x+1;
}
double f2(double x)
{
return 2*x+3;
}
double f3(double x)
{
return pow(2.71828,x)+1;
}
double f4(double x)
{
return pow(x+1,2);
}
double f5(double x)
{
return pow(x,3);
}
double T(double x,double y,int z,double (*fun)(double))
{
double h,Tn;
int i;
h=(x+y)/z;
Tn=((*fun)(x)+(*fun)(y))/2;
for(i=1;i
Tn+=(*fun)(x+i*h);
Tn=Tn*h;
return Tn;
}
double integral(double x,double y,double (*fun)(double))
{
return (4*T(x,y,2*n,fun)-T(x,y,n,fun))/3;
}
3楼
-05-05 03:53
