function [DiasABP]=DiasABP(x)
    
DiasABP=((1-gaussmf(x,[50 154.5])).*(x>154.5)+(1-gaussmf(x,[35 154.5])).*(x<154.5));