******** fig10.38 **********

	/* Standard matrix multiplication.  Arrays start at 0 */

	void
	matrix_multiply( matrix A, matrix B, matrix C, unsigned int n )
	{
		int i, j, k;

		for( i=0; i<n; i++ )	/* Initialization */
			for( j=0; j<n; j++ )
				C[i][j] = 0.0;
		for( i=0; i<n; i++ )
			for( j=0; j<n; j++ )
				for( k=0; k<n; k++ )
					C[i][j] += A[i][k] * B[k][j];
	}