[two_star]

这是判断10000以内成立的代码：
clear;
E = primes(10000);
num = 0;
for x = 2:2:10000
for i = 1:length(E)
if ismember( x - E(i),E) == 1
num = num + 1;
end
end
end