A Polynomial Chaos Method for Dispersive Electromagnetics