1. #include <iostream>
2. #include <cmath>
3.  
4. using namespace std;
5. void print44(double A[4][4])
6. {
7.  
8.         for (int i=0;i<4;i++)
9.  
10.         {
11.                 for (int j=0;j<4;j++)
12.  
13.                 {
14.  
15.                         cout<<A[i][j]<<"     ";
16.  
17.  
18.                 }
19.                 cout<<endl;
20.  
21.         }
22.  
23.  
24.  
25.  
26. }
27. double sumau (int i, int j, double A[4][4])
28.  
29. {
30.         double suma=0;
31.         for (int k=0;k<(i-1);k++)
32.         {
33.  
34.                 suma=suma+A[i][k]+A[k][j];
35.  
36.         }
37.         return suma;
38.  
39.  
40. }
41.  
42.  
43. double sumal (int i, int j, double A[4][4])
44.  
45. {
46.         double suma=0;
47.         for (int k=0;k<(i-1);k++)
48.         {
49.  
50.                 suma=suma+A[j][k]+A[k][i];
51.  
52.         }
53.         return suma;
54.  
55.  
56. }
57. int main ()
58.  
59. {
60.  
61.         double A[4][4]={{3,-4,4,-4},
62.                                         {1.5,-1,2,-2},
63.                                         {1.5,-0.5,0,-3},
64.                                         {4.5,-5.5,4,-9}};
65.         print44(A);
66.         double u[4][4];
67.         double l[4][4];
68.         for (int i=0; i<4;i++)
69.         {
70.  
71.                 for (int j=0;j<4;j++)
72.                 {
73.                         if (i<=j) A[i][j]=A[i][j]-sumau(i,j,A);
74.                         if (i<j) A[j][i]=(A[j][i]-sumal(i,j,A))/A[i][i];
75.                 }
76.  
77.  
78.  
79.         }
80.         cout<<endl;
81.         print44(A);
82.  
83.  
84.         system("Pause");
85.         return 0;
86. }