How to use ode45 on a second-order ODE, but without a function file (just one script)?
2 ビュー (過去 30 日間)
古いコメントを表示
I am given the following second-order ODE and I want to use ode45 and by not creating a function file.
![Capture.PNG](https://www.mathworks.com/matlabcentral/answers/uploaded_files/241590/Capture.png)
I am given the equation above, a timespan t, and initial y and dy values.
I first set up my dydt into a matrix [y1' y2']"
dydt = [y(2); eps*((y(1)^2) - 1)*y(2) - y(1)]
func= dydt(t,y)
and tried to evaluate such ode45(@func, t, [y0, dy0]), but it gives me an error
I also tried setting it up as just one function:
f = @(y1,y2) eps*((y1^2) - 1)*y2 - y1;
Here, it gives me the error that I am using the same function, which was previously used as a variable:
ode45(@f, t, [y0, dy0])
Help please?
0 件のコメント
採用された回答
Star Strider
2019 年 10 月 8 日
You are almost there!
Try this:
dydt = @(t,y) [y(2); eps*((y(1)^2) - 1)*y(2) - y(1)];
ic = [1 1]*eps;
tspan = [0 10];
[T,Y] = ode45(dydt, tspan, ic);
figure
plot(T,Y)
grid
That worked when I tried it. Change the initial conditions vector ‘ic’ to match your initial conditions, and ‘tspan’ to get the result you want.
4 件のコメント
その他の回答 (0 件)
参考
製品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!