朱莉娅的二阶微分方程

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