mu=1:(-0.01):0.01;
x1 = (2*mu+1+sqrt(4*mu.^2+1))/2;
x2 = 1-mu./x1;
plot(x1,x2)
hold on
w1 = mu./x1;
w2 = mu./x2;
plot(w1,w2);
legend('x','w')