[求助]画不出图来！

[me4047]

请现看看我写的程序

function f=ps(M,h,k,sn)
for h=0:13
for sn=1:10
y=0;
for k=0:M-h
pe=erfc(sqrt(sn))/2;
d=nchoosek (M,k).*((1-pe).^(M-k)).*pe.^k;
y=y+d;
end
pm=1-y;


z=0;
for k=0:h
c=nchoosek (M,k);
z=y+c;
end
pf=z.*(1/2).^M;

ps=pm+pf
end
end

我费劲心血写出来的，不要笑话我啊:sweat:

我现在虽然能写出来ps的答案，但是画不出图啊！
请问怎么样才可以画出一个ps和h的图像呢 ？

感谢先