clear all;
close all;

p=0:0.01:1;
s=length(p);
randn(seed,13);
r=0.01*randn(1,s);
q=p+r;
q=(q-min(q))/(max(q)-min(q));

x=-p.*log2(p);
y=-q.*log2(q);

figure;plot(x);grid on;
figure;plot(x);grid on;
hold on;plot(y,color,red);