Using an expanded power series to estimate

1. Nov 4, 2009

NastyAccident

1. The problem statement, all variables and given/known data

(a) Expand f(x) as a power series

$$f(x)=\frac{7}_{\sqrt[4]{1+\frac{x}_{14}}}$$

Which I converted to....

$$7 - \frac{1}{8}*x + \sum^{\infty}_{n=2}7*(-1)^{n}*\frac{1*5*9*...*(4n-3)}{4^{n}*n!}*(\frac{x}{14})^{n}$$

(b) Use part (a) to estimate 7 / (1.1)^(1/4) correct to three decimal places.

2. Relevant equations
Knowledge of Taylor & Maclaurin Series

3. The attempt at a solution

$$7 - \frac{1}{8}*(1.1) + \sum^{\infty}_{n=2}7*(-1)^{n}*\frac{1*5*9*...*(4n-3)}{4^{n}*n!}*(\frac{1.1}{14})^{n}$$
=0.981268

Though, I think this is more on the right track:

1.1=1+x/14
.1=x/14
1.4=x

$$7 - \frac{1}{8}*(1.4) + \sum^{\infty}_{n=2}7*(-1)^{n}*\frac{1*5*9*...*(4n-3)}{4^{n}*n!}*(\frac{1.4}{14})^{n}$$
=0.976454

As always, any guidance will be appreciated and thanked =)

Sincerely,

NastyAccident

2. Nov 4, 2009

clamtrox

I am fairly sure that part a should be more difficult than part b.

3. Nov 4, 2009

lanedance

you have expanded about x=0, so the 2nd part where you subtitute x=1.4 into your formula is the correct way to get
$$\frac{7}{(1+x/14)^{1/4}} = \frac{7}{(1.1)^{1/4}}$$

however i'm not sure how you get your value, you series look reasonable though i haven't checked it 100%...

checking the value with a calculator I get ~6.8. You need to keep adding terms untill the 3rd decimal place is stationary

consider the first 2 terms
7 - (1.4)/8 = 6.825
each term will be close to an order of magnitude less than the previous (as it has (x/14)^n) which should give you a hint as too how many terms are required...

4. Nov 4, 2009

NastyAccident

Hehe, late at night part a. is ten times easier than part b.

7-[7*(1)/(4*1!*14)](1.4) = 6.825

7-[7*(1)/(4*1!*14)](1.4)+[7*(1*5)/(4^2*2!*14^2)](1.4) = 6.83281

7-[7*(1)/(4*1!*14)](1.4)+[7*(1*5)/(4^2*2!*14^2)](1.4) - [7*(1*5*9)/(4^3*3!*14^3)](1.4) = 6.83239

7-(1/8)(1.4) + 7(5/(4^2(2!)*14^2))(1.4) - 7((5*9)/(4^3(3!)*14^3))(1.4) + 7((5*9*13)/(4^4(4!)*14^4))(1.4) - 7((5*9*13*17)/(4^5(5!)*14^5))(1.4)

= 6.83242

So,

6.832

However, that is still being picked up as being wrong.... For some odd reason. I know that the series is correct as it's already showing as correct.

Sincerely,

NastyAccident

Last edited: Nov 4, 2009
5. Nov 4, 2009

TheoMcCloskey

NastyAccident - you're almost there! However, you're not computing the power term correctly. The full term (1.4/14) should be the power term.

Let's review. You've expanded the denominator correctly - Let's put it into a more general form. The expansion of the denominator has the following form:

$$g(1+\xi) = g(1) + \sum_{k=1}^{\infty} \beta_k \, \xi^k$$
where
$$\beta_k = (-1)^k \, \frac{\Pi_{m=1}^{k} (4m-3)}{4^k \cdot k!}$$

and
$$\xi = x /14$$

Of course, the full function approximation to $f(x)$ would be $f(x) = 7 * g(1+\xi)$.

Now, if we compute the first few coefficients:

$$\beta_1 = -1/4 = -0.25$$

$$\beta_2 = \beta_1 \cdot (-5/4) \cdot (1/2) = 5/32 = 0.15625$$

$$\beta_3 = \beta_2 \cdot (-9/4) \cdot (1/3) = -45/384 = -0.1171875$$
Now, for the desired approximatiopn, you've also correctly realized that, here, the value of $\xi=1.4/14 = 0.100$.

