Can Anyone Help Identify This Series and Explain MATLAB's Incorrect Answer?

Click For Summary

Discussion Overview

The discussion revolves around identifying a specific infinite series represented by an equation and understanding why MATLAB's output for the roots of a polynomial derived from this series may not be accurate. The scope includes theoretical exploration of series, polynomial roots, and MATLAB's computational methods.

Discussion Character

  • Exploratory
  • Technical explanation
  • Debate/contested

Main Points Raised

  • One participant presents an equation as a series and seeks help in identifying it as a function or product of functions.
  • Another participant argues that MATLAB's output is incorrect, suggesting that the series is an infinite power sum and resembles the power series for e^-x, which does not equal zero.
  • A different participant speculates that the factorial in the series causes the terms to diminish rapidly, implying that a sufficient number of terms could yield an accurate approximation.
  • Another response suggests that the issue with MATLAB's output may stem from not formulating the problem correctly and mentions the potential ill-conditioning of the companion matrix for truncated series.
  • One participant recommends graphing the series to better understand its behavior, especially for larger values of x.

Areas of Agreement / Disagreement

Participants express differing views on the validity of MATLAB's output and the nature of the series. There is no consensus on the correct identification of the series or the accuracy of the computational results.

Contextual Notes

Participants note limitations related to the formulation of the problem in MATLAB and the behavior of the series for different values of x. The discussion highlights the challenges in working with infinite series and polynomial approximations.

sparklingway
Messages
2
Reaction score
0
I have to find the "second smallest root" of the following equation :

1-x+(x^2)/(2!)^2-(x^3)/(3!)^2+(x^4)/(4!)^2+...=0

Matlab returns quite a satisfactory answer. >> p=[1/518400 -1/14400 1/1576 -1/36 1/4 -1 1]

p =

0.0000 -0.0001 0.0006 -0.0278 0.2500 -1.0000 1.0000

>> roots(p)

ans =

35.5690
-4.6796 +18.5352i
-4.6796 -18.5352i
4.1776 + 3.2154i
4.1776 - 3.2154i
1.4350

But I have been asked to identify this series as well, which I am unable to do. Can anybody help me identify this series as a function or a product of functions? Thanking anybody who answers before hand
 
Physics news on Phys.org
But Matlab does NOT give a solution to that equation! What you have done is enter the first seven terms of the series as if it were a 6th degree polynomial. If I am understanding your question correctly, that is an infinite power sum. In fact, it looks to me like the power series for e-x and that is never 0.
 
Why would MATLAB not give a correct answer. From what I can infer, as the nth terms rises the factorial rises dramatically, therefore the fraction becomes very minuscule in value. Ten terms would be sufficient to give an answer correct to the third or fourth decimal place.

The power series of e^x is 1+x+(x^2)/(2!)+(x^3)/(3!)+...

The problem with this series is the square of the nth factorial is involved.
 
sparklingway said:
Why would MATLAB not give a correct answer.
Two reasons: (1) you're not feeding it the right problem, and (2) the companion matrix is quite ill-formed for a truncated series of order 6 or more.

Try graphing the series. You will need to be a bit creative in how you perform the summation if you want to see behavior for anything but small values of x.
 

Similar threads

Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
1
Views
2K
Replies
8
Views
2K
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K