(答案)2010专：程序设计一.docVIP

文件的命名： 如：2009040108 温玲玲(3) 程序设计题目一 一、试编程实现以下问题: 1.利用 公式求的近似值,直到误差小于。 k=1;s=0; while 4/(2*k-1)10^(-6) s=s+(-1)^(k+1)/(2*k-1); k=k+1; end pai=4*s 结果： pai = 3.14159215358972 2. 应用近似求之值，要求误差小于. 程序：k=1;s=1; while 1/prod(1:k)10^(-6) s=s+1/prod(1:k); k=k+1; end e=s 结果：e = 2.71828152557319 3. 求[100,200]之间的素数。 I.for循环程序： for m=101:2:200 k=fix(sqrt(m)); for i=2:k+1 if rem(m,i)==0 break; end end if i=k+1%除遍了所有（2到K+1） disp(int2str(m)) end end 或者：（以数组的形式显示） u=0; for m=101:2:200 k=fix(sqrt(m)); for i=2:k+1 if rem(m,i)==0 break; end end if i=k+1 u=u+1; a(u)=m; end end b=a II.while循环程序： m=101; while m200 k=fix(sqrt(m)); i=2; while i=k+1 if rem(m,i)==0 break; end i=i+1; end if i=k+1 disp(int2str(m)) end m=m+2; end 或： i=1; x=[]; for j=100:200 if isprime(j)==1%此命令可直接判断j是否为素数 x(i)=j; end i=i+1; end 4.若一个数等于它各个真因子之和，则称该数为完数，如6=1+2+3，所以6就是完数，求[1,500]之间的全部完数。 程序：for m=1:500 s=0; for k=1:m/2 if rem(m,k)==0 s=s+k; end end if m==s disp(m) end end 5.列出所有的水仙花数, 水仙花数是一个三位数,其各位数字立方和等于该数本身. 例如: 程序： 程序1： for a=1:9 for b=0:9 for c=0:9 abc=100*a+10*b+c; if abc==a.^3+b.^3+c.^3 disp(abc) end end end end 程序2： k=1; for a=1:9 for b=0:9 for c=0:9 abc=100*a+10*b+c; if abc==a.^3+b.^3+c.^3 shuixianhua(k)=abc; k=k+1; end end end end shuixianhua 结果： shuixianhua = 153 370 371 407 程序3： for m=100:999 m1=fix(m/100);%求m的百位数字 m2=rem(fix(m/10),10);%求m的10位数字 m3=rem(m,10); if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 d