重庆大学 数学实验 温罗生(2)

图形 5.代码 [x,y]=meshgrid(-5:0.1:5); z=sin(pi*sqrt(x.^2+4*y.^2)); mesh(x,y,z) 图形 6.代码 H=[]; n=1; for a=1:9 for b=0:9 for c=0:9 h=100*a+10*b+c; j=a^3+b^3+c^3; if h==j H(n)=j; n=n+1; end end end end H 结果： H = 153 370 371 407 7.代码: (1) function f=judge(x) f=1; if x<4 else for i=2:x-1 k=mod(x,i); if k==0 f=0; end end end end (2) n=0; H=[]; for i=2:20000 f=judge(i); if f==1 n=n+1; H(n)=i; end end H L=size(H) hist(H,10) 等距区间内质数的个数的直方图（以1000为间距） 8.代码： function [f1,f2,f3]=PI(x) h=[]; H=[3,1,3]; for i=10^(x-1):10^x-1 j=fix(i/3); k=fix(i/3.3); for l=k:j h=[i,l,i/l]; if abs(h(3)-pi)