diff --git a/BayesValidRox/PostProcessing/PostProcessing.py b/BayesValidRox/PostProcessing/PostProcessing.py
index a0a0ca30a4cfcde2241df86782eb9942ad727f7b..36c0baf3a1632f8756327c722e67d3af3f34fac4 100644
--- a/BayesValidRox/PostProcessing/PostProcessing.py
+++ b/BayesValidRox/PostProcessing/PostProcessing.py
@@ -745,9 +745,22 @@ class PostProcessing:
                     else:
                         plotLabel = plot
                         SeqValues = np.nan_to_num(SeqValues)
+                        
+                        # Plot the error evolution for each output
+                        for i in range(SeqValues.shape[1]):
+                            plt.semilogy(x_idx, SeqValues.mean(axis=1), marker=markers[idx],
+                                      ls='--', lw=2, color=colors[idx], alpha=0.25)
+                        
+                        # plt.fill_between(x_idx, SeqValues.mean(axis=1)-2*np.std(SeqValues,axis=1), 
+                        #                  SeqValues.mean(axis=1)+2*np.std(SeqValues,axis=1), color=colors[idx],alpha=0.15)
+                
+                        
                         plt.semilogy(x_idx, SeqValues.mean(axis=1), marker=markers[idx],
                                      ls='--', lw=2, color=colors[idx], label=name.split("_rep",1)[0])
-                
+                        
+                        
+                        
+                        
                 plt.xticks(x_idx)
                 plt.xlabel('Number of runs')
                 plt.ylabel(plotLabel)