通过数值逼近（与符号计算相反）解决数学问题的算法。
我试图应用scipy的solve_bvp来解决以下问题 T'''''(z) = -k^4 * T(z) 在大小为l的域和一些常数A的边界条件下： T(0) = T'''(0) = T''''(l) = 0 T'(l) = A 所以 ...
我试图用欧拉方法找到微分方程在某点的解，但我得到的答案是错误的。请看我下面的代码。#初始条件x0，y0=0，1 h=0. 1 #步长... ...
我试着比较正弦函数的解析导数和近似导数，用后向、前向、中心、3阶后向差分法，它们之间有巨大的差异。它...
我想在C++中实现第一种和第二种复数贝塞尔函数的描述。现在我正在寻找在我的源代码中引入它们的可能性。...
在这个特殊的问题上，我有一个想象中的城市被划分为若干个方块--基本上是一个覆盖城市的MxN方格。M和N可以比较大，所以我有超过40,000平方的案例......。
如何改进Levenberg-Marquardt的多项式曲线拟合方法？
几周前，我开始在Matlab中从头开始编码Levenberg-Marquardt算法。我对数据的多项式拟合很感兴趣，但我还不能达到准确度的水平......
我正在尝试从本文中实现方程式16：上面的方程式16是粘性波动方程式，应该从大开始并随着时间的流逝逐渐消失，最终接近0。但是...
我正在尝试在Python中创建泰勒级数，但我不知道如何消除此错误：x = Symbol（“ x”）def f（x）：return（（math.e）** x） * sin（x）y = f（x）print（y.diff（x））def Taylor（f，x，m，a）：y = ...
[如何在python中创建函数以使它们返回x作为变量？ （插值）
[如何在Python中创建函数，例如：从数学import *从numpy import *从sympy import *从numpy.linalg import * n = int（input（“结数：”））xsolmed = []对于i in ...
我有一个RKF7（8）积分器，我已经用几个简单的测试功能验证了其输出。但是，当我在感兴趣的运动方程式上使用它时，局部截断误差为...
如何使用scipy.brentq函数，使其等效于pascal代码？
我正在尝试进行一些软件考古。尝试用python中的现代脚本替换旧的pascal文件。在某些时候，pascal脚本正在使用如下所示的例程：FUNCTION ...
我正在实现RKF4（5）集成器，但我无法弄清楚我的代码是否有效，并且我不了解本地截断错误，或者我的代码是否无效。我为代码的大小表示歉意...
供参考，我正在使用此页面。我了解原始的rankrank公式，但是我无法理解为什么稀疏矩阵实现是正确的。下面是它们的代码复制：def ...
为什么我的Runge-Kutta解算器无法随着时间步长的减小而收敛？
出于某种原因，我需要在PyTorch中实现Runge-Kutta4方法（因此，不，我将不使用scipy.odeint）。我尝试过，但在最简单的测试用例上却得到了怪异的结果，用x（0）= 1（...
我正在尝试解决问题。运行此代码时，我实现了想要的功能。（。8372e-03 1.8372e-06 1.8372e-09）truncation.error
信息操作系统：Manjaro Linux开发环境：Spyder4（Python 3.7）库：笨拙的问题您好，我按照以下三种方法编写了一些函数来求解线性方程组：LU ...
我在编程方面很烂，所以请多多包涵。我什至不确定我现在需要的概念是什么，所以我不知道要用谷歌搜索或写这篇文章的标题。我的问题是，...
因此，我试图使用Bulirsch-Stoer自适应方法来解决EDO，但是却收到了溢出警告。如何解决呢？从数学导入sin，cos，pi从numpy导入范围，数组，空，...
我有2个函数：ccexpan-用于计算函数f的插值多项 式与第一种基础的Chebyshew多项式中的N个节点的系数。 csum-为...
关于“计算平均值的最佳数字方法是什么”的建议是，计算滚动平均值，即平均值= a [n] / n +（n-1）/ n *平均值可能在数值上更稳定...