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