How to solve this ODE? Urgent help!!

1 次查看(过去 30 天)
Please explain step wise and also I want a plot of m vs. t t varies from [0,0.5] and m(0)=0.4
Thanks

采纳的回答

Star Strider
Star Strider 2014-7-26
This is similar to your previous code, this time with the additional constant (that does not seem to affect the result). There is not much to describe, really. See the documentation on Anonymous Functions and ode45 for those details.
Otherwise:
mdot = @(t,m) -0.31458*sqrt(2E+5.*((1+(0.4-m)./0.5).^-1.4 - 1)); % Anonymous function encoding ODE
[t,m] = ode45(mdot, [0 0.5], 0.4); % Integrate using ode45
figure(1)
plot(t, m, '-*b') % Plot
grid
The ode45 call has as its arguments, the function to integrate ( mdot here ), the time span ( [0 0.5] here ), and the initial conditions ( 0.4 here, and only one because there is only one variable to integrate ). The figure plots m as a function of t plotting m(t) with a blue asterisk and drawing a blue line between them. The grid call adds grid lines to the plot.

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Ordinary Differential Equations 的更多信息

产品

Community Treasure Hunt

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

Start Hunting!

Translated by