diff --git a/assets/default view.png b/assets/default view.png
new file mode 100644
index 0000000..6a77c6e
Binary files /dev/null and b/assets/default view.png differ
diff --git a/assets/tec3pessen.png b/assets/tec3pessen.png
new file mode 100644
index 0000000..bddb90c
Binary files /dev/null and b/assets/tec3pessen.png differ
diff --git a/documentation/PID tuning.md b/documentation/PID tuning.md
new file mode 100644
index 0000000..c3933e6
--- /dev/null
+++ b/documentation/PID tuning.md	
@@ -0,0 +1,73 @@
+# PID Tuning
+
+## Real time plot
+
+When tuning Thermostat PID parameters, it is helpful to view the temperature, PID output and other data in the form of a real time graph.
+
+To use the Python real-time plotting utility, run
+
+```shell
+python pytec/plot.py
+```
+
+![default view](/assets/default%20view.png)
+
+## Temperature Setpoints and Thermal Load
+
+A PID controller with the same set of PID parameters may not work identically across all temperatures, especially when comparing the performance of a TEC module cooling a load versus heating a load. This is due to self ohmic heating of the TEC module aiding efficiency when  heating, but harming efficiency when cooling.
+
+When a PID loop is expected to operate the TEC in both heating and cooling modes, it is important to verify the loop performance in both modes.
+
+For systems expected to operate at a narrow range of temperatures, it is a good idea to tune the PID loop at the temperature region of interest.
+
+The same is also true for controlling loads that are expected to produce heat, e.g. laser cooling blocks. Testing the loop performance across varying amount of thermal load is needed to ensure stability in operation.
+
+## Manual Tuning
+
+Below are some general guidelines for manually tuning PID loops. Note that every system is different, and some of the values mentioned below may not apply to all systems.
+
+1. To start the manual tuning process, set the kp, ki and kd parameters to 0.
+
+2. Begin by increasing kP until the temperature begins to oscillate. Offset between the target temperature and the actual temperature can be ignored for now.
+
+3. Reduce kP by 30%, increase ki until the offset between target and actual temperature is eliminated.
+
+4. Increase kd until the maximum allowable amount of overshoot is observed.
+
+5. Some tweaking will be needed to obtain the desired result, especially when trying to balance between minimizing overshoot and maximizing response speed.
+
+## Auto Tuning
+
+A PID auto tuning utility is provided in the Pytec library. The auto tuning utility drives the the load to a controlled oscillation, observes the ultimate gain and oscillation period and calculates a set of PID parameters.
+
+To run the auto tuning utility, run
+
+```shell
+python pytec/autotune.py
+```
+
+After some time, the auto tuning utility will output the auto tuning results, below is a sample output
+
+```shell
+Ku: 0.7553203471147422
+Pu: 75.93899999999977
+rule: ziegler-nichols
+Kp: 0.45319220826884526
+Ki: 0.011935690706194357
+Kd: 4.301870387965967
+rule: tyreus-luyben
+Kp: 0.3432930977636503
+Ki: 0.0020549280832497956
+Kd: 4.137825730504864
+        .
+        .   
+        .
+```
+
+At the end of the test, the ultimate gain `Ku`, oscillation period `Pu` and a few sets of recommended PID parameters are calculated and displayed.  
+
+Multiple suggested sets of PID parameters based on different calculation rules are displayed. While all sets are expected to work, the different sets trade off response time with overshoot differently, and testing is needed to see which set works best for the system on hand.
+
+With a well designed and constructed setup, the PID parameters calculated by the auto tune utility can provide mK level control stability with no manual tweaking.
+
+![tec3pessen](/assets/tec3pessen.png)