matlab求10000以内的素数clear alldisp('10000以内全体素数:')j=2;count=0;for i=2:10000k=sqrt(i);while jif rem(i,j)==0breakelsej=j+1;endif j>kfprintf('%3d',i)count=count+1;endif rem(count,13)==0fprintf('\n')endendend请问哪里错了.怎么

来源:学生作业帮助网 编辑:六六作业网 时间:2024/12/04 01:38:03
matlab求10000以内的素数clearalldisp(''10000以内全体素数:'')j=2;count=0;fori=2:10000k=sqrt(i);whilejifrem(i,j)==0br

matlab求10000以内的素数clear alldisp('10000以内全体素数:')j=2;count=0;for i=2:10000k=sqrt(i);while jif rem(i,j)==0breakelsej=j+1;endif j>kfprintf('%3d',i)count=count+1;endif rem(count,13)==0fprintf('\n')endendend请问哪里错了.怎么
matlab求10000以内的素数
clear all
disp('10000以内全体素数:')
j=2;
count=0;
for i=2:10000
k=sqrt(i);
while j
if rem(i,j)==0
break
else
j=j+1;
end
if j>k
fprintf('%3d',i)
count=count+1;
end
if rem(count,13)==0
fprintf('\n')
end
end
end
请问哪里错了.
怎么改?
给30分.

matlab求10000以内的素数clear alldisp('10000以内全体素数:')j=2;count=0;for i=2:10000k=sqrt(i);while jif rem(i,j)==0breakelsej=j+1;endif j>kfprintf('%3d',i)count=count+1;endif rem(count,13)==0fprintf('\n')endendend请问哪里错了.怎么
改完的如下:
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