Can anyone tell me why I'm getting an error here? "xold=xnew;'' is giving me an error in this statement and I'm not sure why.

Question

Brendan 2023-2-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1916680-can-anyone-tell-me-why-i-m-getting-an-error-here-xold-xnew-is-giving-me-an-error-in-this-statem

回答： KSSV 2023-2-22

clear all 
close all
clc
% Problem 6.7
while (1)
fx=@(x) x.^3-6*x.^2+11*x-6.1;
% Using the Modified Secant Method 
xold=3.5;
d=0.01;
% Iteration 1
xold=xnew;
xnew=xold-d*xold*fx(xold)/(fx(xold+d*xold)-fx(xold));
ea=abs((xnew-xold)/xnew)*100;
% Iteration 2
xold=xnew;
xnew=xold-d*xold*fx(xold)/(fx(xold+d*xold)-fx(xold));
ea=abs((xnew-xold)/xnew)*100;
% Iteration 3
xold=xnew;
xnew=xold-d*xold*fx(xold)/(fx(xold+d*xold)-fx(xold));
ea=abs((xnew-xold)/xnew)*100;
fprintf('The root using 3 iterations of Modified Secant Method is %8.4f with an approximate relative error of %6.3e\n',xnew,ea)
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Brendan 2023-2-22

I would also like to keep an argument attached to this code.

请先登录，再进行评论。

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

Answer 1

KSSV 2023-2-22

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1916680-can-anyone-tell-me-why-i-m-getting-an-error-here-xold-xnew-is-giving-me-an-error-in-this-statem#answer_1177055

在 MATLAB Online 中打开

These lines:

xold=3.5;
d=0.01;
% Iteration 1
xold=xnew;
xnew=xold-d*xold*fx(xold)/(fx(xold+d*xold)-fx(xold));

will obviously thrwo an error becuase xnew is a new variable and it is not defined. You need to either comment this line or you may need to move to the next line. After calculating xnew then save xnew to xold.