diff --git a/vis_compare_model_size.m b/vis_compare_model_size.m
new file mode 100644
index 0000000..c5ec742
--- /dev/null
+++ b/vis_compare_model_size.m
@@ -0,0 +1,59 @@
+
+f_cp = @(D,N,I) N.*D.*I;
+f_tucker = @(D,N,I) D.^N + N.*D.*I;
+f_btd = @(D,N,I,B) B.*(D./B).^N + B.*N.*(D./B).*I;
+f_tt = @(D,N,I) D.*I.*2 + (N-2).*D.*I.*D;
+
+
+
+figure, 
+l_obs = [10,100,1000];
+l_modes = 3:6
+n_blocks = 2;
+rangeD = 1:50;
+y_max = f_tucker(rangeD(end),l_modes(end),l_obs(end));
+y_max = 10^ceil(log10(y_max));
+for i_obs = 1:length(l_obs)
+    for i_modes = 1:length(l_modes)
+
+    n_obs = l_obs(i_obs);
+    n_modes = l_modes(i_modes);
+    
+    
+    subplot(length(l_obs),length(l_modes), (i_obs-1)*length(l_modes) + i_modes)
+    hold on
+    lw=1.5;
+    plot(rangeD,f_cp(rangeD,n_modes,n_obs),':k','LineWidth',lw)
+	plot(rangeD,f_tucker(rangeD,n_modes,n_obs),'-k','LineWidth',lw)
+	plot(rangeD,f_btd(rangeD,n_modes,n_obs,n_blocks),'--k','LineWidth',lw)
+    plot(rangeD,f_tt(rangeD,n_modes,n_obs),'-.k','LineWidth',lw)
+    set(gca,'YScale','log','YTick',10.^(0:2:ceil(log10(y_max))))
+    ylim([0,y_max])
+    set(gca,'XTick',[min(rangeD),floor(median(rangeD)),max(rangeD)])
+    xlim([min(rangeD),max(rangeD)])
+%     semilogy(rangeD,f_cp(rangeD,n_modes,n_obs),...
+%                  rangeD,f_tucker(rangeD,n_modes,n_obs),...
+%                  rangeD,f_btd(rangeD,n_modes,n_obs,n_blocks),...
+%                  rangeD,f_tt(rangeD,n_modes,n_obs))
+    if i_obs == 1
+        title(sprintf('%i modes',l_modes(i_modes)))
+    end
+
+    if i_obs == length(l_obs)
+        xlabel('Components (D)')
+    else
+        set(gca,'XTick',[])
+    end
+    
+    if i_modes == 1
+        ylabel(sprintf('#Parameters\n (I=%i)',l_obs(i_obs)))
+    else
+        set(gca,'YTick',[])
+    end
+    
+%     set(gca,'Position',get(gca,'Position')-[0.05,0,0.01,0])
+    
+    
+    end
+end
+legend({'CP','Tucker','BTD','TT'},'FontSize',12)
\ No newline at end of file