梯形法求定积分（一）设计梯形法求积分的类模板，梯形法

/*设计梯形法求积分的类模板，梯形法求积分的函数被定义为成员函数，可以求任意函数的定积分，用积分类的模板参数T引入被积函数*/
#include<iostream>
#include<cmath>
#include<iomanip>
using namespace std;
class Mysin{
public:
 double fun(double x)
 {
  return (sin(x));
 }
};
class F1{
public:
 double fun(double x)
 {
  return (1+x+2*x*x);
 }
};
class F2{
public:
 double fun(double x)
 {
  return (1+x+2*x*x+3*x*x*x);
 }
};
template<class T>
class JiFenLei{
private:
 double a,b,h,result;
 int n;
 T cf;
public:
 JiFenLei(double x=0,double y=0,int m=100)
 {
  a=x;
  b=y;
  n=m;
  jifenhanshu();//初始化后就地再类的内部进行积分，
  //但是感觉这样一定要注意构造函数的调用情况，千万不可随意调用构造函数
 }
 void print()
 {
  cout<<"定积分的值为:"<<result<<endl;
 }
 void xiugaitiaojian(double x=0,double y=0,int m=100)
 {
  a=x;
  b=y;
  n=m;
 }
 void jifenhanshu();
};
template<class T>
void JiFenLei<T>::jifenhanshu()
{
  h=(b-a)/n;
  int i;
  result=(cf.fun(a)+cf.fun(b))/2;
  for(i=1;i<n;i++)
   result += cf.fun(a+i*h);
  result = result*h;
}
void main()
{
 JiFenLei<Mysin>J1(0.0,3.0,100);
 J1.print();
 JiFenLei<F1>J2(0.0,3.0,100);
 J2.print();
}