matlab求素数

%任意输入一个正整数m，要找出2~m之间的全部素数，首先在2~m中划去2的倍数
% (不包括2)，然后划去 3的倍数(不包括3)，由于4已被划去，
%再找5的倍数 (不包括5)，…，直到再划去不超过√m的数的倍数，
%剩下的数都是素数。

function retval = sushu (input1, input2)
m = input('m=');
for j=1:m
  a(j)=j;       
  end
for i=2:m
  for k=2:i
  for n=2:i
    if k*n==i
      a(i)=0;
    end
  end
end
end
disp(find (a!=0));
endfunction

 

因为前期学过c语言，所以思想是用c语言来写的