2014自写版Matlab-作业文档.docVIP

1.利用循环语句，计算数列，，，…的极限，要求误差小于10-8。

a=5;

a1=sqrt(5);

while((a-a1)=10^(-8))

a=a1;

a1=sqrt(a1);

n=n+1;

end

vpa(a1,15)

ans=

1.00000000599562

n

n=

画图：

n=1:1:29;

x=n;

y=5.^(0.5.^n);

plot(x,y)

2.利用级数可计算出无理数π的近似值。由于交错数级的部分和数列SN在级数的和上下摆动，截断误差将小于第一个被舍弃的项的绝对值。试分析，若要使截断误差小于10-4或10-8，应取多少项求和？计算π的近似值。

误差为10的-8次方：

n=1;

anspi=0;

while((1/(1+2*(n-1)))=(10^(-8)))

anspi=anspi-4*((-1)^n)*(1/(2*n-1));

n=n+1;

end

vpa(anspi,8)

ans=

3.1415926

n

n=舍去，可使级数误差小于10的-8次方此时π有效值为3.1415926

误差为10的-4次方时：

n=1;

anspi=0;