[NL] Replace order with n_points for correctness.
[NL] Rename order to n_points. Documentation.
[ML] Add Newton-Cotes quadrature rules.

The coefficients are for the closed (including end points) type for the interval [-1, 1].
md gauss legendre correct npoints
md newton cotes