double inprod(doble a[],double b[],int c) {
 return c == 0 ? 0.0 : *a * *b + inprod(a + 1, b + 1, c - 1);
}