用Matlab制作古诺模型企业利润与企业生产数量的三维图像

设价格函数中参数a=100，市场价格p=a-q1-q2,单位成本c=5,企业1生产的数量为q1，企业2生产的数量为q2。

企业1的利润π （q1，q2）=

在matlab中建立以q1和q2为x轴和y轴的平面，取值范围为[0,100]。代码为：[x,y]=meshgrid(1:0.1:100,1:0.1:100)。

再将企业1和企业2的利润设为z和v，并都用x和y表示出来代码为：

z=(-x.^2+(95-y).*x).*(x+y<=100)+(-5*x).*(x+y>100)

v=(-y.^2+(95-x).*y).*(x+y<=100)+(-5*y).*(x+y>100)

然后画出三维图像。

surf(x,y,z)

hold on

surf(x,y,v)

shading interp