#include <iostream.h>
#include <malloc.h>

void sum_matr(const int, const int, const int, double ***, double ***, double ***);

const int DIM1 =  10 ;
const int DIM2 =  15 ;
const int DIM3 =  20 ;

const int DIM11 =  3 ;
const int DIM22 =  4 ;
const int DIM33 =  5 ;

int main()
{
	double*** A = (double***)malloc(sizeof(double**)*DIM1);
	double*** B = (double***)malloc(sizeof(double**)*DIM1);
	double*** C = (double***)malloc(sizeof(double**)*DIM1);

	int i,j,k;

	for(i =  0 ; i < DIM1; i++)
	{
		A[i] = (double**)malloc(sizeof(double*)*DIM2);
		B[i] = (double**)malloc(sizeof(double*)*DIM2);
		C[i] = (double**)malloc(sizeof(double*)*DIM2);

		for(j=  0 ; j < DIM2; j++)
		{
			A[i][j] = (double*)malloc(sizeof(double)*DIM3);
			B[i][j] = (double*)malloc(sizeof(double)*DIM3);
			C[i][j] = (double*)malloc(sizeof(double)*DIM3);

			for(k= 0 ; k < DIM3; k++)
			{
				A[i][j][k] =  1 . 0 ;
				B[i][j][k] =  2 . 0 ;
				C[i][j][k] =  0 . 0 ;
			}
		}
	}

	double*** AA = (double***)malloc(sizeof(double**)*DIM11);
	double*** BB = (double***)malloc(sizeof(double**)*DIM11);
	double*** CC = (double***)malloc(sizeof(double**)*DIM11);

	for(i =  0 ; i < DIM11; i++)
	{
		AA[i] = (double**)malloc(sizeof(double*)*DIM22);
		BB[i] = (double**)malloc(sizeof(double*)*DIM22);
		CC[i] = (double**)malloc(sizeof(double*)*DIM22);

		for(j=  0 ; j < DIM22; j++)
		{
			AA[i][j] = (double*)malloc(sizeof(double)*DIM33);
			BB[i][j] = (double*)malloc(sizeof(double)*DIM33);
			CC[i][j] = (double*)malloc(sizeof(double)*DIM33);

			for(k= 0 ; k < DIM33; k++)
			{
				AA[i][j][k] =  1 . 0 ;
				BB[i][j][k] =  2 . 0 ;
				CC[i][j][k] =  0 . 0 ;
			}
		}
	}

	sum_matr(DIM1, DIM2, DIM3, A, B, C);
	sum_matr(DIM11, DIM22, DIM33, AA, BB, CC);

	for(i =  0 ; i < DIM1; i++)
		for(j=  0 ; j < DIM2; j++)
			for(k =  0 ; k < DIM3; k++)
				cout << "C[" << i << "][" << j << "][" << k << "] = " << C[i][j][k] << endl;

	for(i =  0 ; i < DIM11; i++)
		for(j=  0 ; j < DIM22; j++)
			for(k =  0 ; k < DIM33; k++)
				cout << "CC[" << i << "][" << j << "][" << k << "] = " << CC[i][j][k] << endl;

	for(i =  0 ; i < DIM1; i++)
	{
		for(j=  0 ; j < DIM2; j++)
		{
			free(A[i][j]);
			free(B[i][j]);
			free(C[i][j]);
		}

		free(A[i]);
		free(B[i]);
		free(C[i]);
	}

	free(A);
	free(B);
	free(C);

	for(i =  0 ; i < DIM11; i++)
	{
		for(j=  0 ; j < DIM22; j++)
		{
			free(AA[i][j]);
			free(BB[i][j]);
			free(CC[i][j]);
		}

		free(AA[i]);
		free(BB[i]);
		free(CC[i]);
	}

	free(AA);
	free(BB);
	free(CC);

	return  0 ;
}

void sum_matr(const int ind1, const int ind2, const int ind3, double ***a, double ***b, double ***c)
{
         for(int i =  0 ; i < ind1; i++)
              for(int j=  0 ; j < ind2; j++)
                   for(int k= 0 ; k < ind3; k++)
                        c[i][j][k] = a[i][j][k] + b[i][j][k];
}

SUBROUTINE MATRIX(N,M,K,A)
  INTEGER N,M,K
  DOUBLE PRECISION A( 1 )
  .....
  тело метода
END

DOUBLE PRESISION A( 100 , 200 , 300 )
CALL MATRIX( 100 , 200 , 300 ,A)