diff --git a/MATLAB/Algorithms/IRN_TV_CGLS.m b/MATLAB/Algorithms/IRN_TV_CGLS.m
index 47bf1e63..fbe75978 100644
--- a/MATLAB/Algorithms/IRN_TV_CGLS.m
+++ b/MATLAB/Algorithms/IRN_TV_CGLS.m
@@ -65,7 +65,7 @@
 end
 qualMeasOut=zeros(length(QualMeasOpts),niter);
 resL2 = zeros(1,niter);
-
+remember=0;
 iter=0;
 % for iii = 1:niter_outer
 while iter<niter
diff --git a/Python/demos/d09_Algorithms04.py b/Python/demos/d09_Algorithms04.py
index 5426ae1d..7a0918b5 100644
--- a/Python/demos/d09_Algorithms04.py
+++ b/Python/demos/d09_Algorithms04.py
@@ -189,9 +189,19 @@
     delta=np.array([-0.005]),
 )
 
+
+ 
+#  IRN-CGLS-TV: CGLS with TV regularization
+# ==========================================================================
+# ==========================================================================
+#
+# Should converge fairly fast to a TV regularized solution
+imgCGLSTV=algs.irn_tv_cgls(noise_projections,geo,angles,10,lmbda=5,niter_outer=2)
+
+
 #%% plot results
 
 # plot images
 tigre.plotimg(
-    np.concatenate([imgAWASDPOCS, imgOSASDPOCS, imgASDPOCS, imgOSSART], axis=1), dim="z", step=2
+    np.concatenate([imgAWASDPOCS, imgOSASDPOCS, imgASDPOCS, imgOSSART, imgCGLSTV], axis=1), dim="z", step=2
 )