clear all;
close all;

p=0.01:0.01:0.5;
x_en=-p.*log2(p) - (1-p).*log2(1-p);
figure;plot(p,x_en);grid on;

q1=0.01:0.01:0.99;
p1=ones(1,length(q1));

x_en1=-p1.*log2(q1) - (1-p1).*log2(1-q1);
figure;plot(q1,x_en1);
hold on;plot(0.5,1,o);grid on;