power.m

计算动力学系统的分岔图

function ans = power(x,y)
% X.^Y for adtayl objects.
[m,n,p1] = size(x.tc);
if (fix(y)==y)
     ans=powerin(x,y);    
elseif (x.tc(:,1)==0) && (p1-1<y) 
    for i=1:p1
        x.tc(:,:,i)=0;
    end
    ans=x;    
elseif (x.tc(:,1)==0) && (p1-1>y)
   error('X.^Y for adtayl objects is not defined'); 
else   
 [m,n,p1] = size(x.tc);
 y1=y*log(x);
 ans=exp(y1); 
end