void trmul(a,b,m,n,k,c)
int m,n,k;
double a[],b[],c[];
{ 
	int i,j,l,u;
    for (i=0; i<=m-1; i++)
		for (j=0; j<=k-1; j++)
		{ 
			u=i*k+j; c[u]=0.0;
			for (l=0; l<=n-1; l++)
				c[u]=c[u]+a[i*n+l]*b[l*k+j];
		}
    return;
}