编制求解定积分的通用函数

来源:互联网 发布:什么是网络金融 编辑:程序博客网 时间:2024/05/08 06:38

2.已知两个函数f1(x)和f2(x),求它们在指定区间的定积分,要求编制求解定积分的通用函数(函数中可用函数型参数或函数指针型参数)

#include<stdio.h>

#include<math.h>

double f1(double x)

{

double s;

    s=x;

return s;

}

double f2(double x)

{

double s;

s=x*x;

return s;

}

double dIntegral(double(*p)(double) ,double a,double b,int n)

{

   double x,result=0;

   double h;

   h=(b-a)/n;

   x=a;

   for(int i=0;i<n;i++)

   {

   x+=h;

   result+=(*p)(x)*h;

   }

   return result;

}

void main(){

    double (*p)(double);

    double a=2,b=4,s;

int n=100000;

p=f1;

   s=dIntegral(p,a,b,n);

   printf("f1函数的定积分为:%lf\n",s);

   p=f2;

   s=dIntegral(p,a,b,n);

   printf("f2函数的定积分为:%lf\n",s);

}

原创粉丝点击