вычислить определенный интеграл)

[iris_ka]

подскажите, пожалуйста, как на С++ вычислить интеграл функции x/(x^2+1)^2 на промежутке от -1 до 2?) (с помощью разбиения на 20 прямоугольников)
у меня получилась программка,которая скомпилировалась, но не работает.знаю,что написала бред, но все же:
#include <stdio.h>
#include <math.h>
int main()
{
int i;
float dx, s, h, x, y;
dx=3/20;
s=0;
x=-1+dx*y;
for (i=1; i<21; i++)
{ y=x/((x*x+1)*(x*x+1));
s=s+dx*y;
x=x+dx;
i++;
}
printf("%f", s);
}

[KobolD]

        void Main(string[] args)
        {
            int n = 20;//Количество разбиений
            float a = -1;//Начало интервала
            float b = 2; //конец интервала
            
            float h = (b-a)/n;//Шаг
            double Integral=0;
            for (float x = a; x < b; x += h)
                Integral += Y(x);

        }
        double Y(float X)
        {
            return X / Math.Pow(Math.Pow(X, 2) + 1, 2);
        }

[iris_ka]

спасибо=))

[KobolD]

забыл домножить на шаг, чтобы площадь получилась, поправьте. Должно быть

Integral += Y(x)*h;

[iris_ka]

да, я с этим разобралась)