diff --git a/SFEM/README.md b/SFEM/README.md
index 2c1aef329661313750d6c7a2a010e9ce0b4c0cb0..6c0368ea58abae840d33615c98fe1836fafb0029 100644
--- a/SFEM/README.md
+++ b/SFEM/README.md
@@ -110,7 +110,7 @@ In this example, we consider a multiscale lattice simulation in which the lattic
 
 ### Generate the Gauss-points list and mesh structure for the random field generator
 
-From the [rnnRF](./rnnRF/) directory, run
+This steps requires cm3Libraries (http://www.ltas-cm3.ulg.ac.be/openSource.htm). From the [rnnRF](./rnnRF/) directory, run
 
 ``` bash
 python3 generateGPInfo.py
@@ -157,6 +157,20 @@ gmsh rubics.msh randomFields/RandField_*.msh
 ```
 to vizualize the RF discretization compatible with the geometry.
 
+### Run a MC Stochastic finite element simulation
+
+This steps requires cm3Libraries (http://www.ltas-cm3.ulg.ac.be/openSource.htm). From the [rnnRF](./rnnRF/) directory, run
+
+```bash
+python3 runMCTest.py 
+```
+After execution the ouput files are store in ```'rnnRF/results'```:
+* The script saves the simulations results in ```'rnnRF/results/RandField_And_GP_X.csv/'```
+
+To vizualize the different resulting stress/strain curves
+``` bash
+python3 plotMC.py
+```
 
 ## Disclaimer