Метод трапеции на C++

[sir.andrey]

Уважаемые дамы и господа!!!

Я очень надеюсь на вашу помощь, посмотрите, в чем ошибка метода трапеции:
Надо найти площадь Q способом трапеции!
http://pics.vlfiles.ru/uploads/12840...d9abe9d288.png

double Trapec(){
        double a,b,x,dx,y1,y2,s1,s2,a2,b2;
        int i,n;
        n=10000;
        a=-1; b=0;
        dx=(b-a)/n;
        s1=0; x=a; i=0;
        while (i<n) {
                y1=f1(x)-exp(-x);
				y2=f1(x+dx)-exp(-x+dx);				
				if (y1>=0 && y2>=0) {s1=s1+(y1+y2)*dx;}
                x=x+dx; i++;
        }
                s1*=dx;
        a2=0; b2=1;
        dx=(b-a)/n;
        s2=0; x=a2; i=0;
        while (i<n) {
                y1=f2(x)-exp(-x);
				y2=f2(x+dx)-exp(-x+dx);
                if (y1>=0 && y2>=0) {s1=s1+(y1+y2)*dx;}
                x=x+dx; i++;
        }
                s2*=dx;
        return s1+s2;
 }