MATLAB: Ode45 third order ode

ode45

how to solve

f''' = { [3 * f' * (f'')^2] / [(f')^2 + 1]^(5/2) + 1/f^3 – 1/f^2 + 3} * { [(f')^2 + 1]^(3/2) }

using ode45

with

f(0) = 1.1

f'(0) = 17.1

f''(0) = 144.1

Best Answer

Use ode45(). this ODE can be written as a system of 3 first-order ODEs

odeFun = @(t, y) [y(2);
                  y(3);
                  ((3*y(2).*y(3).^2)./(y(2).^2 + 1).^(5/2) + 1./y(1).^3 - 1/y(1).^2 + 3).*((y(2).^2 + 1).^(3/2))];
              
tspan = [0 1];
ic = [1.1; 17.1; 144.1];
[t, y] = ode45(odeFun, tspan, ic);
plot(t, y);

However, it seems that the ODE is unstable, and the solution diverges to infinity. You may check if the equation is written correctly.

Related Solutions

MATLAB: F=@(t,y) 4exp(0.8t)-0.5*y; tspan = [0 4]; y0 = 2; h = 1; [t,y] = ode_euler(f,tspan,y0,h)

tspan = [0 4];
y0 = 2;
[t1,y1] = ode45(@(t,y) 4*exp(0.8*t)-0.5*y,tspan,y0,h);
figure('name','Reda','color','w');
plot(t1,y1, 'r-','LineWidth',2);grid on;
xlabel('Time');ylabel('Solution y(Time)');
title('And you instead of a blanket substitute solver Euler');

MATLAB: 4th order Runge-Kutta code that can solve for several intial conditions

Here's a version:

% Initial conditions
t0=0;
x0=0:0.2:1;
y0=1:0.2:2;
dt=0.1;
tz=1;
t_range= t0:dt:tz;
X = zeros(numel(x0),numel(t_range));
Y = zeros(numel(y0),numel(t_range));
x_rk=zeros(1,numel(t_range));
y_rk= zeros(1,numel(t_range));
for n = 1:numel(x0)
    x_rk(1)=x0(n);
    y_rk(1)=y0(n);
    for i=1:length(t_range)-1
        [x_rk(i+1), y_rk(i+1)]=rk_method(x_rk(i),y_rk(i),dt);
    end
    X(n,:) = x_rk;
    Y(n,:) = y_rk;
    %% PLots
    subplot(numel(x0)/2,2,n)
    plot (t_range, x_rk,'b-o',t_range, y_rk,'y-o','MarkerSize',5);
    grid
    title([num2str(x0(n)),'   ',num2str(y0(n))])
    legend('x','y')
end
disp('X = '), disp(X)
disp('Y = '), disp(Y)
%% Functions
function [dxdt, dydt]= f(x,y)
dxdt =  5*x-3*y;
dydt = -6*x+2*y;
end 
function [x1, y1]=rk_method(x0,y0,dt)
[k11, k12]=f(x0,y0);     
[k21, k22]=f(x0+0.5*dt*k11,y0+0.5*dt*k12); 
[k31, k32]=f(x0+0.5*dt*k21,y0+0.5*dt*k22); 
[k41, k42]=f(x0+dt*k31, y0+dt*k32); 
x1=x0+dt*((k11/6)+((k21+k31)/3)+(k41/6));
y1=y0+dt*((k12/6)+((k22+k32)/3)+(k42/6));
end

Best Answer

Related Solutions

MATLAB: F=@(t,y) 4*exp(0.8*t)-0.5*y; tspan = [0 4]; y0 = 2; h = 1; [t,y] = ode_euler(f,tspan,y0,h)

MATLAB: 4th order Runge-Kutta code that can solve for several intial conditions

Related Question

MATLAB: F=@(t,y) 4exp(0.8t)-0.5*y; tspan = [0 4]; y0 = 2; h = 1; [t,y] = ode_euler(f,tspan,y0,h)