tens1.m

计算动力学系统的分岔图

function T1=Tens1(mapsf,mapsJ,x,p,n)
% Jacobian of f at x1,x2,....,xN    
  global cds 
    x1=x;
    T1(:,:,1)= feval(mapsJ, 0, x1, p{:});
    for m=2:n
     x1=feval(mapsf,0,x1,p{:});
     T1(:,:,m)= feval(mapsJ, 0, x1, p{:});
    end