Doesn't really seem to match the manual. Maybe the manual is wrong? It states that
teff ~ int b^2 P(b) exp(-b) db
when the code calls for fsum to use a tau of 0.5, meaning that the exponent in the calculated integral is -b/2. So this is strange, please look into