void inivec(int l, int u, float a[], float x)
{
	for (; l<=u; l++) a[l]=x;
}