mpower.m

计算动力学系统的分岔图

function ans = mpower(x,y)
% X^Y for adtayl objects.

[m,n,p1] = size(x.tc);
if (fix(y)==y)
      ans=powerin(x,y);   %  ans=mpowerin(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