Introduction - If you have any usage issues, please Google them yourself
le=zeros(401,2);
i=0;
for u=0:0.01:4
i=i+1;
x=0.5;
y=0;
for j=1:100
x=u*x*(1-x);
end
for k=1:200
x=u*x*(1-x);
y=y+log(abs(1-2*x));
end
le(i,1)=u;
le(i,2)=log(u)+y/200;
end
plot(le(:,1),le(:,2));
hold on;
n=0:0.001:4;
plot(n,0,k);