main_leuc.m

% the stable manifold of the saddle point in yoz

[u0,A,f,kont]=jacob([0;1000;2000]);% the critical point and the Jacobian

f,kont,
[Ps,Pu,L,V]=proj(A);% the eigenvalues and the eigenvectors of the Jacobian,
                    % the projection operators onto the stable eigenspace
work2;