闭。此问题需要调试详细信息。它目前不接受答案。 编辑问题以包括所需的行为、特定问题或错误以及重现问题所需的最短代码。这将帮助其他人回答这个问题。 7天前关闭。 改进此问题 我正在尝试编写代码以...

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

我正在为我的数值方法科目做一项任务，他们要求使用 Runge-Kutta 和 Taylor 方法来估计微分系统（Chua 电路）的解，其场由下式确定： F(x,y,z)=\[2\*(y-0.5\*x...

对于一个学校项目，我正在 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 值返回相同数量的迭代，所以我认为这段代码可能不会使用我希望它使用的所...

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

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