函數f=erti(x)
h =排序(x)
n = 1;
對於i=1:長度(h)
v = PDS szy(h(I));
如果v==1
f(n)= h(I);%
n = n+1;
結束
h
f
函數f=pdsszy(x)
f=1
對於i=2:x-1
n=x/i
if fix(n)==n
f=0
最後兩個上照