通过数值逼近（与符号计算相反）解决数学问题的算法。
f = @（t，y）[y（2）; cos（t）-y（2）-4 * y（1）] a = 0 b = 1 N = 5000 h = 1 y（1）= 0 y（2）= 0.1 t（1）= a对于n = 1：1：N t（n + 1）= t（n）+ hy（n + 1）= y（n）+ h * f（t （n），y（n））；我正在尝试结束图（t，y）...
我正在使用scipy.optimize.newton_krylov解决以下问题。但是，它相当慢（这完成了数千个点）。即使我最初的猜测非常接近解决方案，也可以...
我遇到了涉及包含未知变量的积分的数值估计的优化问题。数值估计积分很简单，只需使用...
我正在尝试为以下方程式实现Laguerre的方法：1 /（99 * x + 1）= 2。 （一般形式更复杂-n次多项式1 /（a * x + 1）+ ... + 1 /（z * x + 1）= res的多项式）其中a，b，... z> = 0和0 ]
计算sqrt（a）的一种方法，a> 0是X（n + 1）=（a +（X（n）* X（n-1））/（X（n）+ X（n -1）），n = 1，2，…，且X0 = 1且X1 = a（即，已知lim n-> Xn的无穷大= sqrt（a）...
以下是曲线f（r）的图，其中r是径向坐标，并针对不同的参数值进行了绘制，如下所示：但是，我不知道曲线的功能形式，我是。 ..
我目前正在尝试使用一些通用的显式Runge-Kutta方法来编写一些Python代码，以解决一阶ODE的任意系统，该方法由值alpha，gamma（两个向量均为...]定义）
程序po隐式无整数:: i，j，nc，nd，ok，iter REAL ::阿尔法，赌注，chi REAL，DIMENSION（:)，可分配:: uexact，...
我正在使用Euler方法在物理引擎中对粒子进行数值积分。浮逆质量; Vector3的位置，速度； Vector3 sumForces，重力，inputForce;空粒子:: Integrate（float dt）...
我将如何使用Octave数值求解以下简单的微分方程组？注意：根据我的理解，我使用限定词“简单”，因为系统首先是...
如何创建特定的数字模式？最后一个数字从第二个模式部分开始。示例：我想跟随/重复第1,5,4节中的模式1,5,4第一模式...
我有一个带两个已知变量x和y的指数函数。输入x时，我需要找到y的值。但是，我的代码无法解决该问题。我的功能和所有...
使用32位浮点值，如果-开始计算时，如果我不知道要有多少个值（在...中），那么最好的（数字上最准确的）平均值计算方法是什么？
我正在尝试使用最初用C语言编写的Python复制DSP算法。诀窍是，我还需要保留C版本中32位定点变量的相同行为，...
我尝试了以下代码，使用增量搜索技术找到方程式的根。当我将精度设置为5或7位数字时，代码可以正常工作。但是当我将精度设置为6 ...
用scipy Solve_bvp解决边界值问题（扩散反应方程）
我正在努力解决以下二阶边值问题：y''+ 2 / x * y'+ k ** 2.0 * F（y）= 0 y（x = 1）= 1，y'（ x = 0）= 0 F（y）= -y或F（y）= -y * exp（AB *（1-y）/（1 + B（1-y）））我以某种方式无法设置。 。
因此该方法本身可以正常工作：double trap_method（double a，double b，double n）{int i = 1;两倍的面积= 0;双倍h =（b-a）/（n-1）;双x = a + h; while（i <= n-2）{area = ...
我正在尝试实现中点公式y [n + 1/2] = y [n] + h / 2 f（x [n]，y [n]）和y [n + 1] = y [n] + h * f（x [n] + h / 2，y [n + 1/2]），因此可以使用中点法求解ODE。我的函数是函数[...
我正在尝试使用角谱方法在数值上传播给定的（电场）。为此，我遵循“傅立叶光学的原理和应用”（Robert K. Tyson）第3章，...
我正在尝试实现光谱聚类算法，以解决图问题中的社区。我有一个非常大的矩阵来计算其特征向量，矩阵> 1Mx1M。 Numpy和Scipy需要...