多态的方法求定积分

#include <iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
using namespace std;
void menu1();
void menu2();
class FUN
{
public:
    virtual double operator()(double x) const=0;
    virtual ~FUN(){}
};
class Function1:public FUN
{
public:
    virtual double operator()(double x) const
    {
        return x;
    }
};
class Function2:public FUN
{
public:
    virtual double operator()(double x) const
    {
        return sin(x);
    }
};
class Function3:public FUN
{
    virtual double operator ()(double x) const
    {
        return exp(x);
    }
};
class Integration
{
public:
    virtual double operator()(double a,double b,double n) const=0;
    virtual ~Integration(){}
};
class Trapzt:public Integration
{
private:
    const FUN &f;
public:
    Trapzt(const FUN &f):f(f){}
    virtual double operator()(double a,double b,double n) const
    {
        double h=fabs(b-a)/n,s=f(a);
        int i;
        for(i=1;i<=n;i++)
            s+=2*f(a+i*h);
        s-=f(b);
        s=s*h/2;
        return s;
    }
};
class Trapzj:public Integration
{
private:
    const FUN &f;
public:
    Trapzj(const FUN &f):f(f){}
    virtual double operator()(double a,double b,double n) const
    {
        double h=fabs(b-a)/n,s=0;
        int i;
        for(i=0;i<n;i++)
            s+=f(a+i*h+h/2);
        s*=h;
        return s;
    }
};
class Trapzs:public Integration
{
private:
    const FUN &f;
public:
    Trapzs(const FUN &f):f(f){}
    virtual double operator()(double a,double b,double n) const
    {
        double h=fabs(b-a)/n,s=(f(a)+f(b))/2;
        int i;
        for(i=1;i<n;i++)
            s+=f(a+i*h);
        s*=h;
        return s;
    }
};
int main()
{
    Function1 f1;
    Function2 f2;
    Function3 f3;
    Trapzt t1(f1);
    Trapzt t2(f2);
    Trapzt t3(f3);
    Trapzj j1(f1);
    Trapzj j2(f2);
    Trapzj j3(f3);
    Trapzj s1(f1);
    Trapzj s2(f2);
    Trapzj s3(f3);
    while(1)
    {
        system("cls");
        menu1();
        cout<<"请输入你想要积分函数的序号:"<<endl;
        int k,l,a,b,flag=1;
        double n;
        char c;
        cin>>k;
        switch(k)
        {
            case 1:while(flag)
            {
                system("cls");
                cout<<"f(x)=x"<<endl;
                menu2();
                cin>>l;
                switch(l)
                {
                case 1:cout<<"梯形法: "<<endl;
                       cout<<"请输入积分上下线和要分的个数"<<endl;
                       cin>>a>>b>>n;
                       cout<<"积分得: "<<t1(a,b,n)<<endl;
                       c=getchar();
                       system("pause");
                       break;
                case 2:cout<<"矩形法: "<<endl;
                       cout<<"请输入积分上下线和要分的个数"<<endl;
                       cin>>a>>b>>n;
                       cout<<"积分得: "<<j1(a,b,n)<<endl;
                       c=getchar();
                       system("pause");
                       break;
                case 3:cout<<"simpson: "<<endl;
                       cout<<"请输入积分上下线和要分的个数"<<endl;
                       cin>>a>>b>>n;
                       cout<<"积分得: "<<s1(a,b,n)<<endl;
                       c=getchar();
                       system("pause");
                       break;
                case 4:flag=0;break;
                }
            }
            break;
            case 2:while(flag)
            {
                system("cls");
                cout<<"f(x)=sin(x)"<<endl;
                menu2();
                cin>>l;
                switch(l)
                {
                    case 1:cout<<"梯形法: "<<endl;
                           cout<<"请输入积分上下线和要分的个数"<<endl;
                           cin>>a>>b>>n;
                           cout<<"积分得: "<<t2(a,b,n)<<endl;
                           c=getchar();
                           system("pause");
                           break;
                    case 2:cout<<"矩形法: "<<endl;
                           cout<<"请输入积分上下线和要分的个数"<<endl;
                           cin>>a>>b>>n;
                           cout<<"积分得: "<<j2(a,b,n)<<endl;
                           c=getchar();
                           system("pause");
                           break;
                    case 3:cout<<"simpson: "<<endl;
                           cout<<"请输入积分上下线和要分的个数"<<endl;
                           cin>>a>>b>>n;
                           cout<<"积分得: "<<s2(a,b,n)<<endl;
                           c=getchar();
                           system("pause");
                           break;
                    case 4:flag=0;break;
                   }
            }
            break;
             case 3:while(flag)
             {
                 system("cls");
                 cout<<"f(x)=e^x"<<endl;
                 menu2();
                 cin>>l;
                 switch(l)
                    {
                        case 1:cout<<"梯形法: "<<endl;
                               cout<<"请输入积分上下线和要分的个数"<<endl;
                               cin>>a>>b>>n;
                               cout<<"积分得: "<<t3(a,b,n)<<endl;
                               c=getchar();
                               system("pause");
                               break;
                        case 2:cout<<"矩形法: "<<endl;
                               cout<<"请输入积分上下线和要分的个数"<<endl;
                               cin>>a>>b>>n;
                               cout<<"积分得: "<<j3(a,b,n)<<endl;
                               c=getchar();
                               system("pause");
                               break;
                        case 3:cout<<"simpson: "<<endl;
                               cout<<"请输入积分上下线和要分的个数"<<endl;
                               cin>>a>>b>>n;
                               cout<<"积分得: "<<s3(a,b,n)<<endl;
                               c=getchar();
                               system("pause");
                               break;
                        case 4:flag=0;break;
                   }
             }
                   break;
              case 4:cout<<endl<<"Thanks for use!!"<<endl<<endl;
                     exit(-1);
        }
    }
    return 0;
}
void menu1()
{
    cout<<"                   求定积分"<<endl<<endl;
    cout<<"       1.f(x)=x"<<endl;
    cout<<"       2.f(x0=sin(x)"<<endl;
    cout<<"       3.f(x)=e^x"<<endl;
    cout<<"       4.退出"<<endl;
}
void menu2()
{
    cout<<"     请输入你要用的方法序号:"<<endl<<endl;
    cout<<"1.梯形法:"<<endl;
    cout<<"2.矩形法:"<<endl;
    cout<<"3.Simpson:"<<endl;
    cout<<"4.返回:"<<endl;
}