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)))<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