#include <iostream>
using namespace std;

double* przejscie(double*, double, double, int);
double licz(double*, double*, double, double, int);
double f(double);

void Calka_Dwuwezlowa();
void Calka_Czterowezlowa();

int main() 
{
	Calka_Dwuwezlowa();
	Calka_Czterowezlowa();

	cout << "\n\n";
	system("pause");
	return 0;
}

void Calka_Dwuwezlowa()
{
	double a, b;
	cout << "Podaj przedzial a: "; cin >> a;
	cout << "Podaj przedzial b: "; cin >> b;

	int n1 = 2;

	double* tabDwuWezlowe = new double[n1];
	tabDwuWezlowe[0] = 0.577350;
	tabDwuWezlowe[1] = -0.577350;

	double* tabAk2 = new double[n1];
	tabAk2[0] = 1;
	tabAk2[1] = 1;

	double* tabT;
	tabT = przejscie(tabDwuWezlowe, a, b, n1);
	cout << "\nCalka dwuwezlowa: " << licz(tabT, tabAk2, a, b, n1);
}

void Calka_Czterowezlowa()
{
	double a, b;
	cout << "Podaj przedzial a: "; cin >> a;
	cout << "Podaj przedzial b: "; cin >> b;

	int n2 = 4;

	double* tabCzteroWezlowe = new double[n2];
	tabCzteroWezlowe[0] = -0.861136;
	tabCzteroWezlowe[1] = -0.339981;
	tabCzteroWezlowe[2] = 0.339981;
	tabCzteroWezlowe[3] = 0.861136;

	double* tabAk4 = new double[n2];
	tabAk4[0] = 0.347855;
	tabAk4[1] = 0.652145;
	tabAk4[2] = 0.652145;
	tabAk4[3] = 0.347855;

	double* tabT;
	tabT = przejscie(tabCzteroWezlowe, a, b, n2);
	cout << "\nCalka czterowezlowa: " << licz(tabT, tabAk4, a, b, n2);
}

double* przejscie(double* tab, double a, double b, int n) 
{
	double tmp;
	for (int i = 0; i<n; i++) 
	{
		tmp = (a + b) / 2 + (b - a) / 2 * tab[i];
		tab[i] = tmp;
	}

	return tab;
}

double licz(double* tab, double* Ak, double a, double b, int n) 
{
	double calka = 0;
	for (int i = 0; i<n; i++) 
		calka += Ak[i] * f(tab[i]);
	
	calka *= (b - a) / 2;
	
	return calka;
}

double f(double x) 
{
	return 1 / x;
}