вычисление интеграла Borland C++

[klykovka]

Даны вещественные числа a, b. Значения функции f(x)=(x+1)2√lg x; a=2; b=10;
записать в массив. Вычислить значение интеграла, используя:
1) Формулу трапеций

I1=h*[f(a)/2+f(a+h)+f(a+2h)+…+f(a+(n-1)h)+f(b)/2]


2) Формулу Симпсона

I2=h/3*(f(a)+f(b)+4*(f(a+h)+f(a+3h)+…+f( a+(n-1)h))+

2*(f(a+2h)+f(a+4h)+…+f(a+(n-2)h)))


h=(b-a)/n, n=100.

Вот такие соображения по алгоритму: расчёт выполняем по формулам (формулы даны) .Объявим массив func[n+1], где n - количество отсчётов функции на интервале от а до b , запмсываем знач-я функции в массив, (шаг вычисляем (b-a)/n ,аргумент ф-ии a+i*h, где i параметр цикла, меняющийся от 0 до n ( т. е. всего n значений). массив используем для вычислений по обеим формулам. в первой формуле сумма знач-й массива от 1 до n-1, плюс граничные func[0]/2+ func[n]/2 и умножаем на шаг(h). Во второй форм. две суммы - нечётные и чётные значения сумм-ся отдельно и потом по формуле находится общая сумма.
Помогите с кодом.

[strider]

#include <stdio.h>
#include <math.h>

//number of intervals
#define N 1000

const double a = 2.0, b = 10.0;
const double h = (b-a)/N;

double func[N+1];

double f(double x)
{
	return (x+1)*2*sqrt(log10(x));
}

double trapecii()
{
	double result = 0.0;
	for(int i=0; i<N; ++i)
		result += (func[i] + func[i+1]) * h/2;
	return result;
}

double simpson()
{
	double result = 0.0;
	result += (func[0] + func[N])*h/3;
	for(int i=1; i<N; i+=2)
		result += 4*h/3*func[i];
	for(int i=2; i<N; i+=2)
		result += 2*h/3*func[i];

	return result;
}

void main()
{
	for(int i=0; i<=N; ++i)
		func[i] = f(a+i*h);


	double I1, I2;
	I1 = trapecii();
	I2 = simpson();

	printf("%lf\n%lf\n", I1, I2);
}

[klykovka]

for(int i=2; i<N; i+=2) выдает ошибку в этой строке:Multiple declaration for
может int убрать?

[Stilet]

Проверил код на VS2010 - Никаких ошибок нет.
Программа работает и выдает какой то результат.
Единственное что добавить:

	printf("%lf\n%lf\n", I1, I2);
 getchar();
}

Чтоб экран ответа не уходил сразу.

[klykovka]

Спасибо все разобрался

Единственное что добавить:
[CODE] printf("%lf\n%lf\n", I1, I2);
getchar();
}
ну эт понятно. у меня была ошибка скорее из-за того что переменные объявлялись как вещественные а потом обозначались как целые. но щас работает все исправил