Trying to use Newton's method until solution converges
1 view (last 30 days)
Show older comments
So this is probably a dumb question, but I'm terrible at Matlab and am really struggling and can't seem to find any examples for my problem. I understand how Newton's method works, but not sure how to apply it in my case. My problem is I'm trying to solve this equation
for delta_sigma . All the other variables are known already.
This is the code I have right now:
% Variables
Eps = 0.011;
E = 73100;
H = 662;
n = 0.07;
maxIter = 100; % # of iterations
i = 1;
x(i) = 500; % Initial guess
while i <= maxIter
Eps = x/E+2.*(x./(2*H)).^(1/n);
dEps = 2./(n.*x).*(x./(2*H)).^(1/n);
y = x(i)-(Eps/dEps);
i = i+1;
x(i) = y;
end
which I think would be correct if I was trying to find the delta_epsilon value. Could someone help me out? Thanks
2 Comments
Accepted Answer
Torsten
on 25 Jan 2023
Edited: Torsten
on 25 Jan 2023
% Variables
Eps = 0.011;
E = 73100;
H = 662;
n = 0.07;
maxIter = 100; % # of iterations
Tol = 1e-8;
i = 1;
error = 1.0;
X(i) = 500; % Initial guess
x = X(i);
while i <= maxIter && error > Tol
f = x/E + 2*(x/(2*H))^(1/n) - Eps;
df = 1/E + 2 * 1/n * (x/(2*H))^(1/n-1) * 1/(2*H);
x = X(i) - f/df;
error = abs(x-X(i));
i = i+1;
X(i) = x;
end
i
X
x/E + 2*(x/(2*H))^(1/n)- Eps
4 Comments
Matt J
on 26 Jan 2023
Wow thanks so much!
@Christina Kersten If so, you should accept-click Torsten's asnwer
More Answers (0)
See Also
Categories
Find more on Interpolation in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!