usage: [bp,wf]=ncrule(m);
This function returns the Newton-Coates base points and weight factors
up to an 8 point Newton-Coates formula.
m -- number of Newton-Coates points (integrates an mth order
polynomial exactly (or an (m+1)th order for even m))