实验一(3)

2）程序 [u,v]=meshgrid([0:0.03*pi:2*pi]) X=(1+cos(u)).*cos(v); Y=(1+cos(u)).*sin(v); Z=sin(u); mesh(X,Y,Z) 图像： 4）程序： [X,Y]=meshgrid([-10:0.5:10]); Z=Y.^2; mesh(X,Y,Z); 图像： 5）程序： [Z,X]=meshgrid([0:0.2:5]); Y=-Z.^2; mesh(X,Y,Z); 图像： 6．建立一个命令M-文件：求所有的“水仙花数”，所谓“水仙花数”是指一个三位数，其各位数字的立333方和等于该数本身。例如，153是一个水仙花数，因为153=1+5+3。 解：思路如下：三位数m=个位z^3+十位y^3+百位x^3, 通过对xyz三个循环穷举，可以找出水仙花数。 程序如下： for x=1:9 for y=0:9 for z=0:9 m=x.*100+y.*10+z; if m==x.^3+y.^3+z.^3; m end end end end 结果如下： 所以水仙花数为153,370,371,407. 7．编写函数M-文件sq.m：用迭代法求x?a的值。求平方根的迭代公式为 xn?1?1a(xn?) 2xn5迭代的终止条件为前后两次求出的x的差的绝对值小于10?。 程序如下： function f=sq(a) x0=1; x1=0.5*(x0+a/x0); n=0; while (abs(x1-x0)>=10^(-5)) x0=x1; x1=0.5*(x0+a/x0); n=n+1; end; x0 x1 n