y" + y = 0
с начальными условиями
y(0) = 3 y'(0) = -0.5
Как я могу это сделать?
измените его на y'' = -y, а затем используйте SecondOrderODEProblem, т.е. https://diffeq.sciml.ai/stable/types/dynamical_types/, через:
y'' = -y
SecondOrderODEProblem
SecondOrderODEProblem((v,u,p,t)->-u,-0.5,3.0,(0.0,1.0))
Дополнительные примеры см. В руководстве по классическим физическим проблемам. :
https://tutorials.sciml.ai/html/models/01-classical_physics.html