diff --git a/src/MATLAB/ssf_vs_cn/goFigure3-C.m b/src/MATLAB/ssf_vs_cn/goFigure3-C.m
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/src/MATLAB/ssf_vs_cn/goPhantomImage.m b/src/MATLAB/ssf_vs_cn/goPhantomImage.m
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/src/MATLAB/ssf_vs_cn/goSSF_CN.m b/src/MATLAB/ssf_vs_cn/goSSF_CN.m
new file mode 100644
index 0000000000000000000000000000000000000000..90c47e34d5b8efa48c841d5b7e7f001645f5af61
--- /dev/null
+++ b/src/MATLAB/ssf_vs_cn/goSSF_CN.m
@@ -0,0 +1,48 @@
+
+sz = [100,100];
+cx = sz./2;
+imCN = zeros(sz);
+imSSF = zeros(sz);
+
+for i = 1 : sz(1)
+    for j = 1 : sz(2)
+
+        cyto = i./sz(1); 
+        nuc = j./sz(2);
+        
+        imCN(i,j) = cyto./nuc;
+        imSSF(i,j) = cyto - nuc;
+        
+    end
+end
+
+clf;imagesc(imCN)
+xt = get(gca,'xtick')
+set(gca,'xticklabels',xt./sz(2))
+yt = get(gca,'ytick')
+set(gca,'yticklabels',yt./sz(1))
+ylabel('cytoplasmic intensity')
+xlabel('nuclear intensity')
+
+% tick_targets = [1,64,128,192,255];
+cb = colorbar()
+cb.Label.String = 'cytonuclear ratio'
+
+exportgraphics(gcf,'imCN.png')
+
+clf;imagesc(imSSF)
+xt = get(gca,'xtick')
+set(gca,'xticklabels',xt./sz(2))
+yt = get(gca,'ytick')
+set(gca,'yticklabels',yt./sz(1))
+ylabel('cytoplasmic intensity')
+xlabel('nuclear intensity')
+
+tick_targets = [-0.99,-0.5,0,0.5,0.99];
+tick_labels = [-1,-0.5,0,0.5,1];
+cb = colorbar('ticks',tick_targets,'ticklabels',tick_labels)
+cb.Label.String = 'SSF'
+
+
+exportgraphics(gcf,'imSSF.png')
+