Fourier Series Approximation

The Fourier Series

Given a function \(f(x)\) defined on \([a,b]\), let \(L=\tfrac{b-a}{2}\) be the half-period and \(c=\tfrac{a+b}{2}\) the midpoint. The Fourier series of \(f\) with \(n\) modes is

The Fourier coefficients are defined by

The constant \(a_0/2\) is the average value of \(f\) over \([a,b]\). Each pair \((a_k,b_k)\) captures the amplitude of the \(k\)-th harmonic. All coefficients are computed numerically using composite 10-point Gauss–Legendre quadrature for high accuracy without requiring a closed-form antiderivative.

Convergence and the Gibbs Phenomenon

By the Dirichlet theorem, if \(f\) is piecewise smooth on \([a,b]\), then \(S_n(x)\to f(x)\) at every point of continuity as \(n\to\infty\). At a jump discontinuity \(x_0\), the series converges to the average of the one-sided limits:

Near a jump you will observe the Gibbs phenomenon: an overshoot of approximately \(9\%\) of the jump magnitude that does not diminish as \(n\to\infty\), but instead concentrates into a narrower and narrower spike. Try the step function or square wave examples to see this clearly.

Relative \(L^2\) Error

The plot title reports the relative \(L^2\) error, measuring how close \(S_n\) is to \(f\) in the root-mean-square sense:

By Parseval's theorem, this error tends to zero as \(n\to\infty\) for any square-integrable \(f\) - even when pointwise convergence fails at jump discontinuities.

Expression Syntax

• Standard arithmetic: +, -, *, /, ^. Always write multiplication explicitly: 2*x not 2x.
• Functions: sin, cos, exp, abs, sqrt, log, and all math.js built-ins.
• Built-in specials: sawtooth(x) and square(x) (period \(2\pi\)).
• Domain accepts plain numbers or \(\pi\)-expressions: -pi, 2*pi, -3*pi/2.

Piecewise Functions

Tick the Piecewise box next to Load example to switch the function input into structured pieces. The default layout has two pieces and one breakpoint; click + Add piece to extend it to three, four, or more, and − Remove last piece to shrink it back. With \(N\) pieces and breakpoints \(b_1

• Piece 1 is used for \(x\le b_1\).
• Piece \(i\) (for \(1
• Piece \(N\) is used for \(x>b_{N-1}\).

Each piece is parsed independently as a normal expression, and each breakpoint accepts a number or a pi-expression (e.g. pi/2, 2*pi). Examples:

• Piece 1 (x+2)^2, breakpoint \(0\), Piece 2 8-x  (quadratic then linear)
• Piece 1 -1, breakpoint \(0\), Piece 2 1  (step function - ideal for Gibbs phenomenon)
• Piece 1 -1, \(-pi/2\), Piece 2 sin(x), \(pi/2\), Piece 3 1  (clamped sine)

Untick Piecewise to enter a single closed-form expression in the function box instead.

How to Use

• Select an example from the dropdown or type your own function. The plot previews immediately.
• Set the domain \([x_\min,x_\max]\) and the number of Fourier modes \(n\) (up to 200).
• Adjust the animation speed (milliseconds per frame).
• Click Start to animate the convergence term by term. Click Stop to pause; Start restarts from the beginning.
• Press Enter in any input field to start immediately.