Some progress tovards cubic spline calculation
Added new constructor to create tridiagonal matrix from supplied arrays
Added new exception
Cubic spline functions now forms tridiaganol matrix and solves it. Haven't tested yet