GaussSeidel.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include <iostream>
	#include "GaussSeidel.h"

	using namespace std;

	int main() {
	int size_system = 5;
	Vector x(size_system);
	for(int j=1;j<=size_system;++j)
	{
	x(j)=1;
	}
	Matrix A(size_system,size_system);

	cout << "La matrice A vaut : " << endl;
	for(int i=1;i<=size_system;++i) {
	for (int j = 1; j <= size_system; ++j) {
	if(i==j) {
	A(i,j)=25;
	cout << A(i,j) << " ";
	}else{
	A(i,j) = 0;
	cout << A(i,j) << " ";
	}
	}
	cout << endl;
	}
	cout << endl;

	Vector b(size_system);
	int l=length(b);

	for(int j=1;j<=l;++j)
	{
	b(j)=2;
	}

	GaussSeidel linearSystem(A,b);

	cout << endl;
	TerminalPrint(linearSystem.Solve(),"La valeur de x est: ");

	return 0;
	}