通过数值逼近(与符号计算相反)解决数学问题的算法。
我正在尝试为有限差分法设计算法,但我有点困惑。所讨论的ODE是y'' - 5y'+ 10y = 10x,y(0)= 0且y(1)= 100。所以我需要一种方法以某种方式获得...
我是R的新手,想学习R如何用数字方法求解方程式的基础知识。我的问题是:当没有解决方案时,fsolve会返回错误。为什么?怎么才能让它回归南...
我想在Matlab中计算匿名函数的二阶导数。我已经知道了这个(数值微分)的一些公式,但它们似乎不起作用。我能算出......
几年前,在20世纪90年代早期,我构建了基于定点算法优化计算的图形包,以及用于sqrt和log的cos,sin和缩放方程的预计算表...
我试图通过在eigen3对象上使用系数方法来消除嵌套for循环。为了达到这个目的,我必须概括一个已经存在的功能,这样我就可以...
特别是,如果在计算机上实现欧拉方法,那么在舍入误差之前可以使用的最小步长是什么导致欧拉近似变得完全不可靠?一世 ...
对于数值方法类,我需要编写一个程序来评估Simpson复合规则的定积分。我已经走到这一步了(见下文),但我的回答并不正确。我正在测试......
我想用这个问题来改进我对计算机工作原理的一般理解,因为我可能永远都没有机会深入而深刻地学习。提前抱歉,如果......
通过“相似矢量”,我定义了一个矢量,该矢量在给定的一个位置上与-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)+ ...的系列,其条件是总和应该停止当。。。的时候 ...