The maths of multi curve stripping

Overview

Multi Dimensional Newton

we are trying to solve $$ f(x) = 0 $$ Where \( f(x) \) is a vector valued function with each element the value of one of the benchmark instruments, and \( x \) is the column vector of curve pillar point values (either discount factors or rates)

This can be solved using Broyden’s method in MathNet or with the QuantSA implementation of a multi dimensional Newton.

The QuantSA version proceeds by approximating \(f(x)\) as linear: $$ f(x+\delta) \approx f(x) + J(x)\delta $$ Where \(J\) is the Jacobian matrix with $$J_{ij} = \frac{\partial f_i}{\partial x_j}$$ since we require $$f(x+\delta)=0$$ we estimate $$\delta = -J(x)^{-1}f(x)$$