private static double ip(double x) {
	        /* Purpose: return the integral part of a double value.     */
	        //  double  tmp;
	        // modf(x, &tmp);
	        //return tmp;
	        double result = round(x);
	        return result;
	    }