clc;clear
x(1)=1;
y(1)=2;
tol=1e-6;
n=1;
while n<=10000
x(n+1)=sqrt(x(n)+y(n));
y(n+1)=0.5*(x(n)+y(n));
dx=x(n+1)-x(n);
dy=y(n+1)-y(n);
if abs(dx)<=tol && abs(dy)<=tol
disp('xn,yn数列均收敛,其收敛值分别为:')
x(n)
y(n)
x'
y'
break
elseif n==10000 && abs(dx)<=tol && abs(dy)>tol
disp('xn数列收敛,yn数列不收敛,其收敛值为:')
x(n)
elseif n==10000 && abs(dy)<=tol && abs(dx)>tol
disp('xn数列不收敛,yn数列收敛,其收敛值为:')
y(n)
elseif n==10000 && abs(dx)>tol && abs(dx)>tol
disp('xn,yn数列均不收敛')
end
n=n+1;
end
结果均收敛于2