[zhd87]

M2=10:1:50;
M1=10:1:50;
。。。。。。
。。。。。。
。。。。。。
for i=0:wl
pe_theo= pe_theo+ 1/2* (-1)^(wl-i).* nchoosek(wl, i)* (q0_w1+ q1_w1*i/wl).^(Ml-1).* (q0_w2+ q1_w2*i/wl).^(Mh);

[M2,M1]=meshgrid(M2,M1);
mesh(M2,M1,pr_theo);
set(gca,'ZScale','log');
grid on;
hold on;




错误提示：
??? Error using ==> mesh at 80
Z must be a matrix, not a scalar or vector.

Error in ==> figure1 at 65
mesh(M2,M1,pr_theo);