isprime.m

来自「关于质数的一些程序」· M 代码 · 共 13 行

function y=isprime(x)
k=length(x);
for m=1:k
    if(x(m)<1) y(m)=inf;disp('only for intergals which are bigger than zero');
    elseif(x(m)==1) y(m)=0;
    else
        y(m)=1;
        for n=2:x(m)-1
            if(rem(x(m),n)==0) y(m)=0;
            end
        end
    end
end