вычислить интеграл

#include <iostream.h>
#include <math.h>
 
ostream& operator<<(ostream &stream,char* s){
    for(char* ps=s; *ps; ps++){
        if(*ps=='?')
            stream<<char(241);
        else if(*ps=='?')
            stream<<char(240);
        else if(*ps>=-64 && *ps<=-17)           
            stream<<char(*ps+64+128);
        else if(*ps<0)
            stream<<char(*ps+64+176);
        else
            stream<<*ps;
    }
    return stream;
}
 
double f(double x)
{
    return exp(-pow(x,2));
}
 
int main()
{
    char ch;
    double a = 0,b = 0,h = 0,x,S = 0;
    cout<<"\tВведите диапазон [a;b] интегрирования\r\n";
    cout<<"a = ";cin>>a;
    cout<<"b = ";cin>>b;
    cout<<"Введите шаг интегрирования\r\n";
    cout<<"h = ";cin>>h;
    x = a;
    do
    {
        while(x < b)
        {
            S = S + 0.5*(f(x + h) + f(x))*h;
            x = x + h;
        }
        cout<<"Интеграл $f(x) равен = "<<S<<"\r\n";
        cout<<"Y - ввести другой диапазон интегрирование\r\n";
        cin>>ch;
    }
    while(ch == 'Y' || ch == 'y');
    return 0;
}