The circular functions occur throughout mathematics. Fourier showed that, under very general assumptions, an arbitrary function can be decomposed into components each of which is a circular function. The functions get their name from their use in defining a circle in parametric form: if

then , the usual equation for a circle in Cartesian coordinates. In the figure, we plot the familiar sinusoid, which has a period of .

**Step Function & Square Wave Function**

Electrical engineers just love step-functions, which can be used to represent processes that are switched on at a specific time. Square-wave functions are also valuable in modelling processes that alternate between two states. One problem with these functions is that they have jump discontinuities, which can be difficult to handle in an analytical context where continuity and smoothness are expected. It is sometimes advantageous to replace them by functions that are continuous or even differentiable.

The step function can be represented by the function . For we put . Alternatively, we may write , again with . This function can be approximated by various elementary analytical functions, for example . This function becomes steeper near the origin, coming closer to the step function as increases.

Considering the square wave, the Fourier series converging to this gives approximations that converge to the correct solution. But the convergence is not uniform and there are spurious oscillations called Gibbs waves that can be problematical.

It is desirable to have a series of functions that converge to a square wave, and that have attractive analytical properties. The square-wave function can be approximated by a sine function. Indeed, if we plot the function we see that it is just a square-wave. To get a smooth approximation to this, we can consider the sine function raised to a small fractional power: . We must be careful: when is negative, this is imaginary, so we use

The first component looks after the sign while the second becomes closer to 1 as increases.

**Elliptic Functions**

The sine functions raised to small fractional powers approach the square wave function, but they are singular at the points where they vanish. We can obtain analytical approximations to the square wave by using Jacobian elliptic functions. The function is defined to be a solution of the ode

with appropriate initial conditions. For we have . As increases from 0 to 1, this function changes character, becoming flatter near its maxima and minima as increases. The notation is often used, with called the modulus.

The function has period , where is the complete elliptic integral of the first kind. To ensure functions having period of , we consider the functions . The results are plotted in the following figure for several values of .

Thus, by using Jacobian elliptic functions we obtain a beautiful sequence of analytical approximations to the square wave.

This is a nice way to square the circular functions!