%-----------------------------------------
%script monster1
% plots the function T(x) = e^x ln(1+e^-x)
% For x>30 a dramatic deviation from the
% Theoretical value is found; for x>36.7,
% the value drops spuriously to 0
x = linspace(33.5,37,333331);
y = exp(x).*log(1+exp(-x));
z = log(1+exp(-x))*10^(15)/3.5*2;
plot(x,y,x,z)
title('plot of f1=log(1+e^-1) with f2 = e^x*f1 ')