我正在用python求解微分方程4（y '）^3-y'= 1/x^2，我熟悉用odeint求解耦合常微分方程和线性常微分方程，但在非线性常微分方程上找不到太多指导，比如我正在努力解决的问题。尝试使用odeint和scipy，但似乎无法正确实现任何想法都非常感谢注意：y是x的函数
7fyelxc51#
问题是，在相空间的每一点（包括双根），你都可以得到3个有效的方向解，但是每个选择准则在双根处都失效了。一种方法是在系统y'=v, 4v^3-v=x^-2上使用DAE求解程序（scipy中不存在第二种方法是对方程求导，得到显式二阶常微分方程y''=-2/x^3/(12*y'^2-1)。两种方法都需要从初始点处的三次根中选择初始方向。
y'=v, 4v^3-v=x^-2
y''=-2/x^3/(12*y'^2-1)
1条答案
按热度按时间7fyelxc51#
问题是，在相空间的每一点（包括双根），你都可以得到3个有效的方向解，但是每个选择准则在双根处都失效了。
一种方法是在系统
y'=v, 4v^3-v=x^-2上使用DAE求解程序（scipy中不存在第二种方法是对方程求导，得到显式二阶常微分方程
y''=-2/x^3/(12*y'^2-1)。两种方法都需要从初始点处的三次根中选择初始方向。