multilinear1ad.m

计算动力学系统的分岔图

function ytayl1 = multilinear1AD(func,q1,x0,p1,n) 

%----------------------------------------------------
% This file computes  A^(n)*q1 using automatic differentiation, where A^(n) is the Jacobian 
% of the map f^(n). 
%----------------------------------------------------

taylorder=1;
y1=Bvv(func,x0,q1,p1,taylorder,n);
ytayl1=tcs(y1);
%---------------------------------------------
function y1= Bvv(mapsf,x0,hc,p1,taylorder,n)
   
   % Convert to "active" variables:
   s = adtayl(0,taylorder);%Base point & Taylor order
   y1= x0 + s.*hc;
   for i=1:n     
     y1 = mapsf(0, y1,p1{:});
   end