tens2.m

来自「计算动力学系统的分岔图」· M 代码 · 共 13 行

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