C++实现计算简单定积分函数

#include "StdAfx.h"
#include <iostream>
#include <cmath>

using namespace std;

//根据最简单的定积分算法，编写一个定积分函数
//实现不同函数的定积分，利用函数指针作为形参
typedef double (*PF)(double x);

//a,b,n,fun 
double Integral(double a,double b,int n, PF fun)
{
	double step=(b-a)/n;
	double sum=0;
	for(int i=0;i<n;i++)
	{
		double x=a+step*i;
		sum+=step*(*fun)(x);
	}
	return sum;
}

//测试
double Sin(double x)
{
	return sin(x);
}

double Linear(double x)
{
	return x;
}

int main(int argc, char*argv[])
{
	cout<<Integral(0,1,10,Sin)<<endl;
	cout<<Integral(1,2,100,Linear)<<endl;
	return 0;
}