【求助】我先求了常微分方程的数值解,怎样再求这个解的根
mathematica吧
全部回复
仅看楼主
level 1
kidukanai 楼主
sol = NDSolve[{2 y'[x] + x y''[x] == -x y[x]^3, y[0.0001] == 1,
y'[0.0001] == 0}, y[x], {x, 0.0001, 9}]%求得数值解
Plot[y[x] /. sol, {x, 0.001, 9}]%画图(和预想一样)
findroot[sol[x], {6, 7}]%这里想求根,根大概是x=6.9,就是上面图像与x轴交点,但是我写的代码没反应
请问应该怎么写求根的代码[泪]救救孩子吧,新人刚接触mathematic期末报告卡在这里了。
如果方便的话,我还想求这个图像与x轴交点处的导数,麻烦了先谢谢啦
问题的背景(虽然估计没啥用)是n=3的Lane-Emden方程,用这个求白矮星的Chandrasekhar Limit质量。
2021年02月08日 20点02分 1
level 9
2021年02月10日 10点02分 3
level 9
2021年02月10日 10点02分 4
明白了,十分感谢🙏
2021年02月10日 10点02分
level 9
上面只是简略的替换写法,其本质是
2021年02月10日 10点02分 5
吧务
level 10
Plot的时候还知道ReplaceAll,咋findroot那行就不用了呢
2021年02月13日 03点02分 6
1