5
$\begingroup$

Is there an analytical solution to the nonlinear ODE $\frac{dx}{d\theta} = -\sqrt{\frac{x^2}{4\cos^2\theta} - \cos^2\theta}$ over $\theta \in [0, \pi/2]$ with initial condition $x(0) = 2$? Using the substitution $y = \sin\theta$, I was able to transform this into $\frac{dx}{dy} = -\sqrt{\frac{x^2-4(1-y^2)^2}{4(1-y^2)^2}},$ but couldn't get much further.

Context: On a typical globe, one can only see half of the earth at a time. Waldo Tobler's slides on "Unusual Map Projections" mention the following amusing solution to this problem: simply wrap the earth around the globe twice (see image).

In trying to find an aesthetically pleasing conformal version of this construction by changing the shape of the globe, I found myself needing to solve the above ODE to determine its cross section; specifically, mapping latitude $\theta$ to distance from the polar axis, $x(\theta)$.

It can be solved numerically, of course, and that's what I did. I'm just curious about whether an analytical solution exists. The numerical solution looks a lot like $x(\theta) = 2\sqrt{\cos\theta}$, but isn't exactly the same.

In case you're curious, here is a render of the final result of my hacking on all this.

  • 0
    @Ragib: I'm not putting too much faith in WolframAlpha here, because if you take out the 4 in the denominator of the first equation and set $x(0) = 1$, there is a closed form solution $x(\theta) = \cos\theta$ which WolframAlpha cannot find either.2011-09-20

1 Answers 1

3

One can apply the following procedure showing that this equation is indeed amenable to an exact analytical treatment. The first thing to note is that it can be rewritten in the following form

$\left(\frac{dx}{d\theta}\right)^2-\frac{x^2}{4\cos^2\theta}+\cos^2\theta=0.$

Now, let us consider the Hamiltonian system

$H=p^2-\frac{x^2}{4\cos^2\theta}+\cos^2\theta.$

It is easy to note that, for $p(0)=0$ and $x(0)=2$, one has $H(\theta=0)=0$. So, a solution of the Hamilton equations (here dot means derivation with respect to $\theta$)

$\dot x=\frac{\partial H}{\partial p}=2p$ $\dot p=-\frac{\partial H}{\partial x}=\frac{x}{2\cos^2\theta}$

gives also a solution to the equation we started from. From Hamilton equations we get

$\ddot x=\frac{x}{\cos^2\theta}$

with the initial conditions $\dot x(0)=0$ and $x(0)=2$. This second order equation has an explicit analytical solution, through hypergeometric functions $\left._2F_1\right.$ and $\cos t$, that can be easily obtained using Mathematica or your other preferred symbolic manipulation software.

  • 0
    There is a sign to be changed in the Hamiltonian: -Cos[x]^2 should be Cos[x]^2, but this is irrelevant. These points correspond to the points where $\cos(x)=0$ and I do not know if it is just a numerical shortcoming or something deeper. I should expect $x(\pi/2)=0$.2011-12-08