急！！！利用函数指针变量编写一个求定积分的通用函数，

答案:4  信息版本：手机版

解决时间 2021-05-05 09:17

已解决

2021-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。 谢谢大神了， 我急用

每一步最好有注释。 谢谢了。

最佳答案

2021-05-05 02:41

//利用Simpson公式来求定积分

#include

#include

double T(double x,double y,int z,double (*fun)(double)) ;

double integral (double a,double b,double (*fun)(double));

double f1(double t);

double f2(double t);

double f3(double t);

double f4(double t);

double f5(double t);

int n=0; //用来记录积分区间划分的间隔数，数量越大，越精确

void main()

{

double a,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);

}

double f1(double t)

{

return 1+t;

}

double f2(double t)

{

return 3+2*t;

}

double f3(double t)

{

return pow(2.71828,t)+1;//自然常数e，取了一个近似值2.71828

}

double f4(double t)

{

return (1+t)*(1+t);

}

double f5(double t)

{

return t*t*t;

}

//辛普森公式：

double T(double x,double y,int z,double (*fun)(double))

{

double h,Tn;

int i;

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);

}

double integral(double x,double y,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楼

2021-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楼

2021-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楼

2021-05-05 03:53