我需要为 y 解决下一个系统:
y''(x) + k(x)y'(x)(y(x)**3/4)+(y'(x)**1/4)=0
是这样的:
i=1: y''(1) + k(1)y'(1)(y(1)**3/4)+(y'(1)**1/4)=0
i=2: y''(2) + k(2)y'(2)(y(2)**3/4)+(y'(2)**1/4)=0
....
i=N: y''(N) + k(N)y'(N)(y(N)**3/4)+(y'(N)**1/4)=0
我基于两个脚本编写了一个 Python 程序。一个主程序和第二个 newton-jacobian 脚本,因为我使用的是牛顿迭代法,导致这个问题在 u 和 u' 处有 2 个边界条件
现在,我可以解决下一个问题:
y''(x) + k*y'(x)(y(x)**3/4)+(y'(x)**1/4)=0 ,
仅当 k=常数时!
在我的问题中,k 取决于 x,在程序中计算,有 N 个值。
有谁知道或可以建议一种方法来解决我插入 k 的上述非线性系统,并且我尝试在第二个牛顿求解器脚本中调用 k_function(),但我没有成功。它不承认它。
任何解决非线性方程模拟系统的人的建议都将受到欢迎!