/*  harjoitus 9 - tehtävä 2 */

#include <stdio.h>

void print_matrix(double a[][3]);
void msum(double a[][3], double b[][3], double sum[][3]);
void mdiff(double a[][3], double b[][3], double diff[][3]);


void main()
{
  double a[3][3] = {{ 1, 2, 3},
		    {-1, 0, 1},
		    { 3, 1, 0}};
  double b[3][3] = {{ 2, 1, 0},
		    { 4,-7, 1},
		    { 0, 4, 1}};
  double c[3][3];

  printf("a = \n");
  print_matrix(a);
  printf("b = \n");
  print_matrix(b);

  msum(a, b, c);
  printf("sum = \n");
  print_matrix(c);

  mdiff(a, b, c);
  printf("diff = \n");
  print_matrix(c);
}

void print_matrix(double a[][3])
{
  int i, j;
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      printf(" %8.3lf", a[i][j]);
    }
    printf("\n");
  }
}

void msum(double a[][3], double b[][3], double sum[][3])
{
  int i, j;
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      sum[i][j] = a[i][j] + b[i][j];
    }
  }
}

void mdiff(double a[][3], double b[][3], double diff[][3])
{
  int i, j;
  for (i = 0; i < 3; i++) {
    for (j = 0; j < 3; j++) {
      diff[i][j] = a[i][j] - b[i][j];
    }
  }
}