function [mu delta] = findmu2(n)





mu = linspace(0.75,1,5000);


delta = zeros(1,length(mu));





for j = 1:length(mu)


    [i x] = iterate(mu(j), 0.5, 2^n+1);


    delta(j) = x(2^n+1)-0.5;


end