Programming with Chebfun. Case study: Richards equation (opens in new tab)

The Chebfun software system is a Matlab extension essentially based on representations of (piece-wise) smooth one-variable functions by expansions in Chebyshev polynomials. One of Chebfun's attractive features is the ability to provide solutions to nonlinear boundary value problems (BVP) with accuracy close to the machine precision. This is done by the chebop class which provides automatic solutions by performing linearizations with a Newton met...