/************************************************/
/*     Function Multiplica						*/
/************************************************/

void Multiplica(COMPLEX **Who,COMPLEX **Global,int N)
{
	int j, i,k;
	COMPLEX **Y;
Y = getComplexMatrix(N);

  for(i=0;i<N;++i){
	   for(j=0;j<N;++j){
	     for(k=0;k<N;++k){

Y[i][j].real +=  ((Who[i][k].real * Global[k][j].real) - (Who[i][k].imag * Global[k][j].imag));
Y[i][j].imag += ((Who[i][k].real * Global[k][j].imag) + (Who[i][k].imag * Global[k][j].real));
		}

	}
}

for(i=0;i<N;++i){
for(j=0;j<N;++j){
	Global[i][j].real=Y[i][j].real;
	Global[i][j].imag=Y[i][j].imag;
  }
  printf("\n");
  }
  free(Y);
  free(Who);


}