Newton Raphson on Mathlab

I was wondering how to solve this equation in matlab? the abc would be replaced by 742. Any help would be deeply appreciated.

Answer 1

[Torsten]

p=[0.5 1 -10 14.742];

r=roots(p);

Best wishes

Torsten.

[Sunny]

Yes I have the root now to solve the equation, the root is -6 how do I now plug this into the equation and solve it?

[Matt Tearle]

What do you mean? The roots in r already solve the equation. The real root near -6 is the first element of r:

>> fprintf('%10.6f\n',r(1))
 -6.083918

You can use the polyval function to evaluate a polynomial at a given value of x. For example:

x = linspace(-8,6);
y = polyval(p,x);
plot(x,y)
grid on

If you evaluate at x = r(1) you'll see that it solves the equation (to within ~10^-13):

polyval(p,r(1))

Answer 2

[Matt Tearle]

Torsten showed how to find the roots of a polynomial. More generally, to find the zeros of any function:

f = @(x) 0.5*x.^3 + x.^2 - 10*x + 14.742;  % define the function
x0 = -4;                                   % initial guess for the solution
xroot = fzero(f,x0)                        % solve

(If you're interested in the details, according to the doc, "the [fzero] algorithm, created by T. Dekker, uses a combination of bisection, secant, and inverse quadratic interpolation methods")

[Sunny]

So to complete this manually x0=-4?

[Matt Tearle]

I don't understand what you mean by "complete this manually". The assignment problem says to solve the equation. One way is to do N-R by hand. The other is "repeat using MATLAB". So I guess the question is what "repeat" entails. Repeat the solving of the equation? Or repeat solving the equation with N-R?

We've shown two ways you can solve the equation in MATLAB: roots (for solving polynomial equations) and fzero (for solving general nonlinear equations), but neither of these use N-R.

If you want to implement Newton-Raphson in MATLAB then that's a bigger issue. That requires knowing the basics of MATLAB programming. Given that this is a homework problem, you'll need to show what you've attempted and ask for some specific pointers.

Answer 3

[Sunny]

I've done the newton raphson manually on paper and I done it on matlab, Can someone help me compare the two on similarities or the differences?

Answer 4

[Torsten]

As Matt already mentioned, neither MATLAB's "fzero" nor MATLAB's "roots" uses N-R.

So it's difficult to compare the two methods.

You should ask your supervisor what he/she means by

"Repeat using MATLAB and compare the two methods.".

Best wishes

Torsten.

Answer 5

[Farhad Sedaghati]

The following link is the matlab code to perform Newton-Raphson's Method:

http://www.mathworks.com/matlabcentral/fileexchange/52362-newton-s-method

Answer 6

[sabik EL YATIM]

hello, Recently, a part of the Matlab code I found on the resolution system of nonlinear equations using the method of Newton-Raphson with the Jacobian matrix (I also left it in my comments). However, although he provides me with the basic code, I can not make it work, no matter how hard I try. I spent many hours trying to present function func, but to no avail, I often did not understand how I used it if you can help me

function [x,F, niter] = newtonsys(Ffun ,Jfun,x0 ,tol ,...
nmax, varargin )
% NEWTONSYS cherche un zéro d’un système non linéaire
% [ZERO ,F, NITER ]= NEWTONSYS(FFUN, JFUN,X0 ,TOL , NMAX)
% tente de trouver le vecteur ZERO, racine d’ un
% système non linéaire défini dans FFUN et dont
% la matrice jacobienne est définie dans la
% fonction JFUN. La racine est cherchée autour
% du vecteur X0.
% La variable F renvoie le résidu dans ZERO
% NITER renvoie le nombre d’ itérations nécessaires
% pour calculer ZERO. FFUN et JFUN sont des fonctions
% MATLAB définies dans des M- files.
niter = 0; err = tol + 1; x = x0;
while err >= tol & niter < nmax
J = feval(Jfun ,x , varargin {:});
F = feval(Ffun ,x , varargin {:});
delta = - J\F;
x = x + delta;
err = norm( delta );
niter = niter + 1;
end
F = norm( feval( Ffun,x, varargin {:}));
if (niter == nmax & err > tol)
fprintf ([’Pas de convergence dans le nombre ’ ,...
’ d’’ iterations imparti \n ’]);
fprintf ([’ La valeur retournée a un résidu ’ ,...
’ relatif de %e\n’],F);
else
fprintf ([’La méthode a convergé à l’’ itération’ ,...
’ %i avec un résidu %e\n’],niter,F);
end
return

Answer 7

[Meysam Mahooti]

https://www.mathworks.com/matlabcentral/fileexchange/73563-newton-raphson-method