专业 姓名 学号 成绩
执行Martin（10，10000，-10000，20000），得：
40003000200010000-1000-2000-3000-4000-5000-4000-3000-2000-1000010002000300040005000
执行Martin（100，1000，-10，20000），得：
5004003002001000-100-200-300-400-500-600-400-2000200400600800
数学实验实验报告
21
专业 姓名 学号 成绩
执行Martin（-1000，17，4，20000），得：
0-200-400-600-800-1000-1200-1200-1000-800-600-400-2000200
2.6 能否找到分式函数
ax?b（其中a,b,c,d,e是整数）,使它产生的
cx2?dx?e迭代序列（迭代的初始值也是整数）收敛到3m(对于3m为整数的学号，请改为求310m)。如果迭代收敛,那么迭代的初值与收敛的速度有什么关系.写出你做此题的体会.
提示：教材54页练习4的一些分析。 若分式线性函数f(x)?的不动点，因此
ax?b的迭代收敛到指定的数2，则2为f(x)cx?d2?a2?b c2?d化简得：(2c?b)?(d?a)2?0。 若a,b,c,d为整数，易见b?2c,d?a。
取满足这种条件的不同的a,b,c,d以及迭代初值进行编。
数学实验实验报告
22
专业 姓名 学号 成绩
解：取m=10000；根据上述提示，取： a=e=1,b=10000,c=1,d=0. 程序如下（初值为1200）： f=inline('(x+926)/(x^2+1)'); x0=1; for i=1:100; x0=f(x0);
fprintf('%g,%g\\n',i,x0); end
运行结果如下： 1,0.00777777 2,9999.4
3,0.000200018 4,10000 5,0.0002 6,10000 7,0.0002 8,10000 9,0.0002 10,10000 11,0.0002 12,10000 13,0.0002 14,10000 15,0.0002 16,10000 17,0.0002 18,10000 19,0.0002 20,10000 21,0.0002 22,10000 23,0.0002 24,10000 25,0.0002
23
数学实验实验报告
专业 姓名 学号 成绩
26,10000 27,0.0002 28,10000 29,0.0002 30,10000 31,0.0002 32,10000 33,0.0002 34,10000 35,0.0002 36,10000 37,0.0002 38,10000 39,0.0002 40,10000 41,0.0002 42,10000 43,0.0002 44,10000 45,0.0002 46,10000 47,0.0002 48,10000 49,0.0002 50,10000 51,0.0002 52,10000 53,0.0002 54,10000 55,0.0002 56,10000 57,0.0002 58,10000 59,0.0002 60,10000 61,0.0002 62,10000 63,0.0002
24
数学实验实验报告
专业 姓名 学号 成绩
64,10000 65,0.0002 66,10000 67,0.0002 68,10000 69,0.0002 70,10000 71,0.0002 72,10000 73,0.0002 74,10000 75,0.0002 76,10000 77,0.0002 78,10000 79,0.0002 80,10000 81,0.0002 82,10000 83,0.0002 84,10000 85,0.0002 86,10000 87,0.0002 88,10000 89,0.0002 90,10000 91,0.0002 92,10000 93,0.0002 94,10000 95,0.0002 96,10000 97,0.0002 98,10000 99,0.0002 100,10000
25
数学实验实验报告