问题标签 [simpsons-rule]

我目前正在尝试制作一个程序，它将使用 matplotlib 绘制一个函数，绘制它，对两个变量之间的曲线下的区域进行着色，并使用辛普森的 3/8 规则来计算阴影区域。但是，当尝试打印我分配给积分最终值的变量时，它会打印一个列表。

首先，这是我的代码的基础：

这个定义定义了我将在这里使用的函数，一个简单的多项式。

这是计算曲线下面积的函数

现在辛普森的规则定义结束了，为了简单起见，我定义了一些变量。

在这一点上，肯定有什么问题发生了，但我将包含其余代码，以防您进一步发现问题。

这里是应该打印积分值的地方：

以下是绘制图表所需的其余代码：

当运行这个程序时，你会得到这个图表，并且在曲线下面积应该是的地方打印了一系列数字

任何帮助在这里表示赞赏。我完全被困住了。另外，对不起，如果这有点长，这是我在论坛上的第一个问题。

我刚开始学习MATLAB。练习的目的是使用 Simpson 的 1/3 规则和 romberg 进行近似/积分。问题是将 x^(1/2) 从 0 积分到 2

当我执行时： simpson(fun,0,2,10)

fun = x^(1/2); 我在第 2 行：或在 simpson 的第 16行收到错误：f = feval(fun,x);

谢谢您的帮助！

这是我的方程式代码：

我的辛普森代码：

我的 romberg 代码：

我编写了以下代码来使用辛普森规则来集成任意函数，在我的例子中是 sin(x)：

但是我在 sum+=f(x1+i*delta_x) 行上收到错误“float object not callable”。有谁知道可能出了什么问题？

谢谢 ：）

我为 Simpson 的规则集成编写了以下代码来近似 sin。方程式在此附件中。我为偶数和奇数项编写了单独的循环，因为它们在附件中分组显示。

但是当我打印这些值时，它给了我非常小的负数。

谁能看到我的代码有什么问题？

我以向量化的方式使用 scipy simpsons 规则来提高性能。（我在时间序列的滚动窗口上计算积分。）

由于某种原因，在有 4GB 可用内存的 64 位机器上对 500MB 浮点数组使用 simps 时出现内存错误。

这是一些示例代码：

我可以看到我的机器的虚拟内存一直在上升，直到使用了 2GB，然后我得到了一个 MemoryError，即使我的空闲内存没有被使用。

所以我很困惑：即使仍有大约 2GB 的未使用内存可用，为什么我会收到 MemoryError？这可以解决吗？

我正在使用 Windows Server 2012 R2。

编辑：为了说明内存使用与输入大小的或多或少的线性缩放，我进行了以下小实验：

这使得 Windows 资源监视器以及其他工具（如 mprof）中的内存消耗增加：

所以说我有来自用户输入的积分的下限值和上限值。我先询问下限，然后检查其有效性。然后为了比较我的上限和我的下限的值，我做了一个嵌套函数，这样我也可以要求用户输入上限，检查它的有效性并确保我的上限大于我的下限（因为你知道集成），如下面的代码所示。

现在这一切都很好，直到我真正需要使用值 a 和 b 因为我需要将这些值放入我设置的积分中。它基本上是由辛普森规则构成的一般积分，现在很容易进入。所以我将函数定义为：

好的，所以我的函数基本上允许用户插入任意函数，上限（a）、下限（b）和 N（辛普森规则中的条带数），它吐出 d，这是对任意函数的积分的评估通过辛普森规则实现上述功能。我现在的问题是，当我尝试打印答案时，我可以将变量 a 取出并放入积分中，但我不能取出变量 b，因为它在函数中！例如，如果我现在打印积分值（例如在 sin(x) 和 N = 20 的情况下）

所以我知道 a 和 b 值在它们自己的函数中是局部的。现在对于 a 的值，我可以很容易地做到这一点来获得值 a

因为因为 k 调用 LowLimCheck() ，它返回我可以放入函数中的 a 的值。但是我怎样才能得到嵌套在第一个函数中的 b 值呢？我基本上想用 b 。有办法解决这个问题吗？

对冗长的问题表示歉意，并提前致谢！

Linklist 的 Simpson 1/3 积分方法。

期望输出为 60.00 但得到 0，不知道为什么？ 尝试通过辛普森的 1/3 积分方法对任何非线性方程进行积分。在代码块IDE中试过这个。应用 Simpson 的 1/3 积分规则的正确逻辑，仍然使 Integrand 值始终为零，不知道我在这段代码中哪里做错了。

我想使用 c 程序绘制与给定曲线相对应的 Simpson 1/3 规则中的积分面积？

这是我的代码：

如何使用 C 编程绘制区域？我可以使用 Python 做到这一点，但想用 c 做到这一点。

问题的解决方案是辛普森 1/3 规则的 1.732400451459101。相反，程序给我的解决方案是 1.73239801

谁能帮我吗？提前致谢。

我想 1. 将辛普森规则表达为 python 中集成的通用函数，2. 用它来计算和绘制函数的傅里叶级数系数​​。

我已经为辛普森规则窃取并修改了这段代码，它似乎可以很好地集成简单的功能，例如, 或

给定 period ，傅立叶级数系数​​计算为：

其中 k = 1,2,3,...

我很难弄清楚如何表达。我知道，因为这个函数很奇怪，但我希望能够为其他函数计算它。

到目前为止，这是我的尝试：

这根本没有给出 0 的正确答案，因为它随着 k 的增加而增加，而且我确信我正在用 for 循环的隧道视觉接近它。特别是我的困难在于说明函数，以便将 k 读取为 k = 1,2,3,...

不幸的是，我遇到的问题没有指定要绘制的系数是什么，但我假设它是针对 k 的。