3
$\begingroup$

What is $x$ in closed form if $2x-\sin2x=\pi/2$, $x$ in the first quadrant?

  • 0
    Yes. But it seems from other posters' answers, there is no known one due to $cosy=y$.2010-11-15

5 Answers 5

8

The solution is given by $\displaystyle x = \pi/4 + D/2$

where $\displaystyle D$ is the root of $\cos y = y$

The root of $\displaystyle \cos y = y$ is nowadays known as the Dottie Number and apparently has no known "closed form" solution. If you consider this number to be part of your constants, then the above can be considered a closed form solution.


For a proof:

If $\displaystyle y = \sin(2x)$

then we have that

$\displaystyle 2x = \pi/2 + y$

$\displaystyle y = \sin 2x = \sin (\pi/2 + y) = \cos y$.

The root of $\displaystyle y = \cos y$ is $\displaystyle y = 0.739085\dots$

Notice that $\displaystyle \pi/2 > x \gt \pi/4$ (as $\displaystyle f(x) = 2x - \sin 2x$ is increasing in $\displaystyle [0,\pi/2]$), so if $\displaystyle x = \pi/4 + z$ then

$\displaystyle \sin(2x) = \sin(\pi/2 + 2z) = \cos 2z = 0.739085\dots$

And thus $\displaystyle z = \dfrac{0.739085\dots}{2}$.

Thus $\displaystyle x \sim \pi/4 + \dfrac{0.739085}{2} \sim 1.154940730005\dots$

See Also: A003957.

2

Repeating the same thing as the other answers, but with a moderately more elegant one:

Let $y = 2x - \frac{\pi}{2}$, then substituting:

$y = sin(y+\frac{\pi}{2})$

$y = cos(y)$

$y = D$

then

$x = \frac{y+\frac{\pi}{2}}{2} = \frac{D}{2} + \frac{\pi}{4}$

2

An analytical form of x can be obtained solving Kepler equation:

$M= E-\epsilon \sin(E)$

with eccentricity=1 and mean anomaly = $\pi/2$ by means of Kapteyn series:

$2x = \frac{\pi}{2}+\sum_{n=1} \frac{2J_n(n)}{n} \sin(\pi n/2)$

where $J_n()$ are the Bessel functions. Simplifying:

$2x = \frac{\pi}{2}+\sum_{n=0} \left( \frac{2J_{4n+3}(4n+1)}{4n+1} - \frac{2J_{4n+3}(4n+3)}{4n+3}\right)$

$x = \frac{\pi}{4}+\sum_{n=0} \left( \frac{J_{4n+1}(4n+1)}{4n+1} - \frac{J_{4n+3}(4n+3)}{4n+3}\right)$

Such series can be numerically evaluated, but it converges slowly and n=10000 terms are required to obtain:

$x = 1.154940317134$

with

$2x-\sin(2x)-\pi/2=-1.38017659479e-006$

In order to improve the convergence, we can employ an acceleration series technique as Levin's acceleration. (See http://en.wikipedia.org/wiki/Series_acceleration)

With only 10 (ten!) terms we obtain:

$x=1.1549406884223$

A simple c++ code, based on gsl library is the following:

#include  #include  #include  #include "gsl_sf.h" #include "gsl_sum.h" using namespace std;  #include   int main(int argc, char* argv[]) {     double PIH = atan(1.)*2;     cout<
  • 0
    I fixed the math. As you can see the code was consistent with the correct formula... I will fix the math also at the Dottie thread2015-03-04
1

To get a graphical impression of the solution have a look here:

http://www.wolframalpha.com/input/?i=solve+2+x+-+sin+(2+x)+%3D+Pi/2+for+x

Edit: That WA doesn't give a closed form is at least a hint that there might be none (very probable)

  • 2
    @Moron & @vonjd: You're both correct. As I said somewhere here previously, "closed-form" is just shorthand for "popular enough to be given a name and notation."2010-11-15
0

As others have noted, $ x = \pi/4 + D/2. $

Based on iteration $ \cos y = y $ there is a fun method to find D using a hand calculator:

Enter any real positive number and press Cos button some 10 or 15 times until none of last digits changes... this is the Dottie number.