diff --git a/jWavelet/morsewave.m b/jWavelet/morsewave.m
index 6ffa24f..fd14526 100644
--- a/jWavelet/morsewave.m
+++ b/jWavelet/morsewave.m
@@ -256,7 +256,10 @@
 X=vswap(X,inf,0);
 
 ommat=vrep(vrep(om,size(X,3),3),size(X,2),2);
-Xr=X.*rot(ommat.*(N+1)/2*fact); %ensures wavelets are centered 
+Xr=X.*(-1).^[0:N-1]'; %ensures wavelets are centered
+if ~mod(N, 2)
+  Xr(N/2+1, :)/=2; %ensures proper wavelet decay and analyticity
+endif
 %figure,plot(vrep(om,size(X,2),2).*(N+1)/2*fact)
 
 x=ifft(Xr);
@@ -267,7 +270,7 @@
 r=(2*be+1)./ga;
 c=r-1;
 L=0*om;
-index=(1:round(N/2));
+index=(1:floor(N/2)+1);
 psif=zeros(length(psizero),1,K);
 
 for k=0:K-1