I can't get this IF statement to work

1 vue (au cours des 30 derniers jours)
Sarfraz Khan
Sarfraz Khan le 3 Mar 2020
Commenté : Rena Berman le 14 Mai 2020
Hi, I have been trying to programme newtons law of itterations, which is an itterative method which converges to the unknown point we are trying to figure out. I have got my code to do the hard part. I would however like my code to break as soon as the same figure begins to repeat itself. Please see below:
If I hadnt included the 1:100 statement the code would run on forever. Please help!!
f = @(x) exp(x)-1.5+atan(x);
g = @(x) exp(x)+1/((x^2)+1);
x0 = 50
for i = 1:100
x1 = x0 - (f(x0)/g(x0)) %Newtons Formula
if x0 == x1
else
x0=x1
end
end
  2 commentaires
Rik
Rik le 3 Mar 2020
Why did you remove your question? That is extremely rude. Now nobody can benefit from your question if they have a similar issue.
Rena Berman
Rena Berman le 14 Mai 2020
(Answers Dev) Restored edit

Connectez-vous pour commenter.

Réponse acceptée

David Hill
David Hill le 3 Mar 2020
f = @(x) exp(x)-1.5+atan(x);
g = @(x) exp(x)+1/((x^2)+1);
x0 = 50;
tol=1e-6;
for i = 1:100
x1 = x0 - (f(x0)/g(x0)); %Newtons Formula
if abs(x0-x1)<tol
break;
else
x0=x1;
end
end
You could also do a while loop:
  2 commentaires
Sarfraz Khan
Sarfraz Khan le 3 Mar 2020
Thanks so much for your previous answer.
How would I use a while loop here? I have attempted but getting similar issues :(
David Hill
David Hill le 3 Mar 2020
root=x-f(x)/fp(x);
while abs(root-x)>tol
x=root;
root=x-f(x)/fp(x);
end

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Programming dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by