通过数值逼近（与符号计算相反）解决数学问题的算法。
我有一个名为 trap_1D 的 Python 函数，它使用复合梯形规则计算单变量函数 f 在区间 [a,b] 上的（近似）积分。现在我想评估一个
我尝试使用牛顿法近似 255 / sqrt(x) 以避免在架构上使用除法或 sqrt() 操作，因为这些操作的成本很高。 我的推导是： y = 255 / 秒...
我这样做是为了使用射击方法解决BVP功能。 我不明白为什么当我的步长 h = 0.1 时，限制 x_end 超出了我在函数solve下的 while 循环中设置的值...
我正在编写代码，我有一个函数可以用满足条件的值来计算不满足条件的值，但是我在管理方面遇到了很多麻烦...
我正在编写代码，我有一个函数可以用满足条件的值来计算不满足条件的值，但是我在管理方面遇到了很多麻烦...
我正在实现隐马尔可夫模型的前向算法（算法见下文）。为了防止上溢/下溢，我使用对数概率，并使用对数和表达式技巧来配合...
我正在尝试使用 SciPy 优化来找到最小化两点之间悬挂绳索势能的函数，即悬链线问题。目标函数定义为 那个...
从十进制数计算分数的算法是已知的，并且有足够的例子。就我而言，我需要找到给定小数的分数，其中 A = N/D 已优化，但是
我正在使用 scipy 的solve_bvp 求解常微分方程，以获得一系列特征态 phi_n 和特征值 omega_n。等式看起来像 y'' = -(omega_n - f(x)) ** 2 * y...
我正在尝试用 Javascript 对弹簧上的质量进行 Runge-Kutta 实现，并使用 D3 对其进行可视化。目的是将其与前向欧拉进行比较并评论差异。我的有限元工作...
我正在寻找一些四边形网格生成的示例代码（任何语言）。然而，这似乎是一项相当艰巨的任务！ 我不挑剔，我想至少对有孔的多边形进行网格划分，什么都没有......
迭代求解器比基于 LU 分解的直接求解器更稳定吗？对于基于 LU 的求解器，我们总是有 cond(A) < cond(L) * cond(U), so factorization amplifies numerical inaccuracy. So...
Burden & Faires 中的浮点问题以及高斯消元法示例
我正在阅读 Burden & Faires 第 9 版第 360 页第 6.2 节。它包含一个示例，旨在说明在没有行交换的情况下高斯消除的数值不稳定性......
有没有一种简单的方法可以判断 varchar 是否是数字？ 示例： abc123 --> 没有号码 123 --> 是的，它是一个数字
有没有一种简单的方法可以判断 varchar 是否是数字？ 示例： abc123 --> 没有号码 123 --> 是的，它是一个数字
我目前正在努力绘制 R 中的下一个函数。 (F1 / F2^m * (R1 - F1)^((1 - s) / s)) - (m * (((R1 - F1))^(1 / s) + ((R2 - F2))^ (1 / 秒)) / (F1^m + F2^m)) = 0 我想将 F2 绘制为
如果我们计算像 A=[1 2 3;4 5 6;7 8 9] 这样的秩亏矩阵的 QR 分解，R 矩阵的一些对角项将非常小。然而，当在 Ju 中进行数值计算时...
使用`solve_ivp`和`LSODA`来求解复杂的ODE
我正在尝试求解复杂的微分方程组。这些方程是刚性的，因此我需要使用一种可以处理复杂的 ODE 和刚度切换的方法。我已经降落在
我需要求解一些常微分方程 $ rac{d y}{d x} = f(x) = x^2 ln(x)$ 并继续在限制 0 之间创建数组 xpt。 <= xpt <= 2. Because I have to be careful ...
我试图用 C 级 LAPACK 替换我的 cython 代码中的 Scipy 的 eigh 例程。本质上我想计算对应于埃尔米特 m 的最大和最小特征值的特征向量...