c1=0:1:150;
c2=0:1:150;
[c1,c2]= meshgrid(c1,c2);
lamda=(-(c1+100)+((c1+100).^2+1200000*c2).^0.5)/600000;
beta1=300000-100-c1-300000*lamda;
beta2=300000*lamda-c2;
mesh(c1,c2,beta1), hold on;
mesh(c1,c2,beta2), hold off;
xlabel('\bfc1');
ylabel('\bfc2');
zlabel('\bf利润')