如何避免随着时间的流逝和“漂移”的累积错误在数字时间集成中？
我正在编码一些研究结构动态问题的程序，并且我正在尝试了解如何随着时间的流逝估计和最大程度地减少累积错误，以进行大规模计算。 一个显式时间集成…
if statements and integate（）在r中函数
如果语句不起作用，请执行吗？我必须做比这更复杂的事情，但是我提供了这个示例，因为它隔离了问题。 Kernel = function(x){ if(abs(x)<1){ w = 1 - a…
如何在R中的网格点上通过函数值给出的任意区域执行2D数值集成？
我想数字评估形式的积分： $$ int_d f（x，y）dx dy $$， 其中$ d $是$ r^2 $和$ f（x，y）$的区域从$ r^2 $到$ r $。 我有一个网格覆盖$ r^2 $的子集…
Scala：递归和堆栈溢出错误
我正在尝试进行以下练习： 定义递归函数 intedef ，该函数 [a ，B] 使用复合梯形规则。该函数将以 eps 为参数，该参数标记了将间隔 [a，b] 划分的子间…
如何在Python中定义一个函数，该函数包含一个确定的积分，其自由参数在两个固定值之间运行？
我想定义一个具有三个参数的Python函数：C，PSI和M，但实际上它具有基于另一个称为ETA的变量的积分。 def func(x: np.float64, C: np.float64, psi: n…
加速 scipy.nquad 计算多重积分
我试图计算这个函数： 这是多重积分的结果。 函数 F_A_、F_B_、F_C_、F_D_ 和 F_FPS_ 是截断的正态函数，包含在设置的间隔之间。我用 scipy.nquad 尝…
我正在尝试计算八度上的二重积分
所以我需要评估八度音程上的二重积分[![question][1]][1] >> x1 = 1; >> x2 = 0; >> y1 = 2; >> y2 = 0; >> f = @(x,y) (y-x)./((x+y).^3); >> integr…
如何从 np.array 获取积分值？
我正在尝试从 np.array 值列表中获取积分值。不是函数下的表面，而是值。我有加速度值并且想要获得速度值。 假设我有一个像这样的数组： a_x = np.arr…
scipy.integrate 和 mpmath.quad 的积分结果不同
我正在集成一个函数来实现本文中的功能，并且我观察到不同的scipy.integrate.quad 和 mpmath.quad 的输出如下： import mpmath as mp from scipy impo…
Scipy 集成返回负值
我有一些 xdata 和 ydata。我已经对 xdata 和 ydata 进行了排序，使得 ydata 不断增加（即 ydata 数组末尾的最高 y 值）。 ydata 还包含一些点 = 0，…
Python：求解自治二阶 ODE
我有以下自治常微分方程系统，其中x和v是向量（初始条件x0和 >v0) dx/dt = v dv/dt = F(x, v) 等价地，可以将其编写为二阶 ODE 系统 d^2 x/ dt^2 = F(…
有没有办法对 NumPy 数组中具有相同值的所有行应用函数？
假设我们有一个矩阵 A，它具有以下值： In [2]: A Out[2]: array([[1, 1, 3], [1, 1, 5], [1, 1, 7], [1, 2, 3], [1, 2, 9], [2, 1, 5], [2, 2, 1], […
尝试将表达式与 lambda 集成时出现问题
我正在尝试集成一个具有实数和复数值的表达式，将其定义为 lambda 表达式。积分变量是 kx，积分的结果解将在 x 和 y 维度上求值，但是在我积分并尝试…
