how can i get the graph of this iterative sequence
t=@(n)(1/(n+2)^2);r=1/5;
g=@(x)((x-1)/(exp(x)));
tolerance=0.000001;
k=100;
x=zeros(1,k);
x(1)=0.3977;
for n=1:k
x(n+1)=(1-(r)*g(x(n)))*(1-t(n))*x(n);
if abs(x(n+1)-x(n))<tolerance
fprintf('The sequence converges after %d iterations for tolerance = %f', n, tolerance)
break;
end
end
x=x(1:n+1); %discard the zeros
format long
disp(x)

 採用された回答

Torsten
Torsten 2023 年 5 月 20 日
編集済み: Torsten 2023 年 5 月 20 日

0 投票

Add the command
plot(1:n+1,x)
at the end of your code.

2 件のコメント

Basha
Basha 2023 年 5 月 21 日
is it possible to dispaly more than one graph by changing x(1)?
Torsten
Torsten 2023 年 5 月 21 日
t=@(n)(1/(n+2)^2);r=1/5;
g=@(x)((x-1)/(exp(x)));
tolerance=0.000001;
k=100;
N=20;
hold on
for NN=1:N
x(1)=rand;
for n=1:k
x(n+1)=(1-(r)*g(x(n)))*(1-t(n))*x(n);
if abs(x(n+1)-x(n))<tolerance
fprintf('The sequence converges after %d iterations for tolerance = %f', n, tolerance)
break;
end
end
plot(1:n+1,x)
end
hold off
grid on

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

ヘルプ センター および File ExchangeApp Building についてさらに検索

製品

リリース

R2023a

タグ

質問済み:

2023 年 5 月 20 日

コメント済み:

2023 年 5 月 21 日

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by