You are here

Answers to self-checks for chapter 8

13 November, 2015 - 11:21

Problem 8.9
First we rewrite the integrand as
\frac{1}{4}(e^{ix}+e^{-ix})(e^{2ix}+2^{-2ix})=\frac{1}{4}(e^{3ix}+e^{-ix}+e^{ix}+e^{-ix})
The indefinite integral is
\frac{1}{12i}(e^{3ix}-e^{-3ix})+\frac{1}{4i}(e^{ix}-e^{-ix})
Evaluating this at 0 gives 0, while at =2 we nd 1/3. The result is 1/3.

Problem 8.8
\begin{align*} sin(a+b)&=(e^{i(a+b)}-e^{-i(a+b)})/2i=(e^{ia}e^{ib}-e^{-ia}e^{-ib})/2i \\&=[(cosa+isina)(cosb+isinb)-(cosa-isina)(cosb-isinb)]/2i \\&=[(cosa+isina)(cosb+isinb)-(cosa-isina)(cosb-isinb)]/2i \\&=cosasinb+isinacosb \end{align*}
By a similar computation, we nd cos(a + b) = cos a cos b - sin a sin b.

Problem 8.10
If z^{3}=1, then we know that |z|= 1, since cubing zcubes its magnitude. Cubing ztriples its argument, so the argument of zmust be a number that, when tripled, is equivalent to an angle of zero. There are three possibilities: 0 ×3 = 0, (2π/3) ×3 = 2π, and (4π/3)×3 = 4π. (Other possibilities, such as (32π/3), are equivalent to one of these.) The solutions are:

z=1, e^{2\pi i/3},e^{4\pi i/3}

Problem 8.11

We can think of this as a polynomial in xor a polynomial in y— their roles are symmetric. Let’s call xthe variable. By the fundamental theorem of algebra, it must be possible to factor it into a product of three linear factors, if the coefficients are allowed to be complex. Each of these factors causes the product to be zero for a certain value of x. But the condition for the expression to be zero is x3 = y3 , which basically means that the ratio of xto ymust be a third root of 1. The problem, then, boils down to finding the three third roots of 1, as in problem 10. Using the result of that problem, we find that there are zeroes when x/y equals 1, e^{2\pi i/3} , and e^{4\pi i/3} . This tells us that the factorization is (x-y)(x-e^{2\pi i/3}y)(x-e^{4\pi i/3}y).

The second part of the problem asks us to factorize as much as possible using real coefficients. Our only hope of doing this is to multiply out the two factors that involve complex coefficients, and see if they produce something real. In fact, we can anticipate that it will work, because the coefficients are complex conjugates of one another, and when a quadratic has two complex roots, they are conjugates. The result is (x-y)(x^{2}+xy+y^{2}).

Problem 8.14

Applying the differential equation to the form suggested gives abx^{b-1}=a^{b+1}x^{b^{2}} . The exponents must be equal on both sides, so b must be a solution of b^{2}-b+1. The solutions areb= (1\pm \sqrt{3}i)/2. For a more detailed discussion of this cute problem, see mathoverflow.net/questions/111066.

Problem 8.15

(a) Let m= 10, 000. We know that integrals of this form can be done, at least in theory, using partial fractions. The ten thousand roots of the polynomial will be ten thousand points evenly spaced around the unit circle in the complex plane. They can be expressed as r^{k}=e^{2\pi k/m} for k= 0 to m1. Since all the roots are unequal, the partial-fraction form of the integrand contains only terms of the form A_{k}/(x-r_{k}). Integrating, we would get a sum of ten thousand terms of the form A_{k}In(x-r_{k}) .

(b) I tried inputting the integral into three different pieces of symbolic math software: the open-source packages Yacas and Maxima, and the web-based interface to Wolfram’s proprietary Mathematica software at integrals.com. Maxima gave a partially integrated result after a couple of minutes of computation. Yacas crashed. Mathematica’s web interface timed out and suggested buying a stand-alone copy of Mathematica. All three programs probably embarked on the computation of the Ak by attempting to solve 10,000 equations in the 10,000 unknowns Ak, and then ran out of resources (either memory or CPU time).

(c) The expressions look nicer if we let w=e^{2\pi /m}, so that r_{k}=w^{k}. The residue method gives

\frac{1}{x^{m}-1}=\sum \frac{1}{(x-w^{k}mw^{k(m-1)})}

Integration gives

\int \frac{dx}{x^{m}-1}=\sum \frac{1}{mw^{k(m-1)}}In(x-w^{k})

(Thanks to math.stackexchage.com user zulon for suggesting the residue mathod, and to Robert Israel for pointing out that for|x|<1 this can also be expressed as a hypergeometric function:

(-x)_{2}F_{1}(\frac{1}{m},1;1+\frac{1}{m};x^{m})