0 个投票

Hi ,
I am trying to make newtons Method work.
I want it to work in a specific way , somthing like this
[s,fc] = mynewton (@somethingfunction,1,10^-4);
function [fv , fpv ] = somethingfunction(x)
% MYFUN Function to test Newton
%
fv = x^2 - 10;
fpv = 2*x;
end
function [x ]= mynewton (somethingfunction,x0 , tol )
fncon=1;
[f(1), g(1)] = myfun(x0);
s(1)=x0 ;%nitial point
for i=1:1000 %it should be stopped when tolerance is reached
[f(i), g(i)] = myfun(s(i));
s(i+1) = s(i) - f(i)/g(i);
if( abs(s(i+1)-s(i)/s(i))>tol) % tolerance
disp(double(s(i+1)));
break
end
end
if fncon
fprintf(' %16.16f \n',s(end))
end
end

请先登录，再回答此问题。

 采纳的回答

Ameer Hamza
Ameer Hamza 2020-6-7

0 个投票

Condition should be
if(abs(s(i+1)-s(i)/s(i))<tol)
% ^ less than operator

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Get Started with MATLAB 的更多信息

产品

版本

R2020a

标签

Community Treasure Hunt

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

Start Hunting!

Translated by