Problem with second order ODE solver
显示 更早的评论
Hello everybody,
I'm trying to solve the following equation: y" = sin(y) - (1/t)*y'
I write down the following code for solving the equation but the output does not show anything, do i miss something here?
syms y(t)
[V] = odeToVectorField(diff(y, 2) == sin(y)-(1/t)*diff(y))
M = matlabFunction(V,'vars', {'t','Y'})
sol = ode45(M,[0 30],[5.71 0])
fplot(@(x)deval(sol,x,1), [0, 30])
I really appreciate your help!
Thanks.
回答（1 个）
Star Strider
2017-5-9
The problem is having ‘t’ to include 0 and having ‘t’ in the denominator. This creates a ‘0/0’ condition that equates to NaN, and tha then propagates throughout the integration of your ODE.
The easiest way to avoid that problem is to ‘cheat’, and use eps instead of 0.
This works:
syms y(t)
[V] = odeToVectorField(diff(y, 2) == sin(y)-(1/t)*diff(y))
M = matlabFunction(V,'vars', {'t','Y'})
[T,Y] = ode45(M,[eps 30],[5.71 0]);
figure(1)
plot(T, Y)
grid
2 个评论
Star Strider
2017-5-10
My pleasure!
There are two outputs (columns) in ‘Y’, ‘Y(:,1)’ (the derivative) and ‘Y(:,2)’, the solved equation. If you only want the solved equation, the plot changes to:
figure(1)
plot(T, Y(:,2))
grid
类别
在 帮助中心 和 File Exchange 中查找有关 Numeric Solvers 的更多信息
2017-5-10
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
