def get_rmg_polynomial(las, has, temps=[200., 1000., 1000., 3000.]):
    """
    Return NASA polynomial as a dictionary in RMG format:

    NASA Polynomial, seven or nine coefficients, Tmax and Tmin = valid temperature range
    polynomials = [{'coeffs':[2.3443,0.00798042,-1.94779e-05,2.0157e-08,-7.37603e-12,-917.924,0.683002], 'Tmin':(200,'K'), 'Tmax':(1000,'K')},
                      {'coeffs':[2.93283,0.000826598,-1.46401e-07,1.54099e-11,-6.88796e-16,-813.056,-1.02432], 'Tmin':(1000,'K'), 'Tmax':(6000,'K')}]
                      Tmin = (200,'K')
                      Tmax = (6000,'K')
                      NASAPolynomial = {'polynomials':polynomials,'Tmin':Tmin,'Tmax':Tmax}
    """
    p = [{'coeffs': las, 'Tmin': (temps[0], 'K'), 'Tmax':(temps[1], 'K')},
         {'coeffs': has, 'Tmin': (temps[2], 'K'), 'Tmax':(temps[3], 'K')}]

    return {'polynomials': p, 'Tmin': (min(temps), 'K'), 'Tmax': (max(temps), 'K')}