南邮matlab数学实验~第二次练习题.doc

第二次练习题 1、 设，数列是否收敛？若收敛，其值为多少？精确到6位有效数字。 f=inline((x+7/x)/2); x1=3; for i=1:20 x1=f(x1); fprintf(%g %g\n,i,x1); end 1 2.66667 2 2.64583 3 2.64575 4 2.64575 5 2.64575 6 2.64575 7 2.64575 8 2.64575 9 2.64575 10 2.64575 11 2.64575 12 2.64575 13 2.64575 14 2.64575 15 2.64575 16 2.64575 17 2.64575 18 2.64575 19 2.64575 20 2.64575 2、设 是否收敛？若收敛，其值为多少？精确到17位有效数字。 注：学号为单号的取，学号为双号的取 s=0;xn=1;n=1; while n50 s=s+xn; fprintf(%g %g\n,n,s); n=n+1; xn=1/n^7; end 1 1 2 1.00781 3 1.00827 4 1.00833 5 1.00834 6 1.00835 7 1.00835 8 1.00835 9 1.00835 10 1.00835 11 1.00835 12 1.00835 13 1.00835 14 1.00835 15 1.00835 16 1.00835 17 1.00835 18 1.00835 19 1.00835 20 1.00835 21 1.00835 22 1.00835 23 1.00835 24 1.00835 25 1.00835 26 1.00835 27 1.00835 28 1.00835 29 1.00835 30 1.00835 31 1.00835 32 1.00835 33 1.00835 34 1.00835 35 1.00835 36 1.00835 37 1.00835 38 1.00835 39 1.00835 40 1.00835 41 1.00835 42 1.00835 43 1.00835 44 1.00835 45 1.00835 46 1.00835 47 1.00835 48 1.00835 49 1.00835 vpa(s,17) ans = 1.0083492773706009 书上习题：(实验四) f=inline((x-1)/(x+1)); x0=3.8; for i=1:20 x0=f(x0); fprintf(%g %g\n,i,x0); end 1 0.583333 2 -0.263158 3 -1.71429 4 3.8 5 0.583333 6 -0.263158 7 -1.71429 8 3.8 9 0.583333 10 -0.263158 11 -1.71429 12 3.8 13 0.583333 14 -0.263158 15 -1.71429 16 3.8 17 0.583333 18 -0.263158 19 -1.71429 20 3.8 2， y1=(x1-1)/(x1+3);y2=(-x2+15)/(x2+1); f1=inline((x1-1)/(x1+3));f2=inline((-x2+15)/(x2+1)); solve(y1-x1) ans = -1 -1 solve(y2-x2) ans = -5 3 f1=inline((x1-1)/(x1+3));f2=inline((-x2+15)/(x2+1)); x10=2.1;x20=4.3; while i201 x10=f1(x10); x20=f2(x20); fprintf(%g %g, %g\n,i,x10,x20); i=i+1; end 1 0.215686 2.01887 2 -0.243902 4.3 3 -0.451327 2.01887 4 -0.569444 4.3 5 -0.645714 2.01887 6 -0.699029 4.3 7 -0.738397 2.01887 8 -0.768657 4.3 9 -0.792642