Формула Гаусса-Лежандра с++

[zubber94]

Здравствуйте! Если кто-то знает, что такое расчет определенного интеграла по формуле Гаусса-Лежандра для случая четырех узлов, проверьте данный код на правильность)

ddouble f(double x)
{
	return sin(x);
}

void output(double);

void solution(float a, float b, double arr[], int n)
{
	double I=0;

	for(int i=0; i<4; i++)

	I+=((b-a)/2)*(f((b+a)/2+(b-a)*arr[i]/2));
	output(I);
}

Переделал для for.

[Smitt&Wesson]

А зделать n через for, религия не позволяет?