for a=0:0.001:0.125
plot(a,sqrt(0.5*(1-2*a+sqrt(1-8*a))),'b.')
hold on
plot(a,sqrt(0.5*(1-2*a-sqrt(1-8*a))),'b.')
end

xlabel('a', 'FontSize',30)
ylabel('b', 'FontSize',30)