diff --git a/docs/make.jl b/docs/make.jl
index 76862401f..a869110f5 100644
--- a/docs/make.jl
+++ b/docs/make.jl
@@ -1,4 +1,5 @@
 # -*- coding: utf-8 -*-
+using Revise
 using Documenter
 using Documenter.DocMeta: setdocmeta!
 using DocumenterCitations
diff --git a/docs/src/Science/04-Continuum-Mechanics.md b/docs/src/Science/04-Continuum-Mechanics.md
index cc5094b5d..178ed3dab 100644
--- a/docs/src/Science/04-Continuum-Mechanics.md
+++ b/docs/src/Science/04-Continuum-Mechanics.md
@@ -386,6 +386,8 @@ Limiting cases are:
 
 ### 021
 
+![](https://www.youtube.com/embed/rrZgH7y15XU)
+
 In fact it is in this lecture that the stripping of dimensions of the equation as presented above in 020 is formalized; *do not even try to solve a problem before making it dimensionless*. The general procedure for making differential operators dimensionless is summarized as follows:
 
 $$
@@ -401,12 +403,16 @@ $$
 
 #bessel-function
 
+![](https://www.youtube.com/embed/Bi9jWgw_aQ8)
+
 *Always start any modeling with the simplest geometry that captures the basic features of the system being modeled; e.g. when expanding a solution of reaction-diffusion into exponential terms, it is worth noticing that they can be replaced by hyperbolic functions, and since $\sinh$ breaks the symmetry, that term may be eliminated already during constant identification from boundary conditions*.
 
 ### 023
 
 #dimless-peclet #plug-flow #convection-diffusion
 
+![](https://www.youtube.com/embed/LBSzqxjQl2k)
+
 Similarly to lecture 019, here we develop convection-diffusion equation instead; in this case the flux of a transported concentration $C$ is given by the following expression: 
 
 $$
@@ -443,4 +449,25 @@ In the limiting case of $\mathrm{Pe}\gg{1}$ diffusion is much slower than convec
 
 ### 024
 
+![](https://www.youtube.com/embed/328868nsFH4)
+
 ### 025
+
+![](https://www.youtube.com/embed/aEe10NkgI98)
+
+### 026
+
+#dimless-peclet #plug-flow 
+
+![](https://www.youtube.com/embed/hdkt3S78f8A)
+
+### 027
+
+### 028
+
+### 029
+
+### 030
+
+### 031
+