/*  
  * Example of multidimensional array, returning an array result 
  * 
  */ 

#include <stdio.h> 
#define NROWS 5
#define NCOLS 5

void add_matrices(double a[][NCOLS], double b[][NCOLS], 
                  double c[][NCOLS], int nrows, int ncols); 

int 
main()
{ 
    double x[NROWS][NCOLS], y[NROWS][NCOLS], xpy[NROWS][NCOLS];
    int n,m,i,j; 

    printf("Enter nrows and ncols = "); 
    scanf("%d%d",&m, &n); 

    printf("Enter %d by %d matrix X = \n",m,n); 
    for (i=0; i < m; i++)
        for (j=0; j <n; j++) 
            scanf("%lf", &x[i][j]);   

    printf("Enter %d by %d matrix Y = \n",m,n); 
    for (i=0; i < m; i++)
        for (j=0; j <n; j++) 
            scanf("%lf", &y[i][j]);   
 
    add_matrices(x,y,xpy,m,n); 

    printf("the sum of X+Y = \n"); 
    for (i=0; i < m; i++){ 
        for (j=0; j <n; j++) 
            printf(" %.1f", xpy[i][j]);   
        printf("\n"); 
    } 

    return 0; 
} 

/*
 *  Adds corresponding elements of matrices a and b, 
 *  storing the result in c. 
 */
void
add_matrices(double a[][NCOLS], double b[][NCOLS], 
             double c[][NCOLS], int nrows, int ncols)
{
    int i, j;

    for (i = 0; i < nrows; i++)
        for (j = 0; j < ncols; j++)
             c[i][j] = a[i][j] + b[i][j];

}