通过数值逼近（与符号计算相反）解决数学问题的算法。
我正在尝试用Python解决一个非常简单的IVP，以便进行错误分析： 将 numpy 导入为 np 从 scipy.integrate 导入solve_ivp def dh_dt_zerodriver(t, h): 返回-2/小时 时间=50 圣...
在 Julia 中，我想编写一个函数，多次提示用户输入矩阵，然后将其输入存储到数组中。到目前为止我已经尝试过以下操作： 函数接受列表（...
我想生成一个矩阵，每一行和每一列中只有一个非零条目，其中每个非零条目的模数为一。朱莉娅有什么办法吗？
如何检测由于浮点加法和乘法舍入而导致的精度损失？ [已关闭]
从计算机系统：程序员的角度： 使用单精度浮点 表达式 (3.14f+1e10f)-1e10f 的计算结果为 0.0：值 3.14 由于舍入而丢失。 表达式 (
我想使用韦尔福德方法来计算运行方差和均值。我在 Python 中遇到了 Welford 方法的实现。然而，在测试时仔细检查它是否会导致......
Sympy 产生 `TypeError: unsupported operand type(s) for *: 'interval' and 'complex'` 对于复杂有理式表达式
我有这个sympy脚本： 从 sympy 导入 * 从 sympy.plotting 导入plot_implicit x, y = 符号('x y', real=True) 阿尔法 = sqrt(2) / 2 expr = 1 + ((1-alpha) * x + y*I) / (1 - alpha * (x + y*...
我的目标是创建一个 C++ 程序来计算几何序列的前三项。提供的是前三项及其乘法的总和。 假设答案...
我正在寻找一种为几何体创建六面体网格的方法，该几何体将通过Python迭代地稍微改变。同时，它会随着更多的迭代而发生显着变化。正如中提到的...
检查 scipy.solve_ivp() 中触发了哪个事件的快速方法是什么？ 作为上下文，我正在模拟双摆。我想检查第一根杆或第二根杆是否已经
我一直在研究一个涉及解决两对耦合颂的问题。建议我们使用 scipysolve_ivp 来执行此操作，但它似乎要么冻结要么运行 ext...
用 Numerov 方法求解一维薛定谔方程 (python)
我目前正在尝试求解一维薛定谔方程。 （与时间无关）使用 Numerov 方法。该方法的推导对我来说很清楚，但我在实现时遇到了一些问题。我尝试...
我正在尝试求解线性泊松方程，其具有齐次诺依曼边界条件，沿 y 方向的区间为 [-1,1]，沿 x 为周期性。最初，我测试了相同的代码......
我想编写一个C++程序来计算几何数列的前三项。给出了什么： 前三项之和， r 它们相乘的结果 老实说我不知道如何开始...
我想编写一个C++程序来计算几何数列的前三项。给出了什么： 前三项之和 它们相乘的结果 老实说，我不知道如何开始...
我正在尝试制作一个可以用于所有情况的休眠王子类型函数（2,3,4.....变量耦合微分方程）。我需要使这个功能适用于所有情况。 所以lotkav是...
所以我有 2 个二阶非线性 ODE，应用状态空间定理后我有 4 个一阶 ODE。 我正在尝试应用 RK4，但我认为我做错了，因为图表存在分歧。 我明白了...
我在MATLAB中有两个数值积分函数，如下： fun1 = @(x) log2(1+x).*(4*exp(2*lambda*(d*sqrt(- d^2 + (a./x).^(2/alp)))-(a ./x).^(2/alp).*... acos(((a./x).^(1/alp)./d).^(-1)))...
在 MATLAB 文档中，他们提供了一个针对具有一个参数 c 的单个变量进行数值积分的示例，如下所示： 乐趣 = @(x,c) 1./(x.^3-2*x-c); q = 积分(@(x)fun(x,5),0,2) 如果我想怎么办...
我试图创建一个等于卷积积分的函数 F(t)，然后通过取值 F(tn) 来计算 theta_n(n+1)。但是，我收到错误“未定义函数 'F' for in...
我创建了一个简单的 simpson_adaptive 方法，该方法使用我自己的 simpson 方法。 我的辛普森方法是正确的，但我的自适应方法似乎不适用于 积分( sin(2*pi*x)² ) 范围从...