showed erickson's newton basis functions and mine are the same
added varargout for point order in calculate newton_basis
tried to plot point selection for 2-d case, still not sure what's wrong
changed colloc_poisson back to script form