I need to use Newton's method to solve f(x). Can any one tell me what's wrong of my code and how to improve it.

1 view (last 30 days)
Chun Kit NG
Chun Kit NG on 9 Feb 2022
Answered: David Hill on 9 Feb 2022

Answers (2)

Yongjian Feng
Yongjian Feng on 9 Feb 2022
what i s your dx?

David Hill
David Hill on 9 Feb 2022
a=1:4;
tiny=1e-10;
x=ones(size(a));%this will not work well for the initial guess
for k=1:length(a)
f=@(x)tan(x)+a(k)*sin(x)./x;
fp=@(x)sec(x)^2-a(k)*sin(x)/(x^2)+a(k)*cos(x)/x;
rat=1;
while rat>tiny
dx=f(x(k))/fp(x(k));
x(k)=x(k)-dx;
rat=abs(dx/x(k));
end
end
If you look at your plot, you will realize that there are lots of crossings. Your initial guess for x must be close to the crossing you want.

Community Treasure Hunt

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

Start Hunting!

Translated by