%-----------------------------------------
%script monster2
% plots the functions L1(x) = 2^x log(1+2^-x)
% (and L2(x) = 2^x log(1-2^-x) )
% For x>50 a dramatic deviation from the
% Theoretical value is found; for x>53 ( 54 for L2 )
% the value drops spuriously to 0 (WHY?)
x = linspace(50,55,333331);
y = 2.^(x).*log(1+2.^(-x));
z = 2.^(x).*log(1-2.^(-x));
plot(x,y,x,z)
title('plot of f1=2^xlog(1+2^-x), f2 = 2^xlog(1-2^-x)')