/*vfun approximate form of Voit function */
double vfun(double x);