diff --git a/doc/Eqs/fix_ttm_blast.jpg b/doc/Eqs/fix_ttm_blast.jpg
new file mode 100644
index 0000000000..d4ffe05f12
Binary files /dev/null and b/doc/Eqs/fix_ttm_blast.jpg differ
diff --git a/doc/Eqs/fix_ttm_blast.tex b/doc/Eqs/fix_ttm_blast.tex
new file mode 100644
index 0000000000..6871ed5cc6
--- /dev/null
+++ b/doc/Eqs/fix_ttm_blast.tex
@@ -0,0 +1,13 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+
+$$ 
+  {\vec F}_i = - \partial U / \partial {\vec r}_i + {\vec F}_{langevin} - \nabla P_e/n_{ion}
+$$
+
+\end{document}
+
+
+
+
diff --git a/doc/Eqs/fix_ttm_blast1.jpg b/doc/Eqs/fix_ttm_blast1.jpg
new file mode 100644
index 0000000000..1106c627cb
Binary files /dev/null and b/doc/Eqs/fix_ttm_blast1.jpg differ
diff --git a/doc/Eqs/fix_ttm_blast1.tex b/doc/Eqs/fix_ttm_blast1.tex
new file mode 100644
index 0000000000..642afe4300
--- /dev/null
+++ b/doc/Eqs/fix_ttm_blast1.tex
@@ -0,0 +1,13 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+
+$$ 
+  \nabla_x P_e = \left[\frac{C_e{}T_e(x)\lambda}{(x+\lambda)^2} + \frac{x}{x+\lambda}\frac{(C_e{}T_e)_{x+\Delta x}-(C_e{}T_e)_{x}}{\Delta x} \right]
+$$
+
+\end{document}
+
+
+
+
diff --git a/doc/Eqs/fix_ttm_ce.jpg b/doc/Eqs/fix_ttm_ce.jpg
new file mode 100644
index 0000000000..4f439e31e6
Binary files /dev/null and b/doc/Eqs/fix_ttm_ce.jpg differ
diff --git a/doc/Eqs/fix_ttm_ce.tex b/doc/Eqs/fix_ttm_ce.tex
new file mode 100644
index 0000000000..f787590b8f
--- /dev/null
+++ b/doc/Eqs/fix_ttm_ce.tex
@@ -0,0 +1,13 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+
+$$ 
+  C_e = C_0 + (a_0 + a_1 X + a_2 X^2 + a_3 X^3 + a_4 X^4) \exp (-(AX)^2)
+$$
+
+\end{document}
+
+
+
+
diff --git a/doc/Eqs/fix_ttm_mod.jpg b/doc/Eqs/fix_ttm_mod.jpg
new file mode 100644
index 0000000000..7cc67200f9
Binary files /dev/null and b/doc/Eqs/fix_ttm_mod.jpg differ
diff --git a/doc/Eqs/fix_ttm_mod.tex b/doc/Eqs/fix_ttm_mod.tex
new file mode 100644
index 0000000000..d0a313e676
--- /dev/null
+++ b/doc/Eqs/fix_ttm_mod.tex
@@ -0,0 +1,15 @@
+\documentclass[12pt]{article}
+
+\begin{document}
+
+$$ 
+  C_e \rho_e \frac{\partial T_e}{\partial t} = 
+  \bigtriangledown (\kappa_e \bigtriangledown T_e) - 
+  g_p (T_e - T_a) + g_s T_a' + \theta (x-x_{surface})I_0 \exp(-x/l_{skin})
+$$
+
+\end{document}
+
+
+
+