gedebahe1.m

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

function [y,num]=gedebahe1(n)
if(rem(n,2)) disp('please input a even integral.');
else 
    m=0;
    for k=1:n/2
        if(isprime(k)&&isprime(n-k)) m=m+1;
            y(m,1)=k;
            y(m,2)=n-k;
        end
    end
end    
num=m