void simplePVect(double* ivec, int* len)
{
  int i;

  for(i=0; i < *len; i++)
    printf("Element %d is %f \n", i, ivec[i]);

}