%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%
%%%%%  Capacity to the Morse code
%%%%%     taudot}=1
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
p=.1:.00001:.9; 
taudot=1;
y=2*p+2*p.^3+p.^7-1;
[yzero,ind]=min(abs(y));
pdot=p(ind)
K=-log(pdot)/log(2)/taudot
