Viscous damping is commonly discussed in beginning differential equations and physics texts but dry friction or Coulomb friction is not despite dry friction being encountered in many physical applications. One reason for avoiding this topic is that the equations involve a jump discontinuity in the damping term. In this article, we adopt an energy…

Quadratic friction involves a discontinuous damping term in equations of motion in order that the frictional force always opposes the direction of the motion. Perhaps for this reason this topic is usually omitted from beginning texts in differential equations and physics. However, quadratic damping is more realistic than viscous damping in many…

In this article we examine 2 x 2 first-order systems of ordinary differential equations and show how to identify separatrices for phase plane portraits when the system has a saddle point critical value. We describe how to use a computer algebra system to generate trajectories from contour plots, when possible, and determine the equation of the…

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…

Finding a periodic solution to a nonlinear ordinary differential equation is in general a difficult task. Only in a very few cases can direct methods be applied to an equation to find initial values leading to a solution of the corresponding initial value problem that is periodic. Oscillatory periodic solutions have such practical importance that…

Through numerical investigations, various examples of the Duffing type forced spring equation with epsilon positive, are studied. Since [epsilon] is positive, all solutions to the associated homogeneous equation are periodic and the same is true with the forcing applied. The damped equation exhibits steady state trajectories with the interesting…

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…

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…

This note discusses the boundary in the frequency--amplitude plane for boundedness of solutions to the forced spring Duffing type equation. For fixed initial conditions and fixed parameter [epsilon] results are reported of a systematic numerical investigation on the global stability of solutions to the initial value problem as the parameters F and…

This paper discusses a result of Li and Shen which proves the existence of a unique periodic solution for the differential equation x[dots above] + kx[dot above] + g(x,t) = [epsilon](t) where k is a constant; g is continuous, continuously differentiable with respect to x , and is periodic of period P in the variable t; [epsilon](t) is continuous…

This paper compliments two recent articles by the author in this journal concerning solving the forced harmonic oscillator equation when the forcing is periodic. The idea is to replace the forcing function by its Fourier series and solve the differential equation term-by-term. Herein the convergence of such series solutions is investigated when…

Described is an approach to the derivation of numerical integration formulas. Students develop their own formulas using polynomial interpolation and determine error estimates. The Newton-Cotes formulas and error analysis are reviewed. (KR)