/* 
Каранин Алексей 08-ПМ
Лабораторная работа 14

  Матрицы

Задание: Дана квадратная матрица NxN. Заменить нулями все элементы матрицы, 
расположенные ниже главной и побочной диагоналей.
Исходную и полученную матрицы вывести на печать.
 */

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

void main()
{
	double **a;
	int n,i,j;
	FILE *fin, *fout;
	
	fin=fopen("matrix14.in","r");
	if (fin==NULL)
	{printf("Can't open file matrix14.in\n"); exit(1);}
	
	fout=fopen("matrix14.out","w");
	if (fout==NULL)
	{printf("Can't open file matrix14.out\n"); exit(1);}
	
	fscanf(fin,"%d",&n);

	a=(double **)malloc(n*sizeof(double *));
	if (a==NULL) {printf("Can't allocate memory\n");exit(1);}

	for (i=0; i<n; i++)
	{
		a[i]=(double *)malloc(n*sizeof(double));
		if (a[i]==NULL) {printf("Can't allocate memory\n");exit(1);}
	}

	for (i=0; i<n; i++)
		for (j=0; j<n; j++)
			fscanf(fin, "%lf", &a[i][j]);

	fprintf(fout, "\nИсходная матрица А размерности %d x %d\n",n,n);
	
	for (i=0; i<n; i++)
	{
	for (j=0; j<n; j++)
		fprintf(fout,"%5.1lf",a[i][j]);
	fprintf(fout,"\n");
	}

	for (i=1; i<n; i++)
	{
	for (j=0; j<i; j++) a[i][j] = 0;
	for (j=n-i; j<n; j++) a[i][j] = 0;
	}

    fprintf(fout, "\nРезультат\n");
	
	for (i=0; i<n; i++)
	{
	for (j=0; j<n; j++)
		fprintf(fout,"%5.1lf",a[i][j]);
	fprintf(fout,"\n");
	}


}