16得票3回答
R中的双重积分

我想知道如何在R中编写代码来进行二重积分。我已经参考了两个类似的问题。 快速计算R中的双重积分 带有额外参数的R中的双重积分 但是我仍然不清楚如何从这些答案中得到我的问题的答案。我的问题如下。 我想在R中编写代码计算以下内容。 根据我的手算和Wolfram Alpha计算,结果为...

11得票3回答
如何在方程输出中保留分数

我一直在使用Python来计算数学方程。例如:from sympy import Symbol, Derivative, Integral x = Symbol('x') d = Symbol('d') Integral(8*x**(6/5)-7*x**(3/2),x).doit() 这将导致...

8得票1回答
SymPy 定积分

我希望将以下方程式整合: f(x) = 积分(E^(-i * omega * t)), 从 -a 到 a. 我写了以下代码: from sympy import * from sympy.abc import a, omega, t init_printing(use_unicod...

19得票2回答
如何增加`scipy.integrate.dblquad`函数中的子区间数量?

我正在使用scipy.integrate.dblquad,但是我收到了这个错误:UserWarning: The maximum number of subdivisions (50) has been achieved. If increasing the limit yields no i...

19得票2回答
在R中快速计算二重积分

我正在寻找一个比<scipy.integrate.nquad>更快的二重积分解决方案。integrate(function(y) { sapply(y, function(y) { integrate(function(x) myfun(x,y), llim, u...

15得票4回答
Python函数:查找数值体积积分?

目标 我想计算数值标量场的三维体积积分。 代码 在本文中,我将使用一个可以精确计算积分的示例函数。因此,我选择了以下函数: 在Python中,我定义了该函数和一组三维点,然后在这些点上生成离散值: import numpy as np # Make data. def fun...

12得票2回答
用积分函数拟合数据

使用Python中的curve_fit函数从scipy.optimize库拟合数据时,首先需要定义拟合函数(例如二次多项式)如下: def f(x, a, b): return a*x**2+b*x 然后进行拟合 popt, pcov = curve_fit(f,x,y) 但...

10得票2回答
如何使用SymPy加速长函数的符号导数计算?

我正在用Python编写程序,使用自由ICI方法(目前是SICI方法...但最终将变成自由ICI)来解决薛定谔方程。如果您对此不熟悉,那是因为这个主题几乎没有信息可供参考,并且完全没有可用的样本代码。 这个过程涉及到迭代地得出偏微分方程的解。在这个过程中,需要执行许多符号求导操作。问题在于,...

7得票5回答
编写一个Python函数,用于整合高斯分布的最佳方法是什么?

在尝试使用scipy的quad方法来积分高斯函数(假设有一个名为gauss的高斯函数),我遇到了一些问题,无法将所需参数传递给gauss并使quad在正确的变量上进行积分。有人有使用带有多维函数的quad的好例子吗? 但是,这引发了一个更大的问题,即通常最好的积分高斯函数的方法是什么。我在s...

11得票4回答
在编程中实现数学方程式遇到问题

我正在为学校实现一个算法,但在理解如何在编程中表示定积分方面遇到了问题。例如,我知道求和公式可以实现为以下示例: 假设y=f(x)if(x==0){ y=x+1; }else{ for(int i = 0; i < n; i++){ y = y + (x...