alpha  = linspace(0.01,0.99);

deltaS = alpha .* log(1./alpha) + (1-alpha) .* log(1./(1-alpha));

plot(alpha,deltaS);