diff --git a/docs/source/images/taguchi_process.png b/docs/source/images/taguchi_process.png
new file mode 100644
index 00000000..176ab19f
Binary files /dev/null and b/docs/source/images/taguchi_process.png differ
diff --git a/docs/source/index.rst b/docs/source/index.rst
index 702a4dc0..02446957 100644
--- a/docs/source/index.rst
+++ b/docs/source/index.rst
@@ -37,6 +37,7 @@ gprMax User Guide
     :caption: User libraries
 
     user_libs_antennas
+    user_libs_opt_taguchi
     user_libs_austinman
 
 .. toctree::
diff --git a/docs/source/references.rst b/docs/source/references.rst
index 264f43b1..ee90ade4 100644
--- a/docs/source/references.rst
+++ b/docs/source/references.rst
@@ -25,6 +25,8 @@ References
 .. [TUR1997] Turcotte, D. L. (1997). Fractals and chaos in geology and geophysics. Cambridge university press. (http://dx.doi.org/10.1017/cbo9781139174695)
 .. [VIA2005] Vial, A., Grimault, A. S., Macías, D., Barchiesi, D., & de La Chapelle, M. L. (2005). Improved analytical fit of gold dispersion: Application to the modeling of extinction spectra with a finite-difference time-domain method. Physical Review B, 71(8), 085416. (http://dx.doi.org/10.1103/physrevb.71.085416)
 .. [WAR2011] Warren, C., & Giannopoulos, A. (2011). Creating finite-difference time-domain models of commercial ground-penetrating radar antennas using Taguchi’s optimization method. Geophysics, 76(2), G37-G47. (http://dx.doi.org/10.1190/1.3548506)
+.. [WEN2007a] Weng, W. C., Yang, F., & Elsherbeni, A. (2007). Electromagnetics and antenna optimization using Taguchi’s method. Synthesis Lectures on Computational Electromagnetics, 2(1), 1-94.
+.. [WEN2007b] Weng, W. C., Yang, F., & Elsherbeni, A. Z. (2007). Linear antenna array synthesis using Taguchi's method: A novel optimization technique in electromagnetics. Antennas and Propagation, IEEE Transactions on, 55(3), 723-730.
 .. [WHI2009] Whittow, W. G., & Edwards, R. M. (2009). Effects of averaging procedures for electrical properties at the interface of dissimilar tissues in the human head with finite-difference time-domain modelling. Science, Measurement & Technology, IET, 3(1), 51-58.
 .. [YEE1966] Yee, K. S. (1966). Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media. IEEE Trans. Antennas Propag, 14(3), 302-307. (http://dx.doi.org/10.1109/TAP.1966.1138693)
 
diff --git a/docs/source/user_libs_antennas.rst b/docs/source/user_libs_antennas.rst
index 5088444b..8017b966 100644
--- a/docs/source/user_libs_antennas.rst
+++ b/docs/source/user_libs_antennas.rst
@@ -1,9 +1,8 @@
-
 User libraries is a sub-package where useful Python modules contributed by users are stored.
 
-***********
-antennas.py
-***********
+********
+Antennas
+********
 
 .. code-block:: python
 
diff --git a/docs/source/user_libs_austinman.rst b/docs/source/user_libs_austinman.rst
index f523deb4..af585653 100644
--- a/docs/source/user_libs_austinman.rst
+++ b/docs/source/user_libs_austinman.rst
@@ -1,5 +1,7 @@
 User libraries is a sub-package where useful Python modules contributed by users are stored.
 
+**UNDER CONSTRUCTION**
+
 ***********************
 AustinMan & AustinWoman
 ***********************
diff --git a/docs/source/user_libs_opt_taguchi.rst b/docs/source/user_libs_opt_taguchi.rst
new file mode 100644
index 00000000..9c82c8ad
--- /dev/null
+++ b/docs/source/user_libs_opt_taguchi.rst
@@ -0,0 +1,46 @@
+User libraries is a sub-package where useful Python modules contributed by users are stored.
+
+**UNDER CONSTRUCTION**
+
+**********************
+Optimisation - Taguchi
+**********************
+
+.. code-block:: python
+
+    # Copyright (C) 2015-2016, Craig Warren
+    #
+    # This module is licensed under the Creative Commons Attribution-ShareAlike 4.0 International License.
+    # To view a copy of this license, visit http://creativecommons.org/licenses/by-sa/4.0/.
+    #
+    # Please use the attribution at http://dx.doi.org/10.1190/1.3548506
+
+The module features an optimisation technique based on Taguchi's method. It allows the user to define parameters in an input file and optimise their values based on a user-defined fitness function.
+
+Taguchi's method
+================
+
+Taguchi's method is based on the concept of Orthogonal Array (OA) and has the following advantages:
+
+* Simple to implement
+* Effective in reduction of experiments
+* Fast convergence speed
+* Global optimum results
+* Independence from initial values of optimisation parameters
+
+Details of Taguchi's method in the context of electromagnetics can be found in [WEN2007a]_ and [WEN2007b]_. The process by which Taguchi's method optimises parameters is illustrated in the following figure.
+
+.. figure:: images/taguchi_process.png
+    :width: 300 px
+
+    Process associated with Taguchi's method.
+
+Implementation
+==============
+
+The module will select from 2 pre-built OAs (http://neilsloane.com/oadir/) depending on the number of parameters to optimise. Currently, up to 7 independent parameters can be optimised, although a method to construct OAs of any size is under testing.
+
+Fitness functions
+-----------------
+
+A fitness function is required
\ No newline at end of file