I am new to matab. I tried to solve a non linear equation with newton raphson and modified newton raphson method but got different solution. could you please help me find the error? please note I tried to run the code in two separate files.
    5 vues (au cours des 30 derniers jours)
  
       Afficher commentaires plus anciens
    
%%Newton raphson method
clc
clear all
close all
%%Initial condition
x0 =0.92;
maxiteration = 50;
tolX = 1e-04;
f = @(x) x^2-2*x+1; %function to be solved
%%Computation
x = x0;
xold = x0;
for i = 1:maxiteration;
    df = 2*x-2;
    x = x-f(x)/df;
    err(i) = norm(x-xold);
    xold = x;
    if (err(i)<tolX)
        break;
    end
end
msg = ['The solution converged in ', num2str(i),'th iterations.'];
disp(msg)
msg1 = ['converged numerical solution(critical values):'];
disp(msg1)
fprintf('x = %.4f\n',x);  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%Modified newton Raphson method
clc
clear all
close all
%%Initial condition
x0 =0.92;
maxiteration = 50;
tolX = 1e-04;
f = @(x) x^2-2*x+1;
%%Computation
x = x0;
xold = x0;
for i = 1:maxiteration;
    df = 2*x0-2;
    x = x-f(x)/df;
    err(i) = norm(x-xold);
    xold = x;
    if (err(i)<tolX)
        break;
    end
end
msg = ['The solution converged in ', num2str(i),'th iterations.'];
disp(msg)
msg1 = ['converged numerical solution:'];
disp(msg1)
fprintf('x = %.4f\n',x)
0 commentaires
Réponses (1)
  darova
      
      
 le 23 Déc 2019
        Plot each iteration to see what it going on
fplot(f,[-5 5])
hold on
for i = 1:maxiteration;
    % computations ...
    yy = df*[-1 1]+f(x);
    h(1) = plot([-1 1]+x,yy);
    h(2) = plot(x,f(x),'.r');
    pause(1)
    delete(h)
end
hold off
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!

