ТАШКЕНТСКИЙ УНИВЕРСИТЕТ
ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ ИМЕНИ
МУХАММАДА АЛ-ХОРАЗМИЙ
Предмет: «Линейная алгебра»
Самостоятельная работа №1
Вариант-78
Выполнил студент факультета
« Праграмный инжинеринг»
Проверил:Абдуллаев Рустам Заирович
Ташкент-2022 г.
#include
#include
using namespace std;
int main(){
long double A[5][5], S[5][5]={0}, B[5], X[5], Y[5];long double V; cout << "V ni kiriting jurnaldagi nomerini = "; cin >> V;
B[0] = 16 * V * V * V + 87 * V * V + 163 * V + 77;
B[1] = V * V * V + 8 * V * V + 61 * V + 40;
B[2] = 14 * V * V * V + 86 * V * V + 90 * V + 31;
B[3] = 18 * V * V * V + 46 * V * V + 93 * V + 57;
B[4] = 19 * V * V * V + 115 * V * V + 125 * V + 56;
//AT*A
A[0][0] = 5 * V * V + 20 * V + 30; A[0][1] = 2 * V * V + 10 * V + 25; A[0][2] = V * V + 6 * V - 5; A[0][3] = 6 * V * V + 11 * V + 8; A[0][4] = 2 * V * V + 7 * V - 18;
A[1][0] = 2 * V * V + 10 * V + 25; A[1][1] = 8 * V * V + 24 * V + 46; A[1][2] = -8 * V * V - 10 * V - 7; A[1][3] = 9 * V * V + 24 * V + 17; A[1][4] = -10 * V * V - 13 * V - 28;
A[2][0] = V * V + 6 * V - 5; A[2][1] = -8 * V * V - 10 * V - 7; A[2][2] = 13 * V * V + 8 * V + 6; A[2][3] = -8 * V * V + 3 ; A[2][4] = 16 * V * V + 15 * V + 15;
A[3][0] = 6 * V * V + 11 * V + 8; A[3][1] = 9 * V * V + 24 * V + 17; A[3][2] = -8 * V * V + 3 * V; A[3][3] = 20 * V * V + 12 * V + 14; A[3][4] = -9 * V * V - 3 * V + 1;
A[4][0] = 2 * V * V + 7 * V - 18; A[4][1] = -10 * V * V - 13 * V - 28; A[4][2] = 16 * V * V + 15 * V + 15; A[4][3] = -9 * V * V - 3 * V + 1; A[4][4] = 20 * V * V + 22 * V + 41;
cout<<"simmetrik matritsa: \n";
for(int i=0;i<=4;i++){
for(int j=0;j<=4;j++){
cout<}
cout<<'\n';
}
cout<<"------------------------------------\n";
//find S
int i, k, j; long double L;
for (i = 0; i <= 4; i++) {
for (k = 0; k <= 4; k++) {
if (i <= k) {
if (i == k) {
L = 0;
for (j = 0; j <= i - 1; j++) {
L += S[j][i] * S[j][i];
}
S[i][i] = pow(A[i][i] - L , 1/2.0);
}
else {
L = 0;
for (j = 0; j <= i - 1; j++) {
L += S[j][i] * S[j][k];
}
S[i][k] = (A[i][k] - L) / S[i][i];
}
}
}
}
//find Y
for (int a = 0; a <= 4; a++) {
L = 0;
for (int b = 0; b < a; b++) {
L += S[b][a] * Y[b];
}
Y[a] = (B[a] - L) / S[a][a];
}
//find X
for (int a = 4; a >= 0; a--){
L = 0;
for (int b = 4; b > a; b--) {
L += S[a][b] * X[b];
}
X[a] = (Y[a] - L) / S[a][a];
}
//cout
for (k = 0; k <5; k++) {
cout <<"X["<}
cout<<"\nA matritsa determinanti: "<return 0;
}
Do'stlaringiz bilan baham: |