1. #include <iostream>
2. #include <fstream>
3.  
4. void functionB(double **& A, double*& B, int n, double **& U, double **& L)
5. {
6.         //wiersz 1. U
7.         for (int j = 1; j <= n; j++)
8.                 U[1][j] = A[1][j];
9.  
10.         //kolumna 1. L
11.         for (int j = 2; j <= n; j++)
12.                 L[j][1] = A[j][1] / U[1][1];
13.        
14.         //>1
15.         for (int i = 1; i <= n; i++){
16.                 double sumU = 0;
17.                 double sumL = 0;
18.                 for (int j = i; j <= n; j++){           //U
19.                         for (int k = 1; k <= i - 1; k++)
20.                                 sumU = L[i][k] * U[k][j];
21.                         U[i][j] = A[i][j] - sumU;
22.                 }
23.                 for (int j = i + 1; j <= n; j++){       //L
24.                         for (int k = 1; k <= i - 1; k++)
25.                                 sumL = L[j][k] * U[k][i];
26.                         L[j][i] = (A[j][i] - sumL) / U[i][i];
27.                 }
28.         }
29. }
30.  
31. int main()
32. {
33.         int n;
34.         double **A = new double* [n + 1];
35.         double *B = new double[n + 1];
36.         double **U = new double*[n + 1];
37.         double **L = new double*[n + 1];
38.  
39.         for (int i = 0; i < n + 1; i++)
40.                 A[i] = new double[n + 1];
41.  
42.         for (int i = 0; i <= n; i++){
43.                 B[i] = 0;
44.                 for (int j = 0; j <= n; j++)
45.                 {
46.                         A[i][j] = 0;                            //zerowanie macierzy
47.                         L[i][j] = 0;
48.                         U[i][j] = 0;
49.                 }
50.         }
51.  
52.         for (int i = 0; i <= n; i++)            //1 na przekątnej L
53.                 L[i][i] = 1;                           
54. }
55.