This post illustrates the working of Ramanujan's generating functions for solving Euler's diophantine equation \(a^3+b^3=c^3+d^3\) as described by Andrews and Berndt in "Ramanujan's Lost Notebook, Part IV", pp 199 - 205 (Section 8.5). The text is available from Springer. Ramanujan's result is that if \[ …

Steffensen's method is based on Newton's iteration for solving a non-linear equation \(f(x)=0\): \[ x\leftarrow x-\frac{f(x)}{f'(x)} \] Newton's method can fail to work in a number of ways, but when it does work it displays qudratic convergence; the number of correct signifcant figures roughly doubling at each step. …

As is well known, tanh-sinh quadrature takes an integral \[ \int_{-1}^1f(x)dx \] and uses the substitution \[ x = g(t) = \tanh\left(\frac{\pi}{2}\sinh t\right) \] to transform the integral into \[ \int_{-\infty}^{\infty}f(g(t))g'(t)dt. \] The reason this works so well is that the derivative \(g'(t)\) dies away at a …

An article by Bailey, Jeybalan and LI, "A comparison of three high-precision quadrature schemes", and available online here, compares Gauss-Legendre quadrature, tanh-sinh quadrature, and a rule where the nodes and weights are given by the error function and its integrand respectively. However, Nick Trefethen …

Note: This blog post is mainly computational, with a hint of proof-oriented mathematics here and there. For a more in-depth analysis, read the excellent article "Gauss, Landen, Ramanujan, the Arithmetic-Geometric Mean, Ellipses, pi, and the Ladies Diary" by Gert Akmkvist and Bruce Berndt, in The American …

The purpose of this post will be to see if we can implement the algorithm in Julia, and thus leverage Julia's very fast execution time. We are working with polynomials defined on nilpotent variables, which means that the degree of any generator in a polynomial term will be 0 or 1. Assume that our generators are indexed …

Introduction We know that there is no general sub-exponential algorithm for computing the permanent of a square matrix. But we may very reasonably ask -- might there be a faster, possibly even polynomial-time algorithm, for some specific classes of matrices? For example, a sparse matrix will have most terms of the …

As I discussed in my last blog post, the permanent of an \(n\times n\) matrix \(M=m_{ij}\) is defined as \[ \text{per}(M)=\sum_{\sigma\in S_n}\prod_{i=1}^nm_{i,\sigma(i)} \] where the sum is taken over all permutations of the \(n\) numbers \(1,2,\ldots,n\). It differs from the better known determinant in having no sign …

Runge's phenomenon says roughly that a polynomial through equally spaced points over an interval will wobble a lot near the ends. Runge demonstrated this by fitting polynomials through equally spaced point in the interval \([-1,1]\) on the function \[ \frac{1}{1+25x^2} \] and this function is now known as "Runge's …

Introduction and the problem The SIR model for spread of disease was first proposed in 1927 in a collection of three articles in the Proceedings of the Royal Society by Anderson Gray McKendrick and William Ogilvy Kermack; the resulting theory is known as Kermack–McKendrick theory; now considered a subclass of a more …