集成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中給「非有限功能值」

回答

相關問題