Вычисление определенного интеграла (Си)

[Alex_Kardo]

Правильно ли составлена программа по вычислению определенного интеграла :

#include <stdio.h>
#include <stdlib.h> 
#include <conio.h>
#include <math.h> 
#include <windows.h> 

/*----Расчет значений функции---*/
double FUNCTION(double x)
{
	double y;
	if (x < -4) y = fabs(x*x - 3);
	else if (x <= 1.5) y = 3 * cos(x + 1);
	else y = (2/pow(x,4))+5;
	return y;
}

/*----Расчет интеграла по методу прямоугольников---*/
double RECTANGLES(double left, double right, double step)
{
	double sum = 0;
	double runner;
	for (runner = left + step * 0.5; runner < right; runner += step)
	sum += FUNCTION(runner);
	sum *= step;
	return sum;
}

int main(void)
{
	unsigned N;
	double x, x2, y, r, dx, integral;

	// Ввод исходных данных
	printf("Programm tabulatest function from x1 to x2, step dx\n");
	printf("Enter x1,x2,dx:\n");
	scanf_s("%lf%lf%lf", &x, &x2, &dx);
	printf ("Programm wait 5 second and clear screen every N*3 points. If you enter 0, the display will be displayed on one point. Enter N: ");
	scanf_s ("%u", &N);
	integral = RECTANGLES(x, x2, dx); // вызов функции расчета интеграла
	r = 1; //счетчик выведенных на экран точек
	system("cls");
	do
	{
		y = FUNCTION(x);
		printf ("f(%5.2lf)=(%5.5lf)\n", x, y);
		if (r > N * 3)
		{
			Sleep(5000); r = 1;
			system("cls");
		} else r++;
		x = x + dx;
	} while (x <= x2);
printf ("The value of the integral is: %lf\n", integral ); // Вывод приближенного значения определенного интеграла на экран
_getch();
}

P.S. Программа работает, но хочу быть уверен, что правильно. Заранее благодарю

[Stilet]

А почему на калькуляторе не проверишь?

[Alex_Kardo]

Меня интересует правильность алгоритма и если честно трудно это все дается, с математикой не все в порядке.