function  f=DESA(x)
%能量分離算法求解|a(n)|
% 算法詳細描述： open TheoryofDEO.doc
%                     2Ψ[x(n)]
%|a(n)|=-------------------------------------
%         sqrt(  Ψ[ x(n+1)-x(n-1) ]  )
U=2*DEO3S(x);
%L=abs(   real(  sqrt(   DEO3S(  circshift(x',-1)'-circshift(x',1)'  )   )  )   );
L=abs(   (  sqrt(   DEO3S(  circshift(x',-1)'-circshift(x',1)'  )   )  )   );

f=abs(U./L);%應去除分母為0的點
A=find(f==inf);
f(A)=0;