集成R中给“非有限功能值”

integral <- function(x) {2.393794315*((3320)*(x/2.24581)^5+(-1613880/11)*(x/2.24581)^4+(171163181/66)*(x/2.24581)^3+(-7563546913/330)*(x/2.24581)^2+(835541173981/8250)*(x/2.24581)+(-2953570085669/16500))*(((483793.161846485)*x^8+(-76823340.9717028)*x^7+(5337025908.822)*x^6+(-211866341077.587)*x^5+(5256530719898.47)*x^4+(-83466263852549.1)*x^3+(828318375700455)*x^2+(-4697211251008830)*x+(11653475160809900))^0.5)} 

integrate(integral, lower=19.538547, upper=20.3245805)

这给了我集成R中给“非有限功能值”

Error in integrate(integral, lower = 19.538547, upper = 20.3245805) : non-finite function value

我不知道该怎么做。如果有人可以将积分输入到Maple这样的软件中以查看它是否有效 - 或者给出如何解决这个错误的答案也会很棒：P

在此先感谢！

编辑：

，我试图围绕x轴的功能是：

Y = 2.393794315 *（（3320）（X/2.24581）^ 5 +（ - 1613880/11）（x/2.24581）^ 4 +（171163181/66）（x/2.24581）^ 3+（ - 7563546913/330）（x/2.24581）^ 2 +（835541173981/8250）*（x/2.24581 ）+（ - 2953570085669/16500））

当以极限绘图{19.538547 < x < 20.3245805}：

如由@StephenWade提供（在评论），以找到表面积的链接，我需要（F'（X））^ 2 + 1。使用钨，α，我计算这个为：

（f'（x））^ 2 + 1 =（483793.161846485）* x^8 +（ - 76823340.9717028）* x^7 +（5337025908.822）* x^6 +（ - 211866341077.587）* x^5 +（5256530719898.47）* X^4 +（ - 83466263852549.1）* X^3 +（828318375700455）* X^2 +（ - 4697211251008830）* X +（11653475160809900）

绘制此，没有限制，给出了：

然而，R上计算这一点，我得到一个错误：

我的输入成R是在下列：

2.393794315 *（（3320）（X/2.24581 （x/2.24581）^ 3 +（ - 7563546913/330）（x/2.24581）^ 2 +（835541173981）^ 5 +（ - 1613880/11）（x/2.24581）^ 4 +（171163181/66）/（8250）（x/2.24581）+（ - 2953570085669/16500））（（（483793.1 61846485）* X^8 +（ - 76823340.9717028）* X^7 +（5337025908.822）* X^6 +（ - 211866341077.587）* X^5 +（5256530719898.47）* X^4 +（ - 83466263852549.1）* X^3 + （828318375700455）* x^2 +（ - 4697211251008830）* x +（11653475160809900））^ 0.5）

我只是需要一个合理的表面积估计...任何想法？

来源

2016-10-22 StopReadingThisUsername

'integral（20）'给出ǸaN。 –

Maple无法整合它。我尝试过Derive并厌倦了等待它终止，但是当我请求近似积分时，它给出了结果会有“可疑准确性”的警告。你需要重新调整。也许改变变量，例如它是在更大范围内具有较小系数的函数的一个整体。您可能需要咨询数值分析专家，因为天真的解决方案（即使它们返回一个值）可能会遇到派生警告的可疑准确性。 –

@JohnColeman感谢您对Maple的尝试。我是一名高中学生，我不知道如何继续。我试图围绕X轴旋转五阶多项式以找到表面积......任何合理的积分估计值都可以工作;你能帮我吗？ – StopReadingThisUsername

我会用polynom包来处理它，它可以为我们计算导数，然后应用Wolfram的公式。

revol_coef <- 2.393794315 * c(-2953570085669/16500, 
           835541173981/(8250 * 2.24581), 
           -7563546913/ (330 * 2.24581^2), 
           171163181/(66 * 2.24581^3), 
           -1613880/(11 * 2.24581^4), 
           3320/2.24581^5) 
y <- polynomial(revol_coef) 
y_d <- deriv(y) 

f <- function(x) { 
    2 * pi * predict(y, x) * sqrt(1 + predict(y_d, x)^2) 
} 


integrate(f, lower = 19.538547, upper = 20.3245805)

我得到的输出是

-45.71118 with absolute error < 0.0017

答案是否定的，由于表面在x-y平面的右下象限被指定（按你的图）。要解决这个问题，请提供所提供数字的负数以获得表面积。

从Wolfram Alpha翻译公式，操作它们，然后将它们放入R可能很容易导致出错。

在这种情况下，我建议尽可能在R中做很多工作，因为可以随时获得所需的计算和简化的工具/函数。

来源

2016-10-23 01:20:20 stephematician

对此我感激不尽！祝你有美好的一天！ :) – StopReadingThisUsername

确保您了解每一步，因为每个人（包括我自己）都可以犯错误。我没有试图验证我的答案。 – stephematician

集成R中给“非有限功能值”

回答

相关问题