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 (최근 30일)
이전 댓글 표시
답변 (2개)
David Hill
2022년 2월 9일
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.
댓글 수: 0
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!