Solution to ODE in matlab or simulink

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.

Thanks in advance.

[Sulaymon Eshkabilov]

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]

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.

Answer 1

[darova]

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;
    [dnt, dnL] = gradient(n,dt,dL);
    n = -V/Q*(dnt+G*dnL);
    surf(t,L,n)
    pause(0.1)
end

The result