plot[ k=18e9 * e^(-13000/t), {t,590,660} ]