通过数值逼近（与符号计算相反）解决数学问题的算法。
在我目前正在阅读的书中，确定性数值方法（例如基于梯度的：牛顿方法，无梯度：Nelder-Mead）和算法方法（例如...
我正在尝试求解包括空气阻力在内的自由落体的常微分方程。 因此我将 ODE 定义为： def f(v, g, k, m): 返回 g - k/m * v**2 我认为这应该代表系统
我有两个 3D 点，例如 a = (100, 100, 10) 和 b = (0, 100, 60)，并且想通过这些点拟合一条线。 我知道，3D 直线方程可以有不同的形状： 矢量形式： ...
尊敬的各位专家、各位老同事： 我想使用提供的代码创建下图。然而，这似乎花费了很长的时间，但没有明显的结果。 这是我的
k=3;m=123 欧拉C_x=[A] eulerC_v=[0] eulerC_t=[0] c=0 对于 np.arange(t_euler,t_max,dt_euler) 中的 z： a=-(k/m)*欧拉C_x[c] eulerC_v.append(eulerC_v[c]+a*dt_euler) eulerC_x.append(eulerC_x[c]+
我需要一个服从麦克斯韦分布的 10,000 个随机数。 对于正态分布，我知道我必须使用 Box-Muller 变换，但我的问题是为什么 正常分布...
我正在寻找一个R代码，它可以帮助我在以数字形式给出pdf时找到分位数。也就是说，假设我的数据是 x = c(0.00,0.05,0.10,0.15,0.20,0.25,0.30,0.35,0.40,0.45,0.50,0.55,0.60,0.65...
为了开始使用 CORDIC for log10，我实现了此 PDF 第 6 页中派生的算法： #包括 #包括 // https://www.mikrocontroller.net/attachment/...
我想计算形状为 (Nmesh, Nmesh, Nmesh) 的 3D 笛卡尔网格的一阶偏导数，例如按四阶精度方案Nmesh=512，即 f'(n) = 2*(f(n+1)-f(n-1))/(3*dh...
我有格式化为一组 x,y 数据的真实数据。我想对这些数据进行梅林变换并找到梅林图像。 有一个python包sumpy积分可以使用...
这是两个非线性方程组： 在对上一个问题的回答中，已经表明该方程组不具有显式解（此处）。但也值得注意的是
我求解金属棒的热方程，一端保持在 100 °C，另一端保持在 0 °C，如下 将 numpy 导入为 np 将 matplotlib.pyplot 导入为 plt dt = 0.0005 dy = 0.0005 k = 10**(-4) y_最大值 = 0.04
我想知道如何使用 CUDA 生成两个以上列表的笛卡尔积。 如何使此代码适用于三个或更多列表？ 它适用于两个列表，但不适用于三个列表，我尝试过......
我需要求解以下非线性方程组： 我正在使用 sympy 的solve()函数来获得方程组的解。通常，solve() 会产生所有解决方案。我...
有人可以找到我的代码中的错误吗？它在特殊点上工作正常，但公差不适合该方法。我的 Errorstepper 非常简单，所以我不认为其中有问题。拜...
我正在使用 JiTCDDE 模拟时滞 Kuramoto 模型，并且使用 get_state() 通过时间演化提取导数时遇到一些问题。 我正在为延时的 Kuma 建模...
我正在尝试使用二分法求解方程。但是，当我尝试运行它时，出现以下错误 “if ((fn(kVec, tVec, b) * fn(kVec, tVec, a) > 0) | (b > a)) { 中的错误：
我有一个（也许不是那么聪明）的想法，应用浅水模型来计算流域的全年流域径流。这意味着使用具有可变坡度和粗糙度的真实地形......
是否存在 Fortunes eigensolve 算法的实现？ 目前我正在尝试将其移植到 python，但到目前为止它没有提供正确的结果。 我唯一发现的是...
目前我的目标函数给出的解大小约为 1e-6。我认为这给我带来了一些精度问题。 处理此类问题的最佳实践是什么？看来...