For numeric calculations we approximate the density by a piecewise constant function:
In order to ensure smoothness of the plasma density we represent it by a natural cubic spline passing through control points and calculate by interpolation at . Fixing and while varying only , we get by a simple matrix multiplication: . In order to obtain the spline matrix we calculate the cubic spline passing through and set .
We finally arrive at a matrix equation relating and the absorbance :