julia - Julia中的二阶微分方程

Question

我是 Julia 编程的新手，我设法解决了一些一阶 ODE，但是当我想转到二阶时，我不知道如何使用求解器来实现所需的方程。

我想解这个方程

y" + y = 0

有初始条件

y(0) = 3
y'(0) = -0.5

我怎样才能做到这一点？

Answer 1

将其更改为y'' = -y然后使用SecondOrderODEProblem，即https://diffeq.sciml.ai/stable/types/dynamical_types/，通过：

SecondOrderODEProblem((v,u,p,t)->-u,-0.5,3.0,(0.0,1.0))

有关这方面的更多示例，请参阅经典物理问题教程：

https://tutorials.sciml.ai/html/models/01-classical_physics.html