Question

有公式：x(n+1)=a/(b+xn).编写程序求迭代的结果，迭代的终止条件为｜x(n+1)-xn｜≤10^(-5),迭代初值x0=1.0

迭代次数不超过500次

Answer 1

x0=1;
a=input('请输入正数a'); b=input('请输入正数b');
if a<=0 | b<=0
disp('输入错误');
end
c=(-b+sqrt(b*b+4*a))/2 d=(-b-sqrt(b*b+4*a))/2
for i=1:1:500
x1=a/(b+x0);
if abs(x1-x0)<=0.00001
break;
x0=x1;
end
end