我需要为一个函数实现一个根搜索算法，该算法可能有一些（不可移除的）分歧点（DP）/奇点，例如f(x)=x**2*np.tan(x)问题是，通常的条件f(x1)*f(x2)<0表示区间[x1,x2]内的根不能应用于不连续函数。如果对上述函数执行此操作，算法会将DP识

我正在尝试求解复杂的微分方程组。这些方程是刚性的，因此我需要使用一种可以处理复杂的ODE和刚度切换的方法。我已经使用odeintw(link)来解决我的问题，但是我还需要事件处理。我尝试使用solve_ivp和|来解决问题||，但是评估花费了不合理的时间。

我需要求解一些常微分方程$\frac{dy}{dx}=f(x)=x^2ln(x)$并继续在限制0之间创建数组xpt。<=xpt<=2因为我必须小心xpt=0，所以我将函数定义如下deff(x):ifx<=1.e-6:return0.else:returnnp.square(x)*np.log(x)我的调用程序读取Np