diff --git a/epod.m b/epod.m
index 5bfee460b3dde79e1df1a4282cab11087f1bf5ed..1c6d01585abdcfb77b639ad31a1213f2038ff94d 100644
--- a/epod.m
+++ b/epod.m
@@ -61,8 +61,8 @@ P_1_2 = a1*Sig_inv1*g_2;
 P_2_1 = a2*Sig_inv2*g_1;
  
 % Obtain extended singular values
-Sig_1_2 = sqrt(diag(P_1_2*P_1_2'));
-Sig_2_1 = sqrt(diag(P_2_1*P_2_1'));
+Sig_1_2 = diag(sqrt(diag(P_1_2*P_1_2')));
+Sig_2_1 = diag(sqrt(diag(P_2_1*P_2_1')));
  
 % Obtain psuedo-inverse
 Sig_inv_1_2 = diag(1./diag(Sig_1_2));