function [K]=K(x)
        
K=(1-gaussmf(x,[1 4.5])).*(x>4.5)+(1-gaussmf(x,[1.5 4.5])).*(x<4.5);