WIP: Ported code from R.
Create sigma matrix in one big swoosh.
Pass knots as constructor parameter.
Create method for calculating coefficients for a single x.