我正在尝试求解以下四阶边值问题
y“"”= K - C*y
我的x变量是一个有100个节点的线性空间。正如你所看到的，K是一个长度为100的向量，它使方程非齐次。但是，当我按下“求解”时，出现了以下错误：

Cell In [11], line 18, in fun(x, y)
     17 def fun(x, y):
---> 18     ans = vector-np.multiply(C,y[0])
     19     return np.vstack((y[1],y[2],y[3],ans))

ValueError: operands could not be broadcast together with shapes (100,) (99,)

为什么解算器突然将y的长度更改了1，我如何修复此错误？
编辑：我必须补充一点，当K不存在时，解算器工作得很好，即方程是齐次的。

from scipy.integrate import solve_bvp
import numpy as np

L = 10 
nodes = 100

A = 1000
B = 1500
C = 0.05

x = np.linspace(0,L,nodes)

vector = np.ones(nodes)

def fun(x, y):
    ans = vector-np.multiply(C,y[0])
    return np.vstack((y[1],y[2],y[3],ans))

def bc(ya, yb):
    return np.array([ya[2], yb[2], ya[3]+A/B, yb[3]])

y_a = np.zeros((4, x.size))

res_a = solve_bvp(fun, bc, x, y_a)

res1 = res_a.sol(x)[0]
res2 = res_a.sol(x)[1]
res3 = B*res_a.sol(x)[2]
res4 = B*res_a.sol(x)[3]