function [WBC]=WBC(x)
WBC=(1-gaussmf(x,[15 9])).*(x>9)+(1-gaussmf(x,[3.5 9])).*(x<9);