float tamvec(int l, int u, int i, float **a, float b[])
{
	int k;
	float s;

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