# Solution to ODE in matlab or simulink

3 views (last 30 days)
Soheil on 4 Oct 2019
Answered: darova on 5 Oct 2019
Hi,
I would be thankful if someone gives a comment on how to solve the following equation in MATLAB/Simulink with ODE45 or ODE15s...It's equation number 3 in this publication [Link]. (Picture of the equation is attached).
It's basically a material balance equation of solid particles in ractor. Sulaymon Eshkabilov on 5 Oct 2019
Hi
Your equation as it does not appear to be correct. Could you please post the correct one?
Just a guess as it looks like a 1st order ODE BVP that means you'd need to employ bvp solver
Soheil on 5 Oct 2019
Thanks for the comment; I will check the BVP solver.
I post a picture of equations from a paper below.
They've mentioned at the end that ODE45 is used.
A general guideline toward solving this would be super helpful. darova on 5 Oct 2019
I tried to solve this equation first What are constants values?
% opengl software
G = 2;
V = 2.4;
Q = 5;
J = 1;
n0 = 1.1;
num1 = 30;
num2 = 40;
n = zeros(num1,num2);
L = linspace(0,30,num1);
t = linspace(0,100,num2);
dL = L(2) - L(1);
dt = t(2) - t(1);
for i = 1:100
n(end,:) = J/G;
n(:,1) = n0;
n = -V/Q*(dnt+G*dnL);
surf(t,L,n)
pause(0.1)
end
The result R2019b

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!