float matmat(int l, int u, int i, int j, float **a, float **b)
{
	int k;
	float s;

	s=0.0;
	for (k=l; k<=u; k++) s += a[i][k]*b[k][j];
	return (s);
}