bmdplotwithgradient function

BMD plot with color gradient