// fichero filcol.c 
// este programa crea una matriz, la lee de disco
//    y elimina la fila y la columna que indica el usuario
// en este programa se reserva memoria dinámica contigua
#include <stdio.h>
#include <stdlib.h>

void main(void)
{
	double **crear_matriz(int, int);
	void readmat(double **, int, int, FILE *);
	void imprimir_matriz(double **, int, int); 
	double **a, **b;
	int m, n, fil, col, i, j, ii, jj;
	char filename[51];
    FILE *pf;
	
	printf("Teclee el nombre del fichero de datos: ");
	scanf(" %[^\n]", filename);
	
	pf = fopen(filename,"r");
	fscanf(pf, "%d%d", &m, &n);
	a = crear_matriz(m, n);
	readmat(a, m, n, pf);
	b = crear_matriz(m-1, n-1);
	
	do {
		printf("\nTeclee la fila que desea eliminar (entre 1 y m): ");
		scanf("%d", &fil);
		printf("\nTeclee la columna que desea eliminar (entre 1 y n): ");
		scanf("%d", &col);
	} while ((fil<1 || fil>m) || (col<1 || col>n));
	
	for (i=0, ii=0; i<m; i++){
		if(i!=(fil-1)){
			for (j=0, jj=0; j<n; j++){
				if(j!=(col-1)){
					b[ii][jj] = a[i][j];
					jj++;
				}
			}
			ii++;
		}
	}
	printf("\nLa matriz original es:");
	imprimir_matriz(a, m, n);
	printf("\n\nLa matriz transformada es:");
	imprimir_matriz(b, m-1, n-1);
    printf("\nYa he terminado. m=%d, n=%d", m, n);
    fclose(pf);
}
  
double **crear_matriz(int m, int n)
{
	double **a;
	int i;
	a = calloc(m, sizeof(double*));
	a[0] = calloc(m*n, sizeof(double));
	for(i=1; i<m; i++) 
		a[i] = a[i-1]+n;
	return a;
}

void readmat(double **a, int m, int n, FILE *pf)
{
	int i, j;
	
	for (i=0; i<m; i++)
		for (j=0; j<n; j++)
			fscanf(pf, "%lf", &a[i][j]);
	return;
}

void imprimir_matriz(double **a, int m, int n)
{
	int i, j;
	for (i=0; i<m; i++) {
		printf("\n");
		for (j=0; j<n; j++)
			printf("  %lf", a[i][j]);
	}
	return;
}