function res=Efun(x,n)
warning off
for k=1:length(x)
 f=@(mi) mi.^(n-2).*exp(-x(k)./mi);
% f=@(mi) mi.^(2-n).*exp(-x(k)./mi);
 res(k)=quad(f,0,1);
end 
warning on

