added spectrum regression to planck.py, EM derivation in writeup