Есть ли скрипт?

Question

[vaselekk]

Пробовал делать сам но вышло не очень

function celoe(val) {
		return (Math.ceil(val * 100) / 100)
	}



var m = [
[1, 3, 2,     5],
[2, 1, 3,     4],
[3, 2, 4,     6],
]


m =	redaction_matrica(m)

	function redaction_matrica(m1) {
		var v1 = []
		var v0 = []
		for (var i of m1) {
			if (i[0] !== 0) {
				v1.push(i)
			}
			if (i[0] == 0) {
				v0.push(i)
			}
		}
		for (var b of v0) {
			v1.push(b)
		}
		
		return v1
	}

	function otnimanie_vectorov(a, b, coeff) {
		var mass = []
		for (var i = 0; i < a.length; i++) {
			mass.push(a[i] - b[i] * coeff)
		}
		return mass;
	}
	m[1] = otnimanie_vectorov(m[1], m[0], m[1][0] / m[0][0]);
	m[2] = otnimanie_vectorov(m[2], m[0], m[2][0] / m[0][0]);
	/* 
	Матрица на текущем этапе должна выглядить так:
    0   0    0
	0   a1   a2
	0   a3   a4 
	*/
	var new_matrica = [
		[m[1][1], m[1][2], m[1][3]],
		[m[2][1], m[2][2], m[2][3]]
	]
	new_matrica = redaction_matrica(new_matrica)


	new_matrica[1] = otnimanie_vectorov(new_matrica[1], new_matrica[0], new_matrica[1][0] / new_matrica[0][0]);

// a*t + a2*v + a3*u = b

	var t = celoe((m[0][3] - (m[0][2] * u) - (m[0][1] * v)) / m[0][0])
	var v = celoe((new_matrica[0][2] - (new_matrica[0][1] * u)) / new_matrica[0][0])
    var u = celoe(new_matrica[1][2] / new_matrica[1][1])


console.dir(t)
console.dir(v)
console.dir(u)

Есть ли скрипт для нахождения неизвестных по методу гауса, для систем уравнений с 3 неизвестными?

Comments

[WbICHA]

Есть

Но на плюсах.

#include <iostream>
#include <stdint.h>
#include <iomanip>
using namespace std;

void showMatrix(double **A, double *B, int n) {
	cout << endl << "==============================================" << endl;

	for (int j = 0; j < n; j++) {
		for (int i = 0; i < n; i++) {
			cout << setw(10) << A[j][i] << " ";
		}
		cout << " = " << setw(10) << B[j] << endl;
	}
}

int main() {
	int n;
	cout << "Enter matrix N*N: ";
	cin >> n;

	//matrix A input
	cout << "Enter matrix elements: \n";
	double** A = new double*[n];
	for (int i = 0; i < n; i++) {
		A[i] = new double[n];
	}
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			cout << "  A [" << i + 1 << "][" << j + 1 << "] = ";
			cin >> A[i][j];
		}
		cout << endl;
	}

	//matrix B input
	cout << "Enter free elements:\n";
	double* B = new double[n];
	for (int i = 0; i < n; i++) {
		cout << "  B [" << "|" << i + 1 << "|" << "] = ";
		cin >> B[i];
	}

	for (int i = 0; i < n - 1; i++) {
		showMatrix(A, B, n);

		int m = i;
		for (int j = i + 1; j < n; j++) {
			if (labs(A[m][i]) < labs(A[j][i])) 
				m = j;
		}
		if (m > i) {
			for (int k = i; k < n; k++) { // k = i ─ не перемещаем "обнулённые" значения; k = 0 ─ перемещаем;
				double c = A[m][k];
				A[m][k] = A[i][k];
				A[i][k] = c;
			}

			double c = B[m];
			B[m] = B[i];
			B[i] = c;
		}

		for (int j = i + 1; j < n; j++) {
			double c = -A[j][i] / A[i][i];
			for (int k = i; k < n; k++) { // k = i + 1 ─ не обнуляем; k = i ─ обнуляет;
				A[j][k] += c * A[i][k];
			}
			B[j] += c * B[i];
		}
	}


	cout << endl << "=== FINAL ===";
	showMatrix(A, B, n);


	double* X = new double[n];

	int m = n - 1;
	X[m] = B[m] / A[m][m];
	for (int i = m - 1; i >= 0; i--) {
		for (int k = i + 1; k < n; k++) {
			B[i] -= X[k] * A[i][k];
		}
		X[i] = B[i] / A[i][i];
	}

	cout << endl;
	for (int i = 0; i < n; i++) {
		cout << "x [" << i + 1 << "] = " << X[i] << "; " << endl;
	}

	cout << endl;
	system("PAUSE");
	return 0;
}

Код писался за один вечер, чтобы на следующий день сдать и забыть о нём, так что код сильно пахнет.)

[Роман]

https://mathjs.org/docs/datatypes/matrices.html