Polytropic EOS, varying gamma

p = [gamma(rho) - 1] rho e
where gamma(rho) is a polynomial:
gamma(rho) = sumi=0order ci rhoi
type name
polytrop
parameters
polynomial_order
For each coefficient (0 to order, i.e. order + 1 coefficients)
   ci