(2) z?sin(xy)0?x?3,0?y?3
>> x=0:0.1:3;y=0:0.1:3; [X Y]=meshgrid(x,y); Z=sin(X*Y); >> mesh(X,Y,Z)
?x?sint(3?cosu)（3）??y?cost(3?cosu)?z?sinu?0?t?2?
0?u?2?ezmesh('sin(t)*(3+cos(u))','cos(t)*(3+cos(u))','sin(u)',[0,2*pi,0,2*pi])
第二次练习题
7?x?(x?)/2n?n?1xn设?，数列{xn}是否收敛？若收敛，其值为多
?x?3?11、
少？精确到6位有效数字。 1.
>> f=inline('(x+917/x)/2');x0=3; >> for i=1:20
x0=f(x0);fprintf('%g,%g\\n',i,x0); end 1,154.333 2,80.1375 3,45.7902 4,32.9082 5,30.3868 6,30.2822 7,30.282
8,30.282 9,30.282 10,30.282 11,30.282 12,30.282 13,30.282 14,30.282 15,30.282 16,30.282 17,30.282 18,30.282 19,30.282 20,30.282
本次计算运行到第三次结果稳定，可得： 数列{xn}收敛，收敛到30.2820 2、
设 xn?1?12p?13p???1np, {xn}是否收敛？若收敛，其值为多少？精确到17位有效数字。 学号为单号，取p?7 >> s=0; for i=1:1:200
s=s+1/i^7;
fprintf('%g,%.17g\\n',i,s); end 1,1 2,1.0078125
3,1.0082697473708275 4,1.0083307825270775 5,1.0083435825270775 ??????????? 181,1.0083492773819187 182,1.0083492773819189 183,1.0083492773819192 184,1.0083492773819194 185,1.0083492773819196 186,1.0083492773819198 187,1.00834927738192 188,1.0083492773819203 189,1.0083492773819205 190,1.0083492773819207
191,1.0083492773819207 192,1.0083492773819207 193,1.0083492773819207 194,1.0083492773819207 195,1.0083492773819207 196,1.0083492773819207 197,1.0083492773819207 198,1.0083492773819207 199,1.0083492773819207 200,1.0083492773819207
运行至第190次后稳定，值为1.00834927738192070
书上习题：(实验四)
1，2，4，7（1），8，12（改为：对例2，取 a?4.5,25,55,120形有什么变化．），13，14 。 练习1 编程判断函数f(x)?x?1x?1的迭代序列是否收敛． >> f=inline('(x-1)/(x+1)'); x0=4; for i=1:20 x0=f(x0);
观察图