# Focal Attraction with Rotation

Focal Attraction with Rotation

This Demonstration uses Mathematica's NDSolve to plot trajectories of the Markus–Wilson model of focal attraction on a rotating medium. This dynamical system is parameterized by four positive parameters, , , , and , which correspond to physical parameters of the model: describes the rotational speed, quantifies the radial drift in the model (the default is zero, i.e., no drift), is the magnitude of the focal attraction whose source is at in Euclidean coordinates, where is the size of the rotating disk. To see how the trajectory develops through time (and hopefully its asymptotic behavior), autorun through the parameter. For large this may result in slow rendering of the smoothed curve.

ω

δ

ν

R

ω

δ

ν

(R,0)

R

T

max

T

max

The following equations govern the time evolution of this system:

dx

dt

δ-+x

2

R

2

x

2

y

2

R

2

x

2

y

ν(R-x)

2

(R-x)

2

y

dy

dt

δ-+y

2

R

2

x

2

y

2

R

2

x

2

y

νy

2

(R-x)

2

y

The solution may not converge rapidly enough to be meaningful for all parameter combinations, in which case warnings will arise. We have attempted to remove most of the offending parameter combinations. However, the parameter space is complicated and one also does not want to remove too much of the content.