tens2.m

计算动力学系统的分岔图

function T2=Tens2(mapsf,mapsH,x,p,n)
%
% Hessian of f at x1,x2, ...,xn

  global cds 
    x1=x;
    T2(:,:,:,1)= feval(mapsH, 0, x1, p{:});
    for m=2:n
     x1=feval(mapsf,0,x1,p{:});
     T2(:,:,:,m)= feval(mapsH, 0, x1, p{:});
    end