To get the solution for an ODE by ODE 45, loop can be used to extract and display the solution at each time step. Please refer to a simple code give below to get the statement:

clear all, close all ,clc

% Define the ODE function

ode = @(t,y)[y(1)*y(2); -y(1)];

% Set the initial conditions and time range

y0 = [1; 0];

tspan = [0, 100];

% Solve the ODE using ode45

[t,y] = ode45(ode, tspan, y0);

% Display the solution at each time step

for i = 1:length(t)

disp(['Time: ', num2str(t(i)), ', Solution: ', num2str(y(i,:))]);

end

% Set the pause time

h = 0.5;

% Display the solution at each time step

for ii = 1:length(t)

figure(1)

clf

plot(y(ii,1), 'ro')

xlabel('y(1)')

ylabel('y(2)')

grid on

hold on

pause(h)

end

thank you for your reply. In fact, i was not clear. the small script you suggest gives the solution after ode45 finishes.

However, i deal with FSI (Fluid Structure Interaction) i.e i should know the solution at each time step because i have to make change in the ODE

Welcome.