专业 姓名 学号 成绩
-(33 x (x - 1))/1010.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
当?=3.5时，上述程序稍加修改，得：
-(7 x (x - 1))/210.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
当?=3.56时，得：
数学实验实验报告
16
专业 姓名 学号 成绩
-(89 x (x - 1))/2510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
当?=3.568时，得：
-(446 x (x - 1))/12510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
数学实验实验报告
17
专业 姓名 学号 成绩
当?=3.6时，得：
-(18 x (x - 1))/510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
当?=3.84时，得：
-(96 x (x - 1))/2510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91
数学实验实验报告
18
专业 姓名 学号 成绩
2.5 对于Martin迭代，取参数a,b,c为其它的值会得到什么图形？参考下表（取自63页练习13）
a m -m -m m/1000 m/1000 m/100 -m/10 b m -m m/1000 m/1000 m m/10 17 c m m -m 0.5 -m -10 4
解：取m=627；迭代次数N=20000； 在M-文件里面输入代码：
function Martin(a,b,c,N)
f=@(x,y)(y-sign(x)*sqrt(abs(b*x-c))); g=@(x)(a-x); m=[0;0]; for n=1:N
m(:,n+1)=[f(m(1,n),m(2,n)),g(m(1,n))]; end
plot(m(1,:),m(2,:),'kx'); axis equal
在命令窗口中执行Martin（10000，10000，10000，20000），得：
数学实验实验报告
19
专业 姓名 学号 成绩
2500020000150001000050000-5000-2-1.5-1-0.500.511.52x 104
执行Martin（-10000，-10000，10000，20000），得：
50000-5000-10000-15000-20000-25000-2-1.5-1-0.500.511.52x 104
数学实验实验报告
20