Introduction - If you have any usage issues, please Google them yourself
% Y = LorenzMap(X, p)
% Given a 3-D vector IC X=[x;y;z] and a parameter vector p=[sigma rho beta]
% LorenzMap returns the next incidence in which zdot=0 while zdot is decreasing
% (i.e., a local max of z).
% If the user does not supply p, the default values of the parameters
% p = [sigma, rho, beta] = [10, 28, 8/3] are used.