Move colorspace code to color package
Relax tolerance of rootfinder for Gauss nodes

Nodes for some high orders could not be found to the previous level of
precision.  New level allows nodes to be found for at least all orders
up to 500.
Import LG element and quadrature weights

This provides support for high-order numerical integration, which may be
used in the future.
Add BlackBody class