Refactor model equations in README to LaTeX

Updated mathematical equations to LaTeX format for clarity.

Author: mironovb

README.md

@@ -49,21 +49,27 @@ Edit parameters in **`scripts/main_ev_aex.py`** (or create your own script):
 * **Programs:** `grad_start`, `grad_end`, `I_load`, `I_elute` and `salt_profile(t)`
 * **Feed:** `Cfeed_total`, `frac`, `feed_profile(t)` (e.g., 5 min load at 1 mL/min)
 
-## What the model solves (brief)
-
-Mobile phase for species *i*:
-[
-\partial_t C_i = -u,\partial_z C_i + D_{\mathrm{ax}},\partial_{zz} C_i - \frac{1-\varepsilon}{\varepsilon},\partial_t Q_i
-]
-Wall phase (competitive Langmuir):
-[
-\partial_t Q_i = k_{a,i} C_i\Big(Q_{\max}-\sum_j Q_j\Big) - k_{d,i}Q_i
-]
-Optional affinity map:
-[
-K_i \equiv \frac{k_{a,i}}{k_{d,i}} = K^\circ \exp!\big(\gamma,|P_i|,f(I)\big)
-]
-*(Use constant (K_i) by setting (\gamma=0).)*
+## What the model solves
+
+**Mobile phase for species \(i\):**
+$$
+\partial_t C_i \;=\; -u\,\partial_z C_i \;+\; D_{\mathrm{ax}}\,\partial_{zz} C_i
+\;-\; \frac{1-\varepsilon}{\varepsilon}\,\partial_t Q_i .
+$$
+
+**Wall phase (competitive Langmuir):**
+$$
+\partial_t Q_i \;=\; k_{a,i}\,C_i\Big(Q_{\max}-\sum_j Q_j\Big) \;-\; k_{d,i}\,Q_i .
+$$
+
+**Optional affinity map:**
+$$
+K_i \;\equiv\; \frac{k_{a,i}}{k_{d,i}}
+\;=\; K^\circ \exp\!\big(\gamma\,|P_i|\,f(I)\big).
+$$
+
+Use constant \(K_i\) by setting \(\gamma=0\).
+
 
 ## Repo layout