**fig10_38.txt**

/* 1*/	// Standard matrix multiplication.
/* 2*/	// Arrays start at zero.

/* 3*/	void
/* 4*/	Matrix_Multiply( Matrix A, Matrix B, Matrix C, unsigned int N )
/* 5*/	{
/* 6*/	    int i, j, k;

/* 7*/	    for( i = 0; i < N; i++ )		// Initialization.
/* 8*/	        for( j = 0; j < N; j++ )
/* 9*/	            C[ i ][ j ] = 0.0;

/*10*/	    for(  i = 0; i < N; i++ )
/*11*/	        for(  j = 0; j < N; j++ )
/*12*/	            for( k = 0; k < N; k++ )
/*13*/	                C[ i ][ j ] += A[ i ][ k ] * B[ k ][ j ];
/*14*/	}