l=0:0.01:1;
for n=[0.1,1,10,100]
Up=l./(1+l.*(1-l)*n)
plot(l,Up)
hold on
end
xlabel('l');
ylabel('Up');