void elmvec(int l, int u, int shift, float a[], float b[], float x)
{
	for (; l<=u; l++) a[l] += b[l+shift]*x;
}