Метод правых прямоугольников! Как исправить?

#include "stdafx.h" 
#include <stdlib.h> 
#include <iostream> 
using namespace std; 
int main() 
{ 
system("cls"); 
setlocale(LC_ALL, "RUS"); 
float a = 1, b = 4, s1, s2 = 0, dx; 
int n; 
//dx -длина отрезка разбиения [A;B] 
//s1- площадь трапеции 
//s2- площадь , ограниченная параболой y=x^2+1 
 
cout « "Введите количество разбиений отрезка [a;b]" 
« "\n(Ctrl + z Enter – выход из программы\n"; 
cin » n; 
while (!feof(stdin))// Начало цикла. Цикл выполняется до тех пор, 
// пока не будет введен символ Ctrl + z. 
{ 
//***** Алгоритмический блок *** 
dx = (b - a) / n; 
s1 = (2 + 1)*(b - a) / 2; 
cout « "\nПлощадь трапеции S1= " « s1 « "\n "; 
// метод прямоугольников по лев.концу отрезка dx 
for (float x = b; x > a; x -= dx) 
{ 
// s2 += 3 * pow(x, 2) - 4 * x + 2; 
s2 += pow(x, 1.5) / 1.5; 
} 
s2 = s2 * dx; 
cout « "\nПлощадь фигуры, образованной параболой S2= " « s2 « "\n "; 
s1 = s2 - s1; 
// ******************************* 
cout « "\nПлощадь искомой фигуры S= " « s1 « "\n "; 
system("pause"); 
cout « "Введи число разбиений отрезка\n(Сtrl + z Enter – выход из программы):\n "; 
cin » n; 
return 0; 
} 
}

#include "stdafx.h" 
#include <stdlib.h> 
#include <iostream> 
using namespace std;
int main()
{
    system("cls");
    setlocale(LC_ALL, "RUS");
    float a = -2, b = 3, s1, s2 = 0, dx;
    int n;
    //dx -длина отрезка разбиения [A;B] 
    //s1- площадь трапеции 
    //s2- площадь , ограниченная параболой y=x^2+1 
 
    cout << "Введите количество разбиений отрезка [a;b]"
        << "\n(Ctrl + z Enter – выход из программы\n";
    cin >> n;
    while (!feof(stdin))// Начало цикла. Цикл выполняется до тех пор, 
                        // пока не будет введен символ Ctrl + z. 
    {
        //***** Алгоритмический блок *** 
        dx = (b - a) / n;
        s1 = (22 + 17)*(b - a) / 2;
        cout << "\nПлощадь трапеции S1= " << s1 << "\n ";
        // метод прямоугольников по лев.концу отрезка dx 
        for (float x = a; x < b; x += dx)
        {
            s2 += 3 * pow(x, 2) - 4 * x + 2;
        }
        s2 = s2 * dx;
        cout << "\nПлощадь параболы S2= " << s2 << "\n ";
        s1 = s1 - s2;
        // ******************************* 
        cout << "\nПлощадь искомой фигуры S= " << s1 << "\n ";
        system("pause");
        cout << "Введи число разбиений отрезка\n(Сtrl + z Enter – выход из программы):\n ";
        cin >> n;
    }
    return 0;
}