Помогите создать програму для вычисления определенного интеграла Методом Симпсона

[Сергей Человек]

Сижу уже второй день никак не могу понять......
задача заключается в седующем: Создать програму в среде С++.
интеграл от (xlnx)^2 dx c приделми от одного до експоненты(2.71).количесво отрезков = 52, з шагом (e-1)/8.
Первообразная от этой функции = (x^3)/27*(9ln^2x-6lnx+2)-2/27
надеюсь на вашу помощь...

[ОДИНОЧЕСТВО В СЕТИ]

// Метод Симпсона
//

#include "stdafx.h"
#include <iostream>
#include <cmath>


typedef double typeOfFloat;

// Вычисление функции f(x)
typeOfFloat f( typeOfFloat x );

// Вычисление определённого интеграла от f(x) на [a,b] разбиением на n отрезков
typeOfFloat integral( typeOfFloat a, typeOfFloat b, int n );

int _tmain(int argc, _TCHAR* argv[])
{

	using namespace std;

	typeOfFloat a,b;
	typeOfFloat E;
	typeOfFloat s;
	typeOfFloat prS = 0;
	int k;

	
	cin >> a;
	cin >> b;
	cin >> E;


	prS = integral(a,b,2);
	s = integral(a,b,4);

	
	for( k = 6; fabs( prS - s ) > E; k+=2 )
	{
		prS = s;

		s = 0;

		s = integral(a,b,k);

	}

	cout.precision(10);  // Точность это, если что... (кол-во цифр при выводе результата)
	cout << s << endl;
	cout << k;


	cin.get(); cin.get(); cin.get();
	return 0;
}


typeOfFloat f( typeOfFloat x ) 
{ 
	return (typeOfFloat)pow( typeOfFloat(4-3*x*x), typeOfFloat(1.0/3.0) ); // вместо 4-3*x*x --- ваша функция а вместо 1.0/3.0 --- ваши пределы
} 

typeOfFloat integral( typeOfFloat a, typeOfFloat b, int n )
{

	typeOfFloat x;
	typeOfFloat h;
	typeOfFloat sum1,sum2;


	h = (b-a)/n;


	sum1 = 0;
	sum2 = 0;

	
for( int i = 1; i < n; i++ )
{

	x = a+i*h;
	
	if( i % 2 == 0 )       //Сумма значений подынтегральной функции на концах четных внутренних отрезков
	{
		sum2 = sum2 + f(x);
	}
	else                   // Сумма значений подынтегральной функции на концах нечетных внутренних отрезков
	{

		sum1 = sum1 + f(x);
	}
}

	return h/3*(f(a)+f(b)+4*sum1+2*sum2);

}

[Сергей Человек]

выбрасывает две ошибки связаные с библиотеками((((что делать ???

[ОДИНОЧЕСТВО В СЕТИ]

stdafx.h--эту просит??

[profi]

Напишите так:

#include<iostream.h>
#include<math.h>

и выкиньте это

stdafx.h

Кстати каким компилятором пользуютесь?

[Сергей Человек]

Borland C 3.1

int _tmain(int argc, _TCHAR* argv[]) и еще здесь гдето ошибка((( просит гдето поставить скобку

ошибка в #include <cmath>

[profi]

#include <iostream.h>
#include <math.h>


typedef double typeOfFloat;

// Вычисление функции f(x)
typeOfFloat f( typeOfFloat x );

// Вычисление определённого интеграла от f(x) на [a,b] разбиением на n отрезков
typeOfFloat integral( typeOfFloat a, typeOfFloat b, int n );

int main()
{


	typeOfFloat a,b;
	typeOfFloat E;
	typeOfFloat s;
	typeOfFloat prS = 0;
	int k;

	
	cin >> a;
	cin >> b;
	cin >> E;


	prS = integral(a,b,2);
	s = integral(a,b,4);

	
	for( k = 6; fabs( prS - s ) > E; k+=2 )
	{
		prS = s;

		s = 0;

		s = integral(a,b,k);

	}

	cout.precision(10);  // Точность это, если что... (кол-во цифр при выводе результата)
	cout << s;
	cout << k;


	cin.get(); cin.get(); cin.get();
	return 0;
}


typeOfFloat f( typeOfFloat x ) 
{ 
	return (typeOfFloat)pow( typeOfFloat(4-3*x*x), typeOfFloat(1.0/3.0) ); // вместо 4-3*x*x --- ваша функция а вместо 1.0/3.0 --- ваши пределы
} 

typeOfFloat integral( typeOfFloat a, typeOfFloat b, int n )
{

	typeOfFloat x;
	typeOfFloat h;
	typeOfFloat sum1,sum2;


	h = (b-a)/n;


	sum1 = 0;
	sum2 = 0;

	
for( int i = 1; i < n; i++ )
{

	x = a+i*h;
	
	if( i % 2 == 0 )       //Сумма значений подынтегральной функции на концах четных внутренних отрезков
	{
		sum2 = sum2 + f(x);
	}
	else                   // Сумма значений подынтегральной функции на концах нечетных внутренних отрезков
	{

		sum1 = sum1 + f(x);
	}
}

	return h/3*(f(a)+f(b)+4*sum1+2*sum2);

}

[Сергей Человек]

ошибка вот здесь cout << s << end; нигде нет end

[profi]

Напиши так:
вместо

cout << s << end;

это

cout << s;

[Сергей Человек]

можешь еще немного разяснить что такое а b и E значения которые мы вводим