diff --git a/demos/krylov_sampling.py b/demos/krylov_sampling.py
index ce7e9913743b96b93f6a46d9c8f4854cbafd4010..21237275e0612f604bef8a012ca2696b378250a7 100644
--- a/demos/krylov_sampling.py
+++ b/demos/krylov_sampling.py
@@ -83,8 +83,8 @@ plt.close()
 pltdict = {'alpha': .3, 'linewidth': .2}
 for i in range(N_samps):
     if i == 0:
-        plt.plot(sky(samps_old[i]).val, color='b', **pltdict, label='Traditional samples')
-        plt.plot(sky(samps[i]).val, color='r', **pltdict, label='Krylov samples')
+        plt.plot(sky(samps_old[i]).val, color='b', **pltdict, label='Traditional samples (residuals)')
+        plt.plot(sky(samps[i]).val, color='r', **pltdict, label='Krylov samples (residuals)')
     else:
         plt.plot(sky(samps_old[i]).val, color='b', **pltdict)
         plt.plot(sky(samps[i]).val, color='r', **pltdict)