Trying to use Newton's method until solution converges
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Christina Kersten
el 25 de En. de 2023
Comentada: Matt J
el 26 de En. de 2023
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 comentarios
Respuesta aceptada
Torsten
el 25 de En. de 2023
Editada: Torsten
el 25 de En. de 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 comentarios
Matt J
el 26 de En. de 2023
Wow thanks so much!
@Christina Kersten If so, you should accept-click Torsten's asnwer
Más respuestas (0)
Ver también
Categorías
Más información sobre Interpolation en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!