FIXED BUG in Freedman-Diaconis bin calculation (was dividing the correct value into the range) and FIXED BUG in quartile method.  I validated both quartile methods now implemented, tukey and moore-mccabe, and I double checked sturges, scott and fd for accuracy