我用MatLab编写了以下代码，

n = 5;
ainit = zeros(1,n); 
    for i=1:ndof
            ainit(i) = 0; 
    end
    
binit = zeros(1,ndof); 
disp(binit)
x0 = [ainit,binit];

tf = 10; step = .001; tspan = 0:step:tf
funky=@(t,x) A*x + BF; %beam(t,x,K,M,wns,psi2,Load);
[t2,q2] = ode23(funky,tspan, x0, []);

其中A是n x n方阵，BF是n x 1向量
我正在尝试将这段代码转换为Python。到目前为止，我得到的信息如下：

y0 = np.zeros((4*n))
Time_Step = np.arange(0, 10, 0.001)
func = lambda x, t: (A*x+BF)
q = scipy.integrate.solve_ivp(func, Time_Step, y0)

这给了我一个错误，说明could not broadcast input array from shape (20,20) into shape (20,)
如何将MatLab的ode23转换为Python？