通过数值逼近(与符号计算相反)解决数学问题的算法。
通过“相似矢量”,我定义了一个矢量,该矢量在给定的一个位置上与-1或1不同。但是,如果给定的元素为0,则只能差1。例子:...
对于方程组,fsolve,brentq和root的使用和精度有何不同?
我问过这个问题,对于任何方程组来说都很好吗?从中得到了满意的答案。我在那里提出的系统x = A * exp(x + y)y = 4 * exp(x + y),只是...
我的想法是计算以下向量场和曲线的线积分:这是我尝试过的代码:从sympy import *导入numpy为np来自sympy import Curve,line_integrate ...
新的一天新问题......我按照Eric所说的方式做到了,我创建了一个加速方法!它看起来像这样:静态Vector2 AccelerationOfTheMoon(Vector2位置){double Mm = 7 ....
我有以下积分(更多细节:https://math.stackexchange.com/questions/3193669/how-to-evaluate-the-line-integral-checking-stokes-theorem)C_3可以用三角函数进行评估...
如何实现自适应步长Runge-Kutta Cash-Karp?
尝试实现自适应步长Runge-Kutta Cash-Karp但失败了这个错误:home / anaconda / lib / python3.6 / site-packages / ipykernel_launcher.py:15:RuntimeWarning:除以零......
我正在尝试使用4阶Runge Kutta方法求解Lorenz系统,其中dx / dt = a *(yx)dy / dt = x(bz)-y dx / dt = x * yc * z因为这个系统没有不依赖于时间,这可能是......
我在C中有一组微分方程(从一个采用xml输入的工具创建),形式如下:#ifdef ODEs dx [0] = Function1(p [67],p [64],p [66], p [65],p [23],x_c [0],x_c [3],p [49]); ...
我试图找出立方体的每个面的表面积和相应的向外单位法线。这个操作是在有限元网格上完成的,所以我改变了每个表面......
Matlab - 为什么有这么多求解线性方程组的方法? [关闭]
这可能听起来像一个简单的问题,但我想知道为什么有这么多不同的内置Matlab方法来求解线性方程组。我已经完成了以下所有定义:......
当使用数值方法(例如Runge-Kutta)时,计算机上浮点数的有限精度会影响解(Brouwer定律)。在本文中,它建议作为一种模拟精确的补救措施......
标题几乎是不言自明的。我试图近似一个可写为f(x)=x¹/(1x2)-x³/(3x4)+ x ^(5)/(5x6)+ ...的系列,其条件是总和应该停止当。。。的时候 ...
设A是一个nxn稀疏矩阵,由一系列形式为(i,j,a)的m元组表示,其中索引为i,j(介于0和n-1之间),a为底层的值a领域F.什么算法......
我使用Python 3进行与数值分析相关的任务。我必须绘制一些源自正弦函数的点。此外,我需要对这些点进行三次插值(三次样条)。 ...
我有一个问题,我的二进制搜索算法找到2的平方根似乎是在一个无限循环并永远运行:num = 2 low = 1 high = num i = 0 while((i ** 2)!= 2 ):#was while(low&...
我正在积分偏微分方程,其中我在$ x $中有一个四阶偏导数,时间上的数值积分给了我荒谬的误差。问题的原因,我......
我在javascript中有大量的数值y。我想通过将它们舍入到最接近的x的倍数来对它们进行分组,并将结果转换为字符串。我怎么绕过烦人的...
如何绘制庞加莱部分? (Duffing Oscillator)
我写了一个程序,成功地显示了Duffing方程的简单极限环。但是,我现在需要为此案例绘制庞加莱部分。我需要拍摄快照......
以数字方式解决问题意味着什么?为什么用数字解决问题很重要?
我想评估N维中蒙特卡罗方法所需的某个总和。在N维我有N个变量与M个数据点,这意味着我基本上有M个N-D向量代表M ...