private double exptLessThan(double b) {
    	double sum=0,interval=(b-12)/1000,val=12;
		for (int i = 0; i <1000; i++) {
			sum+=phi(val, mu, sigma)*val*interval;
			val+=interval;
		}
		return sum/probLessThan(b);
	}