0

我需要为 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(),但我没有成功。它不承认它。

任何解决非线性方程模拟系统的人的建议都将受到欢迎!

4

0 回答 0