GaussSeidel.cpp
1dc59d732f73master

GaussSeidel.cpp
View Options

			#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;
			}