diff --git a/oldMcode/compare.m b/oldMcode/compare.m
new file mode 100644
index 0000000..3383e9b
--- /dev/null
+++ b/oldMcode/compare.m
@@ -0,0 +1,26 @@
+maxSize=6;
+
+mseErrorAutoenc=zeros([maxSize,1]);
+for i=1:maxSize
+    ae=trainAutoencoder(trainingData,i,'ShowProgressWindow',false);
+    predicted=predict(ae,trainingData);
+    mseErrorAutoenc(i)=mse(trainingData-predicted);
+end
+
+fig=figure();
+subplot(1,3,1);
+plot(mseErrorAutoenc)
+
+%PCA
+[COEFF,SCORE,latent] = princomp(trainingData);
+subplot(1,3,2);
+plot(cumsum(latent)./sum(latent))
+
+%NNMF
+D=cell(maxSize,1);
+for i=1:maxSize
+    [W,H,D{i}]=nnmf(trainingData,i);
+end
+
+subplot(1,3,3);
+plot(diff(cell2mat(D)))
\ No newline at end of file
diff --git a/oldMcode/determineSizeOfHiddenLayer.m b/oldMcode/determineSizeOfHiddenLayer.m
deleted file mode 100644
index f90fc51..0000000
--- a/oldMcode/determineSizeOfHiddenLayer.m
+++ /dev/null
@@ -1,8 +0,0 @@
-maxSize=10;
-
-mseError=zeros([maxSize,1]);
-for i=1:maxSize
-    ae=trainAutoencoder(trainingData,i,'ShowProgressWindow',false);
-    predicted=predict(ae,trainingData);
-    mseError(i)=mse(trainingData-predicted);
-end
\ No newline at end of file
diff --git a/papers/2006-tresch-nnmf-for-synergies.pdf b/papers/2006-tresch-nnmf-for-synergies.pdf
new file mode 100644
index 0000000..5c699c7
--- /dev/null
+++ b/papers/2006-tresch-nnmf-for-synergies.pdf
Binary files differ