matlab输出10000以内的素数,matlab求10000以内的素数

改完的如下：

clear all

disp('10000以内全体素数:')

count=0;

for i=2:10000

j=2; % 每次外循环,j都要从头再来,所以要拿到外循环里.

k=sqrt(i);

while j

if rem(i,j)==0

break

else

j=j+1;

end

if j>k

fprintf('%5d',i) % 给你改成%5d了,要不都挨一块分不出个儿来了.

count=count+1;

if rem(count,13)==0 % 这个if块,要放到判断j>k里,要不会多出好多空行

fprintf('\n')

end

break % 这个最重要,j>k时一定要退出内循环,要不然j会继续运行到等于i才会退出的

end

end

end

不过话说回来,matlab本身就带求n以内素数的函数啊,用primes

比如

>> primes(100)

ans =

Columns 1 through 17

2 3 5 7 11 13 17 19 23 29 31 37 41 43 47 53 59

Columns 18 through 25

61 67 71 73 79 83 89 97

如果不可以用primes的话,自己编也不用那么复杂啊

给你重新编了个简单的,不太好理解,不过程序简单.

n=10000;

p = 1:2:n; % 找出奇数,因为除了2的所有素数都是奇数

q = length(p);

p(1) = 2; % 第一个素数为2

for k=3:2:sqrt(n)

if p((k+1)/2) % 如果p((k+1)/2)也就是p中原本应该等于k的那个数不为0的话

p((k+1)/2+k:k:q)=0;

% 这句的作用就是把p中所有能被k整除的数设为0,但不包括k本身.

% (k+1)/2是奇数k在p中的位置,(k+1)/2+k代表p中除了k本身之外的第二个能被k整除的数的位置,依此类推

end

end

p=p(p>0); % 把等于0的,也就是原来位置上的数是合数的去掉

fprintf('%d以内的素数有：\n',n);

q=length(p);

for k=1:13:q

fprintf('%5d',p(k:min(q,k+12))); % 每行输出13个数

fprintf('\n');

end