Metodo de gaus jacobi

/*Método de Gauss-Jacobi:
*/

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>
#include<math.h>
#include <cstdlib>
#include <iostream>

int main(){ int i, j, size, cont=0, n=0; double e, maior=0;; // O usuario deve informar o tamanho da matriz. printf ("\nDigite o tamanho da matriz 1, lembrando que a mesma deve ser quadratica: "); scanf ("%d", &size); // o Usuario deve informar o erro. printf ("\nDigite o criterio de parada, ou seja o erro que se deve atingir: "); scanf ("%lf", &e); double matA[size][size], matB[size], matX0[size], matC[size][size], matG[size], matXN[size]; double matX[size], D1[size], DR, D1MAIOR, XMAIOR, matXpositivo[size],soma[n]; printf("\n\nDigite a matriz A, a matriz B e o chute inicial\n\n"); do{ if (cont == 0){ for (i=0; i<size; i++){ for (j=0; j<size; j++){ //O usuario ira informar a matriz A. printf ("\na[%d][%d]: ", i+1, j+1); scanf ("%lf", &matA[i][j]); } //O usuario ira informar a matriz B. printf ("\nb[%d]: ", i+1); scanf ("%lf", &matB[i]);

//O usuario deverá informar o chute inicial. printf ("\nc[%d]: ", i+1); scanf ("%lf", &matX0[i]); } } else{ for (i=0; i<size; i++) matX0[i] = matX[i]; }

for (i=0; i<size; i++){ for (j=0; j<size; j++){ if (i == j){ //A diagonal principal da matriz C deve ser nula. matC[i][j] = 0; } else{ //Calculo dos elementos da matriz,