multilinear2ad.asv

计算动力学系统的分岔图

function ytayl2 = multilinear2AD(func,q1,q2,x0,p1,n) 

%----------------------------------------------------
% This file computes  B^(n)*q1*q2 using automatic differentiation, where 
%B^(n) is the second order derivatives  of map f^(n). 
%----------------------------------------------------
 taylorder=2; 
 if q1==q2
     
  y1=Bvv(func,x0,q1,p1,taylorder,n);
 else
    
    y11=Bvv(func,x0,q1+q2,p1,taylorder,n),pause
    
    y12=Bvv(func,x0,q1-q2,p1,taylorder,n),pause
    
    y1=1/4.0*(y11-y12),pause
 end
y1
ytayl1=tcs(y1);
%dcoeffs = y1.tc(:,:);
%ytayl2=2*dcoeffs(:,taylorder+1);
%ytayl2=2*y1(:,end),pause
%--------------------------------------------------
function tempvec = Bvv(mapsf,x0,hc,p1,taylorder,n)
   % Convert to "active" variables:
   s = adtayl(0,taylorder);%Base point & Taylor order
   y0=adtayl(x0);
   h=adtayl(hc);
   y1= y0 + s*h;
 for i=1:n       
     y1 = mapsf(0, y1,p1{:});
 end
tempvec=y1;