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.
The following equations govern the time evolution of this system:
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.



  • [Snapshot]
  • [Snapshot]
  • [Snapshot]
  • [Snapshot]


For a more detailed discussion of this model please see This paper by Jorge Sotomayor describes this model and some possible improvements that might allow one to better understand the parameter space. The idea for this model was taken from H. K. Wilson, Ordinary Differential Equations, Reading, MA.: Addison–Wesley, 1971. One may also refer to Jorge Sotomayor, Liçoes de Equaçoes Diferenciais Ordinárias, Rio de Janeiro: IMPA, 1979 (in Portuguese).
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.