how can I plot this differential equation?
显示 更早的评论
I have a complicated differential equation that cannot have explicit solution for m>=3. I am trying to plot the solution so that I can use the data.
%parameters
ed=0.001;
Ev=1000;
n=4000;
m=3;
ed=0.001;
syms y(t)
cond = y(0) == 0;
eq= diff(y,t) == (1/m)*(((Ev*(ed*t-y))/n)^m);
采纳的回答
Star Strider
2018-8-6
编辑：Star Strider
2018-8-6
Add this to the end of your posted code:
[VF,Subs] = odeToVectorField(eq);
DE = matlabFunction(VF, 'Vars',{t,Y});
[t,y] = ode45(DE, [0 10], 0);
figure(1)
plot(t, y)
grid
EDIT — The complete code then becomes:
syms y(t) Y t
%parameters
ed=0.001;
Ev=1000;
n=4000;
m=3;
ed=0.001;
cond = y(0) == 0;
eq = diff(y,t) == (1/m)*(((Ev*(ed*t-y))/n)^m);
[VF,Subs] = odeToVectorField(eq);
DE = matlabFunction(VF, 'Vars',{t,Y});
[t,y] = ode45(DE, [0 10], 0);
figure(1)
plot(t, y)
grid
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
