function [Mg]=Mg(x)
Mg=(1-gaussmf(x,[2.5 2.1])).*(x>2.1)+(1-gaussmf(x,[.8 2.1])).*(x<2.1);