通过数值逼近（与符号计算相反）解决数学问题的算法。
因此，我试图使用Bulirsch-Stoer自适应方法来解决EDO，但是却收到了溢出警告。如何解决呢？从数学导入sin，cos，pi从numpy导入范围，数组，空，...
我有2个函数：ccexpan-用于计算函数f的插值多项 式与第一种基础的Chebyshew多项式中的N个节点的系数。 csum-为...
关于“计算平均值的最佳数字方法是什么”的建议是，计算滚动平均值，即平均值= a [n] / n +（n-1）/ n *平均值可能在数值上更稳定...
如何将这个情节变成动画？我有一个情节，但想以某种方式将其变成动画
在最近的几个晚上，我一直在努力如何在每个时间步长或每个x步长后将此处的波形转换为某种动画。如何修改和编写代码，以便...
我已经在python中使用bvp求解器解决了以下bvp问题。 d4y / dx4 = 0.00033 * V /（0.000001-y）^（2），y（0）= y'（0）= y（1）= y'（1）= 0在上面的等式'V'是一个使用...
因此，我的教授要求我创建一个程序，该程序将使用Runge-Kutta方法（二阶）来解决该问题。但是我得到一个IndexError：列表索引超出范围，这可以回到我的第一个...
您如何在Python中实现基于Householder的QR分解？
我目前正在尝试对矩形矩阵实施基于Householder的QR分解，如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf（第3、4、5页）所述。显然...
如何在Python中为矩形矩阵实现基于Householder的QR分解？
我目前正在尝试对矩形矩阵实施基于Householder的QR分解，如http://eprints.ma.man.ac.uk/1192/1/qrupdating_12nov08.pdf（第3、4、5页）所述。显然...
PDE解决：'NoneType'没有属性'toarray'错误
我正在学习如何在Python中采用有限差分方法，以泊松方程为例。但是，当我运行时，一切都很好，直到A = ASparse.toarray（）逆向代码行...
单面Jacobi SVD和Divide＆Conquer SVD稳定性和误差
我正在MATLAB中研究SVD，尤其是Jacobi SVD和Divide＆Conquer SVD算法。在这些方法的稳定性和错误分析上，我找不到任何东西。有人可以显示...
我正在尝试使用第4阶隐式Runge-Kutta方法（IRK）编写一个能够解决ODES系统的函数，但是我无法正确定义循环。在这里，我们通过...
等式为：d ^ 2 r / dt ^ 2 = -c / m（dr / dt）+ g，其中r是弹丸的位置，c是阻力系数，m是弹丸的质量，g是重力引起的加速度。假设...
我正在尝试求解以下耦合方程：x = 1; y-0.5 * y-0.7 * v = 0; w-0.7 * x-0.5 * x = 0; v =1。（我知道方程= 1似乎是不必要的，但以后我需要它们...
我如何在PYTHON中求解带有空气阻力的弹丸运动的二阶微分方程？
等式为：d ^ 2 r / dt ^ 2 = -c / m（dr / dt）+ g，其中r是弹丸的位置，c是阻力系数，m是弹丸的质量，g是重力引起的加速度。假设...
我正在使用斜率方法计算函数的导数。我需要放大图表以查看实际情况，但是数据点太低而无法给出平滑的曲线。所以当我...
我正在解决我的宏观经济学课上的问题。考虑以下方程式：在此，k是固定的，而c（k）是通过Matlab中的``interp1''函数定义的。这是我的代码：beta = 0.98; ...
Brent方法的停止标准是，如果abs（m）<= tol或fb == 0.0，则//找到根（间隔足够小）发现：= true;但是，如果abs（m）达到上述容差以下，但...
嗨，我被要求使用MATLAB中的fsolve命令来求解SIR模型，而Euler 3向后指向。我对如何继续感到很困惑，请帮忙。这是我到目前为止所拥有的。我创建了一个函数...
n阶ODE的四阶Runge-Kutta积分器的索引列表超出范围误差
问题参数：您要构建的函数需要使用自变量的当前值（x）因变量的当前值（y）
我的泰勒级数图的代码cos（2x）代替了Matlab中的cos（x）
我正在编写一个计算任何函数的泰勒级数的函数。 syms x y = cos（x）; y0 = 0; a = 0;对于i = 0：25 diff（y，i）; %%给出导数公式...