我们被问到以下问题： 编写伪代码以在 f 上实现 bisection 方法，该方法在 问题 1.您应该包括一个基于巴纳赫定理证明中使用的参数的停止标准。 问题 1 的给出如下：考虑函数 f（x） =...

来自 Jenkins-Traub 迭代的维基百科文章建议，这些步骤也可以通过三种不同形式的逆向量迭代来实现。 它建议一个矩阵来执行此迭代。 我的问题是，当根源未知时，如何创建这个矩阵。因此，难道不应该...

我正在研究一个用 JAX 编写的微分方程求解器。我遇到的一个常见工作流程是这样的： import jax.numpy as jnp from jax import jit # Function t...

def hermite_interpolation(x, y, yp, xi): n = len(x) - 1 result = 0.0 for j in range(n + 1): ter...

我想在scipy函数scipy.optimize.curve_fit（）背后获得一些直觉。下面，我写了一个简单的版本，试图最小化最小二乘函数的总和（权重都是 1），模型是 sin（omega*t）。我...

对于一个学校项目，我正在 Python 中实现一种有限差分方法，以数值求解以下反应扩散偏微分方程组： 到目前为止，我已经实现了以下目标： import numpy as np import mat...

我有一个目标函数，我想使用 MLE 来查找模型的参数。这是我到目前为止的代码。此外，我的returnarr是股票的日志返回数组 以下是我的 returnarr 的前 30 个值，仅供参考： retu...

错误是：- 程序接收信号 SIGSEGV：分段故障 - 无效的内存引用。 该程序用于梯形积分方法，其代码是 法典：- program tapizoidal implicit none real...

我有一个Matlab代码代码 model = createpde(); geometryFromEdges(model,@lshapeg); applyBoundaryCondition(mode...

我正在尝试评估（使用数值方法求解超越方程）当我改变初始猜测 M0 的值时，迭代次数会变化多少。我想我不应该期望它为如此不同的 M0 值返回相同数量的迭代，所以我认为这段代码可能不会使用我希望它使用的所...