View difference between Paste ID: qy6yrNRx and

-
+Program Sistema_Linear;
 
 	type vetor = array[1..20] of extended;
 		matriz = array[1..20,1..20] of extended;
 
 	var Constantes: vetor;
 		Mampliada: matriz;
 		Mauxiliar: matriz;
 		Cons_auxiliar: vetor;
 		i, j, N, Num_escolha: integer;
 		Escolha: char;
 
 function compara_vetor (Vetora, Vetorb: vetor): extended;
 	var Subtrai: vetor;
 		Max: extended;
 	begin
 		max := 0;
 		for i := 1 to N do
 			begin
 				Subtrai[i] := Vetora[i] - Vetorb[i];
 			end;
 		for i := 1 to N do
 			begin
 				if abs(Subtrai[i]) > Max then
 					begin
 						Max := abs(Subtrai[i]);
 					end;
 				compara_vetor := max;
 			end;
 	end;
 
 procedure escalona_sup (MatrizA: matriz; Vetorconst: vetor);
 
 	var i, j, k: integer;
 		aux: extended;
 	begin
 		for k := 1 to (N - 1) do
 			begin
 				for j := (k + 1) to N do
 					begin
 						aux := MatrizA[k, j]/MatrizA[j, j];
 						for j := (k + 1) to N do
 							begin
 								MatrizA[k, i] := MatrizA[k, i] - aux * MatrizA[j, i];
 							end;
 						Vetorconst[k] := Vetorconst[k] - aux * Vetorconst[j];
 					end;
 			end;
 	end;
 
 procedure escalona_inf (MatrizA: matriz; Vetorconst: vetor);
 	var i, k : integer;
 		s: extended;
 	begin
 		for k := 1 to (N - 1) do
 			begin
 				for i := (k + 1) to N do
 					begin
 						s := MatrizA[i, k]/MatrizA[k, k];
 						MatrizA[i, k] := 0;
 						for j := (k + 1) to N do
 							begin
 								MatrizA[i, j] := MatrizA[i, j] - s * MatrizA[k, j];
 							end;
 						Vetorconst[i] := Vetorconst[i] - s * Vetorconst[k];
 					end;
 			end;
 	end;
 
 procedure gauss (MatrizAmp: matriz; Constante: vetor);
 	var i, k: integer;
 		Resposta: vetor;
 		aux: extended;
 	begin
 		writeln(' ');
 		writeln('Gauss');
 		Mauxiliar := MatrizAmp;
 		Cons_auxiliar := Constante;
 		escalona_inf (Mauxiliar, Cons_auxiliar);
 		Resposta[N] := Cons_auxiliar[N]/Mauxiliar[N, N];
 		for i := (N - 1) downto 1 do
 			begin
 				aux := 0;
 				for j := (i + 1) to N do
 					begin
 						aux := aux + Mauxiliar[i, j] * Resposta[j];
 						Resposta[i] := (Mauxiliar[i] - aux)/Mauxiliar[i, i];
 					end;
 				for i := 1 to N do
 					begin
 						writeln('X', i, ' = ', Resposta[i]);
 						writeln(' ');
 					end;
 			end;
 	end;
 
 begin
 	repeat
 	clrscr;
 	write('Quantidade de variáveis = ');
 	readln(N);
 	for i := 1 to N do
 		begin
 			writeln(' ');
 			writeln('Fornecer os coeficientes da equacao ', i, ':');
 			for j := 1 to N do
 				begin
 					write('Coeficiente de X', j, ' = ');
 					readln(Mampliada[i,j]);
 				end;
 			write('Valor da constante da equacao ', i,' = ');
 			readln(Constantes[i]);
 	end;
 	writeln(' ');
 	writeln('Escolha o metodo de resolucao:');
 	writeln(' ');
 	writeln('1 - Eliminacao de Gauss');
 	writeln('2 - Metodo de eliminacao de Gauss-Jordan');
 	writeln('3 - Metodo de Jacobi');
 	writeln('4 - Metodo de Gauss-Seidel');
 	writeln(' ');
 	repeat
 	write(' Digite o numero do metodo escolhido (1-4): ');
 	readln(Num_escolha);
 	case num of
 		1: gauss(Mampliada, Constantes);
 		2: jordan(Mampliada, Constantes);
 		3: jacobi(Mampliada, Constantes);
 		4: seidel(Mampliada, Constantes);
 	end;
 	write('Escolher outro metodo? (S/N): ');
 	readln(Escolha);
 	Escolha := upcase(Escolha);
 	until Escolha = 'N';
 	write('Sair do programa? (S/N): ');
 	readln(Escolha);
 	resposta := upcase(Escolha);
 	until Escolha = 'S';
 end.