我需要为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程序。一个主程序，以及第二个牛顿雅可比脚本，因为我正在使用牛顿迭代方法，导致这个问题在你和你有两个边界条件
现在，我可以解决下一个问题：
y''(x) + k*y'(x)(y(x)**3/4)+(y'(x)**1/4)=0 ,
仅当k =常数！
在我的问题中，k依赖于x，在程序中计算，具有N个值。
有没有人知道或可以提出解决上述非线性系统的方法 我插入k，我试图在第二个牛顿求解器脚本中调用k_function（），但我没有成功。它没有认出来。
任何已经解决了非线性方程模拟系统的人的建议都将受到欢迎！