p1=10; c2=12; c3=10; M=10000; b1=40; b2=30; k=2000; w=0.2; p0=20;a=0.3;
m=0:0.001:1;
p20=(k*w*(M+b1*(c3+p0))-b2^2*(c2-p0))/2/b1/k/w;
p21=-(b2^2*p1 - b2^2*c2 + M*k*w - b2^2*c3*m - b2^2*m*p1 + b1*c3*k*w + b1*k*p1*w)./(m*b2^2 - 2*b1*k*w);
ph1=plot(m,p20);hold on
ph2=plot(m,p21);xlabel('m'),ylabel('p2');hold on
lh=legend([ph1 ph2],{'无契约M=10000';'有契约M=10000'});


为什么总显示错误