在python中使用scipy.integrate.solve_ivp时，是否有一种方法来绑定变量的值？

问题描述

我想用两个变量x和u来解决python中的IVP，但是我需要u的值在0到1之间。现在，它为我提供了一个负值的解决方案。这是我的代码。

def model_case_3(t,z,Kmax,k,b,list_Kmax,sigma):
 a=1.5
 x,u= z
 m=1
 r = list_Kmax[0][0]
 dxdt =  (x)*(r*(1-a*u**2)*(1-x/(Kmax*(1-0.999*u**2)))-m/(k+b*u)-0.05)
 dudt = sigma*((-2*a*(b**2)*r*(u**3)+4*a*b*k*r*(u**2)+2*a*(k**2)*r*u-b*m)/((b*u+k)**2))
 return [dxdt,dudt]

 sol = solve_ivp(fun=model_case_3,t_span=[scaled_days[i][j][0],scaled_days[i][j][-1]],y0=[scaled_pop[i][j][0],list_u[0][0][0]],t_eval=scaled_days[i][j],args=(list_Kmax[0][0],k0,b0,list_b,sigma0))

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）