clear all 

for n=[50 100 150];
n
h=1./n   
p1=log2((e_num(h)-e_num(h/2))./(e_num(h/2)-e_num(h/4))) 
end