Home MathLab Newton Raphson Method Matlab

# Newton Raphson Method Matlab

The aim of this session is to use a basic example to illustrate how to use the Newton Raphson method in Matlab.

We will make use of the ode45 solver, and use boundary conditions in the following example.

Here come the exercise

With boundary conditions

Here is the route we will take

## Newton Raphson method in Matlab

The code below solve this initial value problem (IVP) using the function ode45

it used the Newton Raphson method in the iteration process to approach the exact solution and finally end the iteration when y(1) is accurately converged up to the third decimal.

In the beginning of the problem we divide the ODE (ordinary differential equation) to a set of first order equations and we use 1 as initial guess for y'(0)

Here is the code

```ua=0;
s(1)=1; %First guess for the derivative
s(2)=1.1; %second guess for the derivative

tb=1;  %Second time
target_ub=0;  % target

f=@(t,y) [y(2); 1./(1+y(1))^2]; % set of 1st order ODE
rb=@(v) (v-target_ub)  % boundary conditions at b

[t,y]=ode45(f,[0 tb],[ua; s(1)]);
ub(1)=y(end,1);

for (j=2:50) % stop after 50 iteration even if the accuracy is not reached
[t,y]=ode45(f,[0 tb],[ua; s(j)]);
ub(j)=y(end,1);
if abs(rb(ub(j)))&lt;0.001 % quit the for loop if accuracy reached
break
end
s(j+1)=s(j)-rb(ub(j))*(s(j)-s(j-1))/(rb(ub(j))-rb(ub(j-1)));
end

plot(t,y(:,1)) % plot the last solution
xlabel('t')
ylabel('y')
grid on```

The Plot

#### Related Posts

This website uses cookies to improve your experience. We'll assume you're ok with this, but you can opt-out if you wish. Accept Read More