diff --git a/compute_mca.py b/compute_mca.py
new file mode 100644
index 0000000000000000000000000000000000000000..4807c071a442b6b8bbb41f2c237181af9c728308
--- /dev/null
+++ b/compute_mca.py
@@ -0,0 +1,24 @@
+import numpy as np
+from supersvd import supersvd
+import matplotlib.pyplot as plt
+
+ts = np.fromfile('ts.std', dtype=np.float32).reshape(1147, 28, 31)
+ps = np.fromfile('ps.std', dtype=np.float32).reshape(1147, 28, 31)
+
+svd = supersvd(ts, ps, 4)
+
+first_mode_ts = svd.x_vect[0]
+
+plt.matshow(first_mode_ts)
+plt.savefig('first_mode_ts.png')
+
+first_ps_timeseries = svd.x_coeff[0]
+
+plt.figure()
+plt.plot(first_ps_timeseries)
+plt.savefig('first_ps_timeseries.png')
+
+print('Correlation coefficient between first_ts_timeseries. and first_ps_timeseries:', svd.corrcoeff[0])
+
+
+