> moved resizing of A and b to initialize()
> added loop in initialize() to set all terms of A and b to zero
> added global variable 'order', which generalizes the least-squares matrices to using n-dimensional polynomials
> completely restructured summation() to allow for generality in matrix definitions