%M函数
function z=testfun(a,b,c)
z=(prod(1:b)/(prod(1:a)*prod(1:(b-a))))*(c^a)*((1-c)^(b-a));
end

程序代码

clc;clear;
syms x y
N=40;
M=10;
c=1;
L=0.3;
A=N*(y/c+(1-L/M)*1.5*y/c)*testfun(M-1,N-1,x);
g=1;
for i=0:M-1
    g=g+((i/M)-1)*testfun(i,N,x);
end
fun=(A-sum(g));
ezplot(fun,[0,1,0,10]);
hold on

以上是绘出x为横坐标，y为纵坐标的图像，请问各位怎样绘出x为横坐标，N为纵坐标的图像，其中N为自然数，y可取值为2时。
多谢！！