p_secant.m

来自「这是国外用的研究分岔的完整的M程序」· M 代码 · 共 31 行

function secant=p_secant(secant,norm_point)

% function n_secant=p_secant(secant,norm_point)
% INPUT: 
%	secant unnormalized secant
%	norm_point norm of computed point
% OUTPUT:
%	n_secant normalized secant

% (c) DDE-BIFTOOL v. 2.00, 21/10/2001 

secant=p_axpy(norm_point/p_norm(secant),secant,[]);

if secant.kind=='psol' | secant.kind=='hcli'
  m=secant.degree;
  ll=(size(secant.profile,2)-1)/m;
  if isempty(secant.mesh)
    mesh=0:1/(ll*m):1;
  else
    mesh=secant.mesh;
  end;
  dmesh(1)=(mesh(2)-mesh(1))/2;
  dmesh(2:ll*m)=(mesh(3:ll*m+1)-mesh(1:ll*m-1))/2;
  dmesh(ll*m+1)=(mesh(ll*m+1)-mesh(ll*m))/2;
  for k=1:size(secant.profile,2)
    secant.profile(:,k)=secant.profile(:,k)*dmesh(k);
  end;
end;

return;