diff --git a/ecgsyn-1.0.0/Matlab/ecgsyn.m b/ecgsyn-1.0.0/Matlab/ecgsyn.m
index 22f1748..b24fa6a 100644
--- a/ecgsyn-1.0.0/Matlab/ecgsyn.m
+++ b/ecgsyn-1.0.0/Matlab/ecgsyn.m
@@ -181,7 +181,7 @@ dw2 = w-w2;
 Hw1 = sig1*exp(-0.5*(dw1/c1).^2)/sqrt(2*pi*c1^2);
 Hw2 = sig2*exp(-0.5*(dw2/c2).^2)/sqrt(2*pi*c2^2);
 Hw = Hw1 + Hw2;
-Hw0 = [Hw(1:n/2); Hw(n/2:-1:1)];
+Hw0 = [Hw(1:n/2+1); Hw(n/2:-1:2)];
 Sw = (sfrr/2)*sqrt(Hw0);
 
 ph0 = 2*pi*rand(n/2-1,1);