diff --git a/src/pynwb/misc.py b/src/pynwb/misc.py
index 14c2e08d1..9989b3b5f 100644
--- a/src/pynwb/misc.py
+++ b/src/pynwb/misc.py
@@ -137,11 +137,42 @@ class Units(DynamicTable):
         'resolution'
     )
 
-    waveforms_desc = ('Individual waveforms for each spike. If the dataset is three-dimensional, the third dimension '
-                      'shows the response from different electrodes that all observe this unit simultaneously. In this'
-                      ' case, the `electrodes` column of this Units table should be used to indicate which electrodes '
-                      'are associated with this unit, and the electrodes dimension here should be in the same order as'
-                      ' the electrodes referenced in the `electrodes` column of this table.')
+    __waveforms_desc = (
+        """
+        Individual waveforms for each spike. If the dataset is three-dimensional, the third dimension
+        shows the response from different electrodes that all observe this unit simultaneously. In this
+        case, the ``electrodes`` column of this Units table should be used to indicate which electrodes
+        are associated with this unit, and the electrodes dimension here should be in the same order as
+        the electrodes referenced in the ``electrodes`` column of this table.
+
+        Example usage::
+
+            waveforms_list = [
+                [  # unit 1
+                    np.array([  # electrode 1
+                        [1, 2, 3, 4, 5],  # spike time 1 [sample 1, sample 2, ...]
+                        [2, 3, 4, 5, 6],  # spike time 2
+                    ]),
+                    np.array([  # electrode 2
+                        [3, 4, 5, 6, 7],  # spike time 1 [sample 1, sample 2, ...]
+                    ]),
+                ],
+                [  # unit 2
+                    np.array([  # electrode 1
+                        [10, 20, 30, 40, 50],  # spike time 1 [sample 1, sample 2, ...]
+                    ]),
+                ],
+            ]
+            electrodes_list = [[1, 2], [3]]
+            for unit_id in range(2):
+                nwbfile.add_unit(
+                    id=unit_id,
+                    electrodes=electrodes_list[unit_id],
+                    waveforms=waveforms_list[unit_id]
+                )
+        """
+    )
+
     __columns__ = (
         {'name': 'spike_times', 'description': 'the spike times for each unit', 'index': True},
         {'name': 'obs_intervals', 'description': 'the observation intervals for each unit',
@@ -151,7 +182,7 @@ class Units(DynamicTable):
         {'name': 'electrode_group', 'description': 'the electrode group that each spike unit came from'},
         {'name': 'waveform_mean', 'description': 'the spike waveform mean for each spike unit'},
         {'name': 'waveform_sd', 'description': 'the spike waveform standard deviation for each spike unit'},
-        {'name': 'waveforms', 'description': waveforms_desc, 'index': 2}
+        {'name': 'waveforms', 'description': __waveforms_desc, 'index': 2}
     )
 
     @docval({'name': 'name', 'type': str, 'doc': 'Name of this Units interface', 'default': 'Units'},
@@ -195,7 +226,7 @@ def __init__(self, **kwargs):
              'default': None},
             {'name': 'waveform_sd', 'type': 'array_data', 'default': None,
              'doc': 'the spike waveform standard deviation for each unit. Shape is (time,) or (time, electrodes)'},
-            {'name': 'waveforms', 'type': 'array_data', 'default': None, 'doc': waveforms_desc,
+            {'name': 'waveforms', 'type': 'array_data', 'default': None, 'doc': __waveforms_desc,
              'shape': ((None, None), (None, None, None))},
             {'name': 'id', 'type': int, 'default': None, 'doc': 'the id for each unit'},
             allow_extra=True)