请教求解二阶ODE方程# Computation - 科学计算
t*o
1 楼
y'' + f(x) y' + g(x) = 0
以前一直用4阶的R-K方法
z1 = y
z2 = y'
原方程等价于方程组
z1' = z2
z2' = -f(x)*z2 - g(x)
但现在要解类似 xy'' + f(x) y' + g(x) = 0
这样的话
z2' = (-f(x)*z2 - g(x))/x
当 x = 0时， 就会出问题
请教如何解决？ 谢谢！！！
另外求解此方程需要两个边界条件y(a) y'(a)， 如果不知道边界条件， 尤其是y'(a)，
有什么办法估算初值吗？
thanks a lot!!!
以前一直用4阶的R-K方法
z1 = y
z2 = y'
原方程等价于方程组
z1' = z2
z2' = -f(x)*z2 - g(x)
但现在要解类似 xy'' + f(x) y' + g(x) = 0
这样的话
z2' = (-f(x)*z2 - g(x))/x
当 x = 0时， 就会出问题
请教如何解决？ 谢谢！！！
另外求解此方程需要两个边界条件y(a) y'(a)， 如果不知道边界条件， 尤其是y'(a)，
有什么办法估算初值吗？
thanks a lot!!!