Alpha=0:0.01:1;
Beta=0:0.01:1;
[alpha,beta]=meshgrid(Alpha,Beta);
PL=165500./((110+beta.*(1-alpha)).*(7+beta.*(1-alpha)));
surf(alpha,PL,beta)
hold on
PU=165500./((110-beta.*(1-alpha)).*(7-beta.*(1-alpha)));
surf(alpha,PU,beta)
Alpha=0:0.01:1;
Beta=0:0.01:1;
[alpha,beta]=meshgrid(Alpha,Beta);
PL=165500./((110+beta.*(1-alpha)).*(7+beta.*(1-alpha)));
surf(alpha,PL,beta)
hold on
PU=165500./((110-beta.*(1-alpha)).*(7-beta.*(1-alpha)));
surf(alpha,PU,beta)