2334.sagews

CWimax =  73.200;
CWifi =  11;
theta1,theta2 = var('theta1 theta2');
Utility=((CWimax*theta1)/(1+theta2))+(1-theta1)/((2-theta2)/CWifi)+((CWimax*theta2)/(1+theta1))+(1-theta2)/((2-theta1)/CWifi);

plot3d(Utility,(theta1,0,1),(theta2,0,1), mesh=True)

hess=Utility.hessian()

valprop=hess.eigenvalues()

plot3d(valprop[0],(theta1,0,1),(theta2,0,1), mesh=True)

plot3d(valprop[1],(theta1,0,1),(theta2,0,1), mesh=True)