[proportion]

先做一个判断合质数的函数：
function k=zhs(x)
k=1;
if x<=1
disp('error!')
k=-1;
return
end
if x>=4
for i=2:sqrt(x)
if mod(x,i)==0
k=0;
break
end
end
end