diff --git a/README.md b/README.md
new file mode 100644
index 0000000..5129a7f
--- /dev/null
+++ b/README.md
@@ -0,0 +1,15 @@
+# Simulink Model Usage
+
+This is the Simulink model used in:
+
+Chambers JD, Bornstein JC, Thomas EA (2011) Multiple Neural Oscillators and Muscle Feedback Are Required for the Intestinal Fed State Motor Program. *PLoS ONE* 6(5):e19597. doi:10.1371/journal.pone.0019597
+
+To run the model, in MATLAB click **File -> Open**, then select the file **"Chambers et al 2010.mdl"**. Then press play.
+
+The numerical values are set to produce the "control response" described in the paper. Changing the appropriate variables will produce all model figures in the paper.
+
+When the model is run it saves matrices (**Circular_Muscle, Excitatory_Motor_Neurons, Inhbitory_Motor_Neurons, Sensory_Neurons, Serotonin_release**) in the MATLAB workspace. Each matrix has two columns, the first column is the simulation time for each data point and the second column is the numerical value of the activity in the population of neurons, or muscle etc. These matrices can be plotted using the MATLAB plot function to reproduce the figures in the paper. The **Circular_Muscle** matrix can be passed to the **muscleresponse** function (e.g. `muscleresponse(Circular_Muscle)`) and this function will calculate the duration of the activity episodes and quiescent periods as reported in the paper.
+
+---
+
+2025-06-02: Converted README to Markdown.
\ No newline at end of file
diff --git a/README.txt b/README.txt
deleted file mode 100644
index 1f5c2b5..0000000
--- a/README.txt
+++ /dev/null
@@ -1,25 +0,0 @@
-This is the Simulink model used in:
-
-Chambers JD, Bornstein JC, Thomas EA (2011) Multiple Neural
-Oscillators and Muscle Feedback Are Required for the Intestinal Fed
-State Motor Program. PLoS ONE 6(5):e19597.
-doi:10.1371/journal.pone.0019597
-
-To run the model, in MATLAB click File->Open, then selected the file
-"Chambers et al 2010.mdl". Then press play.
-
-The numerical values are set to produce the "control response"
-described in the paper. Changing the appropriate variables will
-produce all model figures in the paper.
-
-When the model is run it saves matrixes (Circular_Muscle,
-Excitatory_Motor_Neurons, Inhbitory_Motor_Neurons, Sensory_Neurons,
-Serotonin_release) in the MATLAB workspace. Each matrix has two
-columns, the first column is the simulation time for each data point
-and the second column is the numerical value of the activity in the
-population of neurons, or muscle etc. These matrices can be plotted
-using the MATLAB plot function to reproduce the figures in the
-paper. The Circular_Muscle matrix can be passed to the muscleresponse
-function (e.g. muscleresponse(Circular_Muscle)) and this function will
-calculate the duration of the activity episodes and quiescent periods
-as reported in the paper.