#include <iostream>
#include <fstream>

void functionB(double **& A, double*& B, int n, double **& U, double **& L)
{
	//wiersz 1. U
	for (int j = 1; j <= n; j++)
		U[1][j] = A[1][j];

	//kolumna 1. L
	for (int j = 2; j <= n; j++)
		L[j][1] = A[j][1] / U[1][1];
	
	//>1
	for (int i = 1; i <= n; i++){
		double sumU = 0;
		double sumL = 0;
		for (int j = i; j <= n; j++){		//U
			for (int k = 1; k <= i - 1; k++)
				sumU = L[i][k] * U[k][j];
			U[i][j] = A[i][j] - sumU;
		}
		for (int j = i + 1; j <= n; j++){	//L
			for (int k = 1; k <= i - 1; k++)
				sumL = L[j][k] * U[k][i];
			L[j][i] = (A[j][i] - sumL) / U[i][i];
		}
	}
}

int main()
{
	int n;
	double **A = new double* [n + 1];
	double *B = new double[n + 1];
	double **U = new double*[n + 1];
	double **L = new double*[n + 1];

	for (int i = 0; i < n + 1; i++)
		A[i] = new double[n + 1];

	for (int i = 0; i <= n; i++){
		B[i] = 0;
		for (int j = 0; j <= n; j++)
		{
			A[i][j] = 0;				//zerowanie macierzy
			L[i][j] = 0;
			U[i][j] = 0;
		}
	}

	for (int i = 0; i <= n; i++)		//1 na przekątnej L
		L[i][i] = 1;				
}