function [HCT]=HCT(x)
    
HCT=((1-gaussmf(x,[12 38])).*(x>38)+(1-gaussmf(x,[10 38])).*(x<38));