19.01.2015 Views

My title - Departamento de Matemática da Universidade do Minho

My title - Departamento de Matemática da Universidade do Minho

My title - Departamento de Matemática da Universidade do Minho

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

6 FLOWS 43<br />

Counter-example.<br />

Both the curves x(t) = 0 and x(t) = t 3 solve the equation<br />

ẋ = 3x 2/3<br />

with initial condition x(0) = 0. The problem here is that the velocity field v(x) = 3x 2/3 , although<br />

continuous, is not differentiable and not even Lipschitz at the origin. You may notice that the<br />

solution starting, for example, at x 0 = 1 reaches (or better comes from) the singular point x − = 0<br />

in finite time, since<br />

∫ 0<br />

t(x − ) − t(x 0 ) =<br />

1<br />

1 3 y−2/3 dy<br />

= −1 .<br />

One-dimensional Newtonian motion in a time in<strong>de</strong>pen<strong>de</strong>nt force field. The onedimensional<br />

motion of a particle of mass m subject to a force F (x) that <strong>do</strong>es not <strong>de</strong>pend on time<br />

is <strong>de</strong>scribed by the Newton equation<br />

mẍ = −U ′ (x) ,<br />

where the potential U(x) = − ∫ F (x)dx is some primitive of the force. The total energy<br />

E (x, ẋ) = 1 2 mẋ2 + U(x)<br />

(which of course is <strong>de</strong>fined up to an arbitrary additive constant) of the system is a constant of<br />

the motion, i.e. is constant along solutions of the Newton equation. In particular, once a value E<br />

of the energy is given (<strong>de</strong>pending on the initial conditions), the motion takes place in the region<br />

where U(x) ≤ E, since the kinetic energy 1 2 mẋ2 is non-negative. Conservation of energy allows to<br />

reduce the problem to the first or<strong>de</strong>r ODE<br />

ẋ 2 = 2 (E − U(x)) ,<br />

m<br />

which has the unpleasant feature to be quadratic in the velocity ẋ. Meanwhile, if we are interested<br />

in a one-way trajectory going from some x 0 to x, say with x > x 0 , we may solve for ẋ and find the<br />

first or<strong>de</strong>r autonomous ODE<br />

√<br />

2<br />

ẋ = (E − U(x)) .<br />

m<br />

There follows that the time nee<strong>de</strong>d to go from x 0 to x is<br />

t(x) =<br />

∫ x<br />

x 0<br />

dy<br />

√<br />

.<br />

2<br />

m (E − U(y))<br />

The inverse function of the above t(x)<br />

√<br />

will give the trajectory x(t) with initial position x(0) = x 0<br />

2<br />

and initial positive velocity ẋ(0) =<br />

m (E − U(x 0)), at least for sufficiently small times t.<br />

6.3 Exponential<br />

The exponential. The exponential function, according to Walter Rudin “the most important<br />

function in mathematics” ([Ru87], 1st line of page 1), is the unique solution of the (autonomous)<br />

differential equation<br />

ẋ = x<br />

with initial condition x(0) = 1. Actually, it is convenient to complexify time, i.e. take z = t+iθ ∈ C<br />

with t, θ ∈ R, and <strong>de</strong>fine the exponential as the power series<br />

exp(z) = 1 + x + x2<br />

2 + x3<br />

6 + x4<br />

24 + · · · = ∑ z n<br />

n!<br />

n≥0<br />

Since lim sup n→∞ (1/n!) 1/n = 0, the radius of convergence is ∞, hence the power series <strong>de</strong>fines an<br />

entire function, i.e. a holomorphic function exp : C → C. Deriving each term of the series, we

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!