Through numerical investigations, we study examples of the forced quadratic spring equation [image omitted]. By performing trial-and-error numerical experiments, we demonstrate the existence of stability boundaries in the phase plane indicating initial conditions yielding bounded solutions, investigate the resonance boundary in the [omega]…

We report on a study of the forced van der Pol equation x + [epsilon](x[superscript 2] - 1)x + x = F cos[omega]t, by solving numerically the differential equation for a variety of values of the parameters [epsilon], F and [omega]. In doing so, many striking and interesting trajectories can be discovered and phenomena such as frequency entrainment,…

We discuss the boundary in the Poincare phase plane for boundedness of solutions to spring model equations of the form [second derivative of]x + x + epsilonx[superscript 2] = Fcoswt and the [second derivative of]x + x + epsilonx[superscript 3] = Fcoswt and report the results of a systematic numerical investigation on the global stability of…

We investigate the pendulum equation [theta] + [lambda][squared] sin [theta] = 0 and two approximations for it. On the one hand, we suggest that the third and fifth-order Taylor series approximations for sin [theta] do not yield very good differential equations to approximate the solution of the pendulum equation unless the initial conditions are…

The phenomenon of nonlinear resonance (sometimes called the "jump phenomenon") is examined and second-order van der Pol plane analysis is employed to indicate that this phenomenon is not a feature of the equation, but rather the result of accumulated round-off error, truncation error and algorithm error that distorts the true bounded solution onto…

Given three points in the plane, interest is in the locus of all points for which the sum of the distances to the given points is a prescribed constant. These curves turn out to be sixth degree polynominals in x and y , and thus are complicated. However, it turns out that often there is a point, within the triangle formed by the three given…

Non-linear second-order differential equations whose solutions are the elliptic functions "sn"("t, k"), "cn"("t, k") and "dn"("t, k") are investigated. Using "Mathematica", high precision numerical solutions are generated. From these data, Fourier coefficients are determined yielding approximate formulas for these non-elementary functions that are…