Add Polynomial class

[carlos-adir]

Currently the heavy.py calculates a 3D matrix with the coefficients of Spline Basis.

Adding a Polynomial class will encapsulate the horner method and hence the Spline Basis will store a 2D matrix of Polynomials