#include float boost(plft) float plft; { if(plft < .5) plft = 1. - plft; return((float)sqrt((double)(1./plft))); }