Dastur tanasi

Masalaga ko’ra tenglamalar sistemasini Gauss usulida yechish uchun noma’lumlar oldidagi koeffitsientlardan hosil bo’lgan massiv diagonali quyi qismini 0 larga aylanitirishimiz kerak. Buni men qism dastur sifatida yozdim. Bunda esa yuqorida qayd qilingan operatorlar bilan birga FOR takrorlanish operatorini ham qo’lladim. Bu operator ma’lum o’zgaruvchi ma’lum qiymatga yetguncha tarkorlanishlarni bajaradigan operator hisoblanadi. Bunga misolni siz asosiy dasturda ko’rishingiz mumkin.
Dastur tanasi.
Avvalo har qanday dasturni tuzishdan oldin qilinishi kerak bo’lgan ishlar algoritmi tuzib olinadi.
1 - chizma

Yuqoridagi blok sxema dasturni ishga tushirish uchun zarur buyruqlar tartibidir. Lekin quyidagi programmada esa buyruqlar blok sexmadagidan ko’p. chunki ko’p buyruqlar porgramma naijasini ko’rishdagi dizayn, natijani chiqarishda qulaylik yaratish uchun yozilgan.

#include

#include
#include

using namespace std;

void print(vector< vector > A) {
int n = A.size();
for (int i=0; ifor (int j=0; jcout << A[i][j] << "\t";
if (j == n-1) {
cout << "| ";
}
}
cout << "\n";
}
cout << endl;
}

vector gauss(vector< vector > A) {

int n = A.size();

for (int i=0; i
// Search for maximum in this column
double maxEl = abs(A[i][i]);
int maxRow = i;
for (int k=i+1; kif (abs(A[k][i]) > maxEl) {
maxEl = abs(A[k][i]);
maxRow = k;
}
}

// Swap maximum row with current row (column by column)

for (int k=i; kdouble tmp = A[maxRow][k];
A[maxRow][k] = A[i][k];
A[i][k] = tmp;
}

// Make all rows below this one 0 in current column

for (int k=i+1; kdouble c = -A[k][i]/A[i][i];
for (int j=i; jif (i==j) {
A[k][j] = 0;
} else {
A[k][j] += c * A[i][j];
}
}
}
}

// Solve equation Ax=b for an upper triangular matrix A

vector x(n);
for (int i=n-1; i>=0; i--) {
x[i] = A[i][n]/A[i][i];
for (int k=i-1;k>=0; k--) {
A[k][n] -= A[k][i] * x[i];
}
}
return x;
}

int main() {

int n;
cin >> n;

vector line(n+1,0);

vector< vector > A(n,line);

// Read input data

for (int i=0; ifor (int j=0; jcin >> A[i][j];
}
}

for (int i=0; i
cin >> A[i][n];
}

// Print input

print(A);

// Calculate solution

vector x(n);
x = gauss(A);

// Print result

cout << "Result:\t";
for (int i=0; icout << x[i] << " ";
}
cout << endl;
}
7 - chizma.

8 - chizma


Shunday qilib, oldimizga qo’yilgan masala o’z yechimini topdi. Lekin natijalar tahliliga keladigan bo’lsak, so’ralgan Z qiymat har doim har hil qiymatda chiqadi. Chunki bizga A va B matritsalar qiyamti ihtiyoriy ekani ma’lum qilingan. Uning bunday katta qiymat olganiga to’xtaladigan bo’lsak, x