function Lyapunov(x)
if length(x)<400
   helpdlg('數據量需要大于400，請重新選定時間','提示信息');pause(3),
end   
dm=26;tao=5;set_error=0.001;
for i=0:(dm-1)
    for j=1:(length(x)/1.2)
       k=j+i*tao;
           Y(i+1,j)=x(k);
    end
end 

l=1; i=1; p=1;                      %  i為Y的一列，j為Y的一列 
while i<length(Y)
     for j=p:length(Y)               %求兩點之間距離
        if j>i
        e_normal=max(abs(Y(:,i)-Y(:,j)));
        end
     end    
    flag=e_normal-set_error;
            if flag>0
              d(l)=e_normal;l=l+1;
              %i=i+1;
                 for j=1:length(Y)             %找離點最小距離的點  
                    l_normal(j)=max(abs(Y(:,i)-Y(:,j)));
                 end
              p=find((l_normal-min(l_normal))==0);            %找到最近點Y(:,j)
           else 
              j=j+1;i=i+1;                    %范數小于set_error，Y(:,i)和Y(:,j)繼續向前一步  
           end
end
 
LE1=0;
for i=1:l
    f1=d(i);f2=d(1);
 LE1=log10(f1/f2)+LE1;
end
LE1=1/(l*0.014)*LE1,