function y=mu1mu2(t)
y=(1-1./(1+exp(-14*(t-pi/8))))./(1+exp(-14*(t+pi/8)));
