How to solve for the derivative using ode solver

Question

Noya Linder 2023-7-12

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1994958-how-to-solve-for-the-derivative-using-ode-solver

评论： Noya Linder 2023-7-12

采纳的回答： Jayant

I have the following differential equation

function dRdt = odefun(t,R)

c = 29979245800;

Estart = 10^52;

ni = 10^(-2);

mp = 1.67262192*10^(-24);

T3 = 17*Estart/(8*pi*ni*mp*c^2);

dRdt = sqrt((T3-R^3)*c^2/T3);

end

and I solve it using ode78

[t, R, te, Re, ie] = ode78(@odefun, tspan, r0, options)

thus obtaining R(t). I want to solve the same equation but instead for the derivative (I want to get an array of dRdt and it's corresponding t). How should I go about that?

Thank you in advance!

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jayant 2023-7-12

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1994958-how-to-solve-for-the-derivative-using-ode-solver#answer_1271918

在 MATLAB Online 中打开

Hello @Noya Linder,

To obtain an array of dRdt and its corresponding t values correctly, you can modify the code as follows:

function dRdt = odefun(t, R)
    c = 29979245800;
    Estart = 10^52;
    ni = 10^(-2);
    mp = 1.67262192*10^(-24);
    T3 = 17*Estart/(8*pi*ni*mp*c^2);
    dRdt = sqrt((T3-R^3)*c^2/T3);
end
tspan = [t_start, t_end];
r0 = initial_condition;
options = [];
% Solve the differential equation
[t, R] = ode78(@odefun, tspan, r0, options);
% Calculate dRdt for each time point
dRdt = arrayfun(@odefun, t, R);

This modified code will correctly calculate dRdt for each time point by using the arrayfun function to apply the odefun to each t and R value obtained from the ode78 solver.

You may refer this documentation of arrayfun for reference.

Hope that helps.