#include <stdio.h>
int DIM1  =  2 ;
int DIM2 =  3 ;
int DIM3 =  4 ;
void sum_matr(int DIM1,int DIM2, int DIM3, double ***a, double ***b, double ***c)
{
         for(int i =  0 ; i < DIM1; i++)
              for(int j=  0 ; j < DIM2; j++)
                   for(int k= 0 ; k < DIM3; k++)
                        c[i][j][k] = a[i][j][k] + b[i][j][k];
};
int main()
{
    int i,j,k;
    double a[DIM1][DIM2][DIM3];
    double b[DIM1][DIM2][DIM3];
    double c[DIM1][DIM2][DIM3];
    for (i= 0 ; i < DIM1; i++) 
	for(j= 0 ;j<DIM2; j++) 
	    for (k= 0 ; k < DIM3; k++) 
    {
         a[i][j][k] = i+j + k;
         b[i][j][k] = i-j - k;
    }
    for (i= 0 ; i < DIM1; i++) 
    {
	for(j= 0 ;j<DIM2; j++) 
        {
	    for (k= 0 ; k < DIM3; k++) 
    	    {
                printf("a[%d][%d][%d]=%f; ", i, j, k, a[i][j][k]);
            }
	    printf("\n");
        }
        printf("\n");
    };
    for (i= 0 ; i < DIM1; i++) 
    {
	for(j= 0 ;j<DIM2; j++) 
        {
	    for (k= 0 ; k < DIM3; k++) 
    	    {
                printf("b[%d][%d][%d]=%f; ", i, j, k, b[i][j][k]);
            }
	    printf("\n");
        }
        printf("\n");
    };
    sum_matr(DIM1,DIM2,DIM3,a,b,c);
    for (i= 0 ; i < DIM1; i++) 
    {
	for(j= 0 ;j<DIM2; j++) 
        {
	    for (k= 0 ; k < DIM3; k++) 
    	    {
                printf("a[%d][%d][%d]=%f; ", i, j, k, c[i][j][k]);
            }
	    printf("\n");
        }
        printf("\n");
    };
    return  0 ;
};