matlab实验四答案详解.docVIP

实验四 循环结构程序设计 实验内容 程序代码 （1）%循环结构 y=0; n=input(n=); for i=1:n y=y+1/i/i; end pi=sqrt(6*y) （2）%向量方法 n=input(n=); i=1./(1:n).^2; s=sum(i); pi=sqrt(6*s) 运行结果 （1） f n=100 pi = 3.1321 （2） f n=10000 pi = 3.1415 程序代码 y=0;n=0; while y3 n=n+1; y=y+1/(2*n-1); end y n if y3 n=n-1; end n 运行结果 f y = 3.0033 n = 57 n = 56 程序代码 a=input(请输入正数a=); b=input(请输入正数b=); x=1.0;n=0; while abs(x-a/(b+x))0.00001 if n500 x=a/(b+x); n=n+1; end end x r(1)=(-b+sqrt(b^2+4*a))/2 r(2)=(-b-sqrt(b^2+4*a))/2 s=r-x 运行结果 f 请输入正数a=1 请输入正数b=1 x = 0.6180 r = 0.6180 r =