用mathematica求 farey 数列程序
1个回答

Mathematica 程序:

farey[x_] :=

Module[{m,n,jieguo},

jieguo = Union[Sort[Flatten[Table[m/n,{n,1,x},{m,1,n}]]]];

jieguo

]

farey[12]

运行结果:

{1/12,1/11,1/10,1/9,1/8,1/7,1/6,2/11,1/5,2/9,1/4,3/11,2/7,3/10,1/3,4/11,3/8,2/5,5/12,3/7,4/9,5/11,1/2,6/11,5/9,4/7,7/12,3/5,5/8,7/11,2/3,7/10,5/7,8/11,3/4,7/9,4/5,9/11,5/6,6/7,7/8,8/9,9/10,10/11,11/12,1}