运行结果如下: K>> %循环结构计算pi值 y=0; n=input('n='); for i=1:n y=y+1/i/i; end pi=sqrt(6*y) n=100 pi =
3.1321 n=1000 pi = 3.1406 n=10000 pi = 3.1415 %向量方法计算Pi值 n=input('n='); i=1./(1:n).^2; s=sum(i); pi=sqrt(6*s) n=100 pi = 3.1321 n=1000 pi = 3.1406 n=10000 pi = 3.1415
2. 根据y?1?111,求: ????352n?1
(1) y<3时的最大n值。
(2) 与(1)的n值对应的y值。 解:M—文件如下:
运行结果如下: K>> y=0;n=0; while y<3 n=n+1; y=y+1/(2*n-1); end y n if y>3 n=n-1; end n y = 3.0033
n = 57 n = 56
3. 考虑以下迭代公式:
xn?1?a b?xn其中a、b为正的学数。
(1) 编写程序求迭代的结果,迭代的终止条件为|xn+1-xn|≤10-5,迭代初值x0=1.0,迭代次数不超过500次。
?b?b2?4a(2) 如果迭代过程收敛于r,那么r的准确值是,当(a,b)的值取(1,1)、
2(8,3)、(10,0.1)时,分别对迭代结果和准确值进行比较。
解:
M文件如下:
运算结果如下; 请输入正数a=1 请输入正数b=1 x = 0.6180 r = 0.6180 -4.7016 r = 0.6180 -1.6180 s =