Estudiante

#include <stdio.h>
#include <conio.h>
#include <stdlib.h>

double** LeerDatos(char *nombrefichero,int *dimension)
{

double **matriz;
int i,j;
FILE *f;if(f=fopen(nombrefichero,"r"),!f){
printf("Problema con la lectura del fichero\n");
return matriz;
}
fscanf(f,"%d\n",&(*dimension));

/* Asignacion de memoria dinamica*/matriz=(double**)malloc(*dimension*sizeof(double*));
for(i=0;i<*dimension;i++)
matriz[i]=(double*)malloc((*dimension+1)*sizeof(double));

for(i=0;i<*dimension;i++)
for(j=0;j<*dimension;j++)fscanf(f,"%lf\n",&(matriz[i][j]));

for(i=0;i<*dimension;i++)
fscanf(f,"%lf\n",&(matriz[i][*dimension]));

fclose(f);
return matriz;
}

void liberaMatriz(double **matriz,int dimension){
int i;for(i=0;i<dimension;i++)
free(matriz[i]);
free(matriz);
}
double **multiplica(double **A,double **B,int dimension){
int i,j,k;
double **C;
C=(double**)malloc(dimension*sizeof(double*));for(i=0;i<dimension;i++)
C[i]=(double*)malloc((dimension+1)*sizeof(double));

for (i=0;i<dimension;i++)
for (j=0;j<dimension;j++)
for (k=0;k<dimension;k++)C[i][j]=C[i][j]+A[i][k]*B[k][j];
return C;
}
double *matrizPorVector(double **matriz,double *vector,int dimension){
double *resultado,mult;
int i,j;

resultado=(double*)malloc(dimension*sizeof(double));for(i=0;i<dimension;i++){
mult=0;
for(j=0;j<dimension;j++)
mult+=(matriz[i][j]*vector[j]);
resultado[i]=mult;
}

return resultado;
}
double *A,B;
double *suma_vector(double *A,double *B,intdimension){
int i;
double *C;

C=(double*)malloc(dimension*sizeof(double));

for (i=0;i<dimension;i++)
C[i] = A[i]+B[i];

return C;
}

void jacobi2 (double **m, double *v, double*aprox, int dimension, int iter){
double **D1,**LR, *C;
int i,j;

D1=(double**)malloc(dimension*sizeof(double*));
for(i=0;i<dimension;i++)
D1[i]=(double*)malloc((dimension)*sizeof(double));...