使用近似技术在一个或多个维度上集成函数的算法，而不是使用符号代数和微积分的精确闭合形式解决方案。包括自适应求积法，蒙特卡罗方法，有限元分析，马尔可夫链等概念。
如何解决scipy.integrate.quad中特定被积函数的舍入误差问题？
我正在做下面的数值计算。 将 numpy 导入为 np 从 scipy.integrate 导入四边形 def E_k（间隙，q，mu）： return np.sqrt(abs(gap) * abs(gap) + (q * q - mu) * (q * q - mu)) 四边形（...
我必须在 3D 空间上进行大量集成。我对所述空间进行了（网格）离散化（例如，在每个方向上等距 30 个点），我知道我想要的函数值
我想定义一个函数，它接受一定长度的数组输入，然后为每个值找到一个值的最佳参数。 它看起来像这样： 这里 a, b, c 和 measu...
我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长，t= np。排列 (0, 3600,10)。这是正确的吗？
我想通过求解运动方程来跟踪一个小时的粒子路径。对于时间步长， t= np。排列 (0, 3600,10)。 它是否正确？ 持续时间是否取决于时间步长。如果我...
我正在尝试模拟流体（气体或液体，例如空气或水）中的可操纵车辆（飞机、气垫船、船）。 我的 JavaScript 实现是基于 3D 刚体飞机模拟...
我正在尝试加速以下计算代码： 我只需要为 x > y 从 0 到 1 计算此函数（但需要非常高的离散化，如 dt = 0.001）。我已经矢量化了...
odeint -(TypeError: can't convert expression to float ) - 将表达式导入到函数中执行 odeint。
这是我得到的错误 回溯（最近一次调用）。文件 "C:\Users\user.spyder-py3\Numerical Methods Problems\FreeFall.py"，第40行，在 "C:\Users\user.spyder-py3 "中。 ans=odeint(vel,0,t) 文件 "C:\..."。
在使用C语言的函数时向scipy.LowLevelCallable传递参数。
我试图在SciPy中使用C定义的函数进行数值积分。这里给出的例子(SciPy文档)可以正常工作。在我的例子中，testlib.c文件是 * testlib.c * #include &...。
为什么这个太阳系的数值积分一直在运行？(MIT-scme SCMUTILS)
我想做一个太阳系的数值积分。我以前用普通的Scheme做过，现在我想用麻省理工学院非常有趣的SCMUTILS-library来做。我做了什么。我把太阳能...
经过大量的编辑，研究，我决定做正确的（我希望）。所以，我想从命令行中得到一个浮点数，下面是代码： format PE Console entry start include 'INCLUDEwin32ax......'。下面是代码： format PE Console entry start include 'INCLUDEwin32ax......'。
我用MATLAB和这个函数：函数[x,w,P]=lglnodes(N) % Truncation + 1 N1=N+1; x=cos(pi*(0:N)N)'; % Legendre Vandermonde Matrix P=zeros(N1,N1); xold=2; while max(abs(x-xold))>eps ...
我想用sympy解一个积分。但它给了我一个错误的解决方案，为什么？为什么？ import sympy from sympy import Integral, exp, oo x, y = sympy.symbols("x y", real=True) b, u, l, t = sympy.symbols("...")
我有一个RKF7（8）积分器，我已经用几个简单的测试功能验证了其输出。但是，当我在感兴趣的运动方程式上使用它时，局部截断误差为...
我有一个有效的Integrator类，它将计算单个变量的基本函数的确定积分。我已经测试了一些基本功能的集成，并且似乎可以正常使用...
我正在实现RKF4（5）集成器，但我无法弄清楚我的代码是否有效，并且我不了解本地截断错误，或者我的代码是否无效。我为代码的大小表示歉意...
为什么用Python重写Simpson的积分规则会得到不同的结果？
以下源代码是辛普森积分规则的实现。 C＃源代码。使用系统；公共类Simpson {private double Function（double x）{return 1 ....
[我想知道是否有人知道基于numpy / scipy的python软件包，以数字方式将复杂的数值函数集成在细分的域上（在我的特定情况下，二维域以...为边界）]
我在编程方面很烂，所以请多多包涵。我什至不确定我现在需要的概念是什么，所以我不知道要用谷歌搜索或写这篇文章的标题。我的问题是，...
因此，我试图使用Bulirsch-Stoer自适应方法来解决EDO，但是却收到了溢出警告。如何解决呢？从数学导入sin，cos，pi从numpy导入范围，数组，空，...
并希望每个人都安全。我想使用Simpson规则来评估以下双重积分：通常，如果两个限制都是数字，那么我只需要从scipy中应用simpy就可以避免....