bmdplot function

BMD plot optionally with confidence intervals on BMD