clc
x=0:1/100:1;
y=0.02:1000/100:1000;
[x,y]=meshgrid(x,y);
t=-0.8./x+log(1-x)+0.8.*(1./(1-x)-1).*(1-1./x);
z=1./(1+exp(0.8-t)./y);
plot3(x,y,z)
%x,y,z分别代表α，β，E(α）