www.gusucode.com > smith圆图matlab源码程序 > smith/N1.m
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')