t=linspace(-.2,.2,1000);
q=-.5*(sqrt(1/2)+t).^2 +...
                (1/2).*(sqrt(1/2)+t)-(1/2)*sqrt(1-(sqrt(1/2)+t).^2);
plot(q)
sym x
lam1=-2;
lam2=0;
mu=1;
b1=-(lam1+mu)*x;
b2=-(lam2+mu)*sqrt(1-x);