It is here that lanedance was trying to give you a hint to look at the value of the various terms, considering the powers of $\xi$ and coefficients, to determine the numbers of terms that we might need for 3 digit accuracy. Since $\xi=0.1$, each succesive power will decrease the significance of the term by approximately one decade.

For example, consider the magnitude of $\beta_2 \xi^2$ term: the contribution to the overall sum for f(x) from this term will be approximately
$$|7 \cdot \beta_2 * \xi^2 |= 7 \cdot |\beta_2| \cdot (0.1)^2 = 7 \cdot (1.56 \, \cdot 10^{-1}) \cdot 10^{-2} = 7 \cdot (1.56 \, \cdot 10^{-3}) = 1.09 \, \cdot 10^{-2}$$.​

(Note the multiplication by 7). Thus we see that terms beyond this need to be taken into account. What about the next term?

$$|7 \cdot \beta_3* \xi^3|= 7 \cdot |\beta_3| \cdot (0.1)^3= 7 \cdot (1.17 \, \cdot 10^{-1}) \cdot 10^{-3} = 7 \cdot (1.17 \, \cdot 10^{-4}) = 8.20 \, \cdot 10^{-4}$$.​
We see that the k=3 term will be less than 0.001 (least significant digit for three place accuracy) and will not effect the third decimal place on truncation (it will, however, effect the result upon rounding).

This is sufficient for the result to be correct to three decimal places. If we wanted to be assured that the rounded result would be accurate to three decimal places, than we would have more work to do (see below).

Hence, to be correct to three decimal places, two terms (k =1, 2) plus the "zero" term (g(1)) is needed in the expansion

Computing this requires the following:
$$y = 7 \cdot (1+\beta_1 \cdot \xi + \beta_2 \cdot \xi^2 )$$

which yields
$$y = 7 \cdot (1 - \frac{1}{4} \xi + \frac{5}{32} \xi^2) = 7 \cdot (1 - 0.25 \cdot 10^{-1} + 0.15625 \cdot 10^{-2})$$

which computes as
$$y = 7 \cdot (0.976563)$$
$$y = 6.835938$$

The actual computed value is
$$\frac{7}{\sqrt[4]{1.1}} = 7 \cdot (0.976454) = 6.835179$$

Notice that the approximation is correct to three decimal places, however, the rounded approximation to three decimals places would be in error in the last digit.

As an aside, consider the rounded to three decimal places condition (not part of OP problem).
From prior analysis with three $\beta$ terms, we would expect (hope!) that further terms are not needed since $\xi$ will influence the next term by a factor of ten, but what about the coefficient? A quick calculation shows that the next coefficient is $\beta_4 = \beta_3 \cdot (-13/4)*(1/4)$ and its magnitude is less than $\beta_3$.

Hence, a total of three terms (k = 1,2,3) plus the "zero" term (g(1) = 1) should be used in the expansion for proper approximation to three digits.

Computing this requires the following:
$$y = 7 \cdot (1+\beta_1 \cdot \xi + \beta_2 \cdot \xi^2 + \beta_3 \cdot \xi^3)$$

which yields
$$y = 7 \cdot (1 - \frac{1}{4} \xi + \frac{5}{32} \xi^2 - \frac{45}{384} \xi^3) = 7 \cdot (1 - 0.25 \cdot 10^{-1} + 0.15625 \cdot 10^{-2} - 0.1171875 \cdot 10^{-3})$$

which computes as
$$y = 7 \cdot (0.976445)$$
$$y = 6.835117$$

The actual computed value is
$$\frac{7}{\sqrt[4]{1.1}} = 7 \cdot (0.976454) = 6.835179$$

The above approximation is now accurate to three decimal places.

6. Nov 4, 2009

NastyAccident

I, literally, just realized that I wasn't using the nth power of x, but was just using x. Such a silly and idiotic error on my part through off the calculation big time.

The new calculation:

7-(1/8)(1.4) + 7(5/(4^2(2!)*14^2))(1.4)^2 - 7((5*9)/(4^3(3!)*14^3))(1.4)^3 + 7((5*9*13)/(4^4(4!)*14^4))(1.4)^4 - 7((5*9*13*17)/(4^5(5!)*14^5))(1.4)^5 = 6.83518

Thank you folks! I truly do appreciate it!

Sincerely,

NastyAccident.