python scipy库函数solve用法_对带有额外参数的函数使用scipy solve_ivp

我试图用solve-ivp解决一个初始值问题。问题是我的函数f(x，y)=dy/dx包含额外的参数。我试着用这个：

但它总是给我错误。代码如下：from numpy import arange

import math

import numpy as np

from scipy.integrate import solve_ivp

from numpy import pi

sigmav = 1.0e-9

Mpl = 1.22e19

ms=100.0

gg=100.0

gs=106.75

def Yeq(x):

return 0.145*(gg/gs)*(x)**(3/2)*np.exp(-x)

def ss(x,m_dm):

return (2/45)*((pi)**2)*gs*(m_dm**3/x**3)

def hubb(x,m_dm):

return 1.66*(gg**(1/2))*(m_dm**2/Mpl)*(1/x**2)

def derivada(x,yl,m_dm,σv):

return -(σv*ss(x,m_dm)*((yl**2)-(Yeq(x)**2)))/(x*hubb(x,m_dm))

xx=np.logspace(np.log10(3),3,100)

y00 = Yeq(xx[0])

x00 = xx[0]

sigmav = 1.7475e-9

ms=100.0

solucion1 = solve_ivp(fun=lambda x, y: derivada(x, y, args=(ms,sigmav),[np.min(xx),np.max(xx)],[y00],method='BDF',rtol=1e-10,atol=1e-10))

当我尝试为solucion1运行单元时，它返回以下结果：

^{pr2}$

这一定很简单，但我刚刚开始使用python，所以请帮助我。在