function [BUN]=BUN(x)
    
BUN=(1-gaussmf(x,[55 15.26])).*(x>15.26)+(1-gaussmf(x,[13 15.26])).*(x<15.26);