sqrt.m

计算动力学系统的分岔图

function ans = sqrt(x)

% sqrt for adtayl objects.

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