deletions | additions       diff --git a/sectionAnalysis__sub.tex b/sectionAnalysis__sub.tex  index e8a1b98..e7c2c7b 100644  --- a/sectionAnalysis__sub.tex  +++ b/sectionAnalysis__sub.tex 
...
B_{mod} =B_0 \frac{2 h \nu^{3+\beta}}{c^3}\frac{1}{e^{\frac{h\nu}{kT}}-1}  \end{equation*}  where {\bf T is temperature, $\nu$ is frequency and}  $\beta${\bf, the emissivity index,} is assumed to be 2 \citep{Desert2008}. B$_0$ and T were taken as free parameters and a Levenberg{\textendash}Marquardt algorithm was used to find the best-fit. The fitting was done to the flux density in Jy, so B$_0$ carries these units. The total column density{\bf, N$_{tot}$,} was calculated following \citet{Miettinen_2010}. Briefly, we used the following relations: \begin{eqnarray*}  N_{tot} = \frac{I}{B_\nu \mu m_H \kappa R_d}\nonumber\\