MATLABのセカント根発見関数の問題

Question

私はこの関数をしばらく使用していませんでした。私が持っている問題を回避する方法を見つけることができません。

私のコードは次のとおりです。私は、コードを実行すると、私は手で行うよう

function [x,i] = secant(f, x0, x1, tol, maxiters) 
%SECANT Secant method 
% [X,I] = SECANT(F, X0, X1, TOL, MAXITERS) performs the secant 
% method with F(x), starting at x_0 = X0 and x_1 = X1, and continuing 
% until either |X_i+1 - X_i| <= TOL, or MAXITERS iterations have 
% been taken. The number of iterations, I, is also returned. 
% An error is raised if the first input is not a function handle. 
% A warning is raised if the maximum number of iterations is reached 
% without achieving the tolerance. 

if ~isa(f, 'function_handle') 
error('Your first input was not a function handle') 
end 

i = 0; % initialise iterate counter 
x = x1; 
x_old = x0; 
while abs(x - x_old) > tol && i < maxiters 
x_old = x0; 
x = x - f(x)*(x - x_old)/(f(x) - f(x_old)); % compute the new x 
i = i + 1; % increase our counter by one 

end 
if abs(x - x_old) > tol 
warning('Maximum number of iterations reached without achieving tolerance.') 
end 

が、私は同じ答えを得ることはありません。現在、私は、私は私の問題は下while abs(x - x_old) > tol && i < maxiters .Iあるラインx_old = x0;から来ていると信じて、

f = @(x) 0.9cos(x) - sqrt(x) 

x0 = 0 

x1 = 1 

tol = 1e-08 

maxiters = 4 

を使用しています、私は関数を呼び出すたびに、それが戻っx0からxを設定だと思います。

どうすればこの問題を回避できますか？

Answer 1

問題は前述したとおりです。 x_oldは実際にはすべての反復iに対してx_i-2に設定する必要があるときは常に0に設定されます。これは、変更する必要がある関連コードです

i = 0; % initialise iterate counter 
%Store the initial values in a different variable 
x_0 = x0; 
x_1 = x1; 
while abs(x - x_old) > tol && i < maxiters 
    x = x_1 - f(x_1) * (x_1 - x_0)/(f(x_1) - f(x_0)); % compute the new x 
    % Update the previous values 
    x_0 = x_1; 
    x_1 = x; 
    i = i + 1; % increase our counter by one 
end