/* fichero antisim.c */

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

void main(void)
{
   double **a;
   double **sim;
   double **antisim;
   int    i, j, n;

   printf("Este programa calcula la descomposición en suma de simétrica\n");
   printf("y antisimétrica de una matriz cuadrada cualquiera.\n");
   printf("\n");
   printf("Introduzca la dimensión de la matriz a descomponer:\n");
   printf("Dimensión: ");
   scanf("%d", &n);
   printf("Introduzca los valores de la matriz:\n");

   a = malloc(n*sizeof(double*));
   sim = malloc(n*sizeof(double*));
   antisim = malloc(n*sizeof(double*));
   for (i=0; i<n; i++) {
      a[i] = malloc(n*sizeof(double));
      sim[i] = malloc(n*sizeof(double));
      antisim[i] = malloc(n*sizeof(double));
   }

   for (i=0; i<n ; i++)
      for ( j = 0; j < n; j++) {
         printf("%d,%d: ", i+1, j+1);
         scanf("%lf", &a[i][j]);
      }

   for (i=0; i<n ; i++)
      for (j=0; j<n ; j++) {
         sim[i][j] = (a[i][j] + a[j][i])*0.5;
         antisim[i][j] = (a[i][j] - a[j][i])*0.5;
      }

   printf("La matriz simétrica es:\n");
   for (i=0; i<n; i++) {
      for (j=0; j<n; j++)
         printf("%lf  ", sim[i][j]);
      printf("\n");
   }
   printf("\n");

   printf("La matriz antisimétrica es:\n");
   for (i=0; i<n; i++) {
      for (j=0; j<n; j++)
         printf("%lf   ", antisim[i][j]);
      printf("\n");
   }
   printf("\n");

   printf("La matriz original debe ser:\n");
   for (i=0; i<n; i++) {
      for (j=0; j<n; j++)
         printf("%lf   ", sim[i][j] + antisim[i][j] );
      printf("\n");
   }
   printf("\n");
}