Low-pass Filtering and the Remarkable Integrals of Borwein and Borwein

In last week's post we looked at aspects of puzzles of the form ``What is the next number''. We are presented with a short list of numbers, for example $latex {1, 3, 5, 7, 9}&fg=000000$ and asked for the next number in the sequence. Arguments were given indicating why any number might be regarded as … Continue reading Low-pass Filtering and the Remarkable Integrals of Borwein and Borwein

The Power of the 2-gon: Extrapolation to Evaluate Pi

  Richardson's extrapolation procedure yields a significant increase in the accuracy of numerical solutions of differential equations. We consider his elegant illustration of the technique, the evaluation of $latex {\pi}&fg=000000$, and show how the estimates improve dramatically with higher order extrapolation. [This post is a condensed version of a paper in Mathematics Today (Lynch, 2003).] … Continue reading The Power of the 2-gon: Extrapolation to Evaluate Pi

Real Derivatives from Imaginary Increments

The solution of many problems requires us to compute derivatives. Complex step differentiation is a method of computing the first derivative of a real function, which circumvents the problem of roundoff error found with typical finite difference approximations. For finite difference approximations, the choice of step size $latex {h}&fg=000000$ is crucial: if $latex {h}&fg=000000$ is … Continue reading Real Derivatives from Imaginary Increments

Does Numerical Integration Reflect the Truth?

Many problems in applied mathematics involve the solution of a differential equation. Simple differential equations can be solved analytically: we can find a formula expressing the solution for any value of the independent variable. But most equations are nonlinear and this approach does not work; we must solve the equation by approximate numerical means. The … Continue reading Does Numerical Integration Reflect the Truth?

Dimension Reduction by PCA

We live in the age of ``big data''. Voluminous data collections are mined for information using mathematical techniques. Problems in high dimensions are hard to solve --- this is called ``the curse of dimensionality''. Dimension reduction is essential in big data science. Many sophisticated techniques have been developed to reduce dimensions and reveal the information … Continue reading Dimension Reduction by PCA

The Rambling Roots of Wilkinson’s Polynomial

Finding the roots of polynomials has occupied mathematicians for many centuries. For equations up to fourth order, there are algebraic expressions for the roots. For higher order equations, many excellent numerical methods are available, but the results are not always reliable. James Wilkinson (1963) examined the behaviour of a high-order polynomial $latex \displaystyle p(x,\epsilon) = … Continue reading The Rambling Roots of Wilkinson’s Polynomial

Bouncing Billiard Balls Produce Pi

There are many ways of evaluating $latex {\pi}&fg=000000$, the ratio of the circumference of a circle to its diameter. We review several historical methods and describe a recently-discovered and completely original and ingenious method. Historical Methods Archimedes used inscribed and circumscribed polygons to deduce that $latex \displaystyle \textstyle{3\frac{10}{71} < \pi < 3\frac{10}{70}} &fg=000000$ giving roughly … Continue reading Bouncing Billiard Balls Produce Pi

Marvellous Merchiston’s Logarithms

Log tables, invaluable in science, industry and commerce for 350 years, have been consigned to the scrap heap. But logarithms remain at the core of science, as a wide range of physical phenomena follow logarithmic laws  [TM103 or search for “thatsmaths” at irishtimes.com]. The method of logarithms was first devised by John Napier, 8th Laird … Continue reading Marvellous Merchiston’s Logarithms

French Curves and Bézier Splines

A French curve is a template, normally plastic, used for manually drawing smooth curves. These simple drafting instruments provided innocent if puerile merriment to generations of engineering students, but they have now been rendered obsolete by computer aided design (CAD) packages, which enable us to construct complicated curves and surfaces using mathematical functions called Bézier … Continue reading French Curves and Bézier Splines