n1.m

smith圆图

x=factorial(16);

n=2:9;
a=1;
for t=1:n
a=a*t;
end 
y=a;

b=1;
for p=1:(16-n)
    b=b*p;
end
z=b;

v=x/(y*z)
plot(n,v,'ok')