FIX Options
Detection
 detection method of bifurcations : When the value of a system parameter varies, a situation that passes through a bifurcation value may occur. These are the options to control the program's behavior.
I  when detecting a bifurcation value, keep on the calculations without halting.  II  when detecting a bifurcation value, halt the program.  III  when detecting a bifurcation value, adjust the step size for a variational parameter and automatically compute the closest values when the bifurcation occurs. 
 SeigTol_val : If you selected the above detection method "III", the parameters are automatically adjusted based on the eigenvalues conditioning the saddlenode bifurcation. This is the parameter that will fix how close to the exact eigenvalue we want to set the precision.
 Saddlenode bifurcation of a equilibrium point
 If the distance between the current eigenvalue and the bifurcation condition () reached within the range of the SeigTol value, the program halts.
 Hopf bifurcation of a equilibriumpoint
 If the distance between the current eigenvalue and the bifurcation condition ( and ) reached within the range of the SeigTol value, the program halts.
 Saddlenode bifurcation of a limit cycle
 If the distance between the current characteristic multiplier (eigenvalue) and the bifurcation condition () reached within the range of the SeigTol value, the program halts.
 Perioddoubling bifurcation of a limit cycle
 If the distance between the current characteristic multiplier (eigenvalue) and the bifurcation condition () reached within the range of the SeigTol value, the program halts.
 NeimarkSacker bifurcation of a limit cycle
 If the distance between the current characteristic multiplier (eigenvalue) and the bifurcation condition () reached within the range of the SeigTol value, the program halts.
 Pitchfork bifurcation of a limit cycle
 If the distance between the current characteristic multiplier (eigenvalue) and the bifurcation condition () reached within the range of the SeigTol value, the program halts.
Graphic display
 Graphic display: The switch to display eigenvalues obtained from the characteristic equation.
 On : display the graph while computing
 Off : do not display the graph while computing
 single figure : see below.
 On : Display the results on one window for all computations.
 Off : Create a new window every computation. To be used for results comparison, for instance.
Newton's method
 parameters for Newton's methods： The precise position of an equilibrium and a fixed/periodic point of the limit cycle is computed by the Newton method. Each of the following items is a control parameter of this Newton method.
gmax  If the sum of all the absolute value of the derivatives with respect to the variables exceeds this value, the program considers it diverged.  emax  If the sum of all the absolute value of the system variables used in the Newton method exceeds this value, the program considers it diverged.  eps  If the sum of all the absolute value of the differential values of the system variables used in the Newton method goes under this value, the program considers it converged.  feps  If the sum of all the absolute value of the system variables used in the Newton method goes under this value, the program considers it converged.  iter  Limit amount of iterations for the Newton method to converge. If it exceeds this number, the program halts.  nnn  Extrapolating the curve already computed yields initial values greatly accelerating the Newton method's convergence speed. But such method cannot be used for the 2 first points of the curve. For the first and the second points, parameters must be partitioned. nnn controls the number of partitions. 
ODE solver
 parameters for ODE solver： Selection and detailed settings of numerical integration.
 Solver：selection of ODE solver built in MATLAB is possible.
ode45  This is based on an explicit RungeKutta formula, the DormandPrince pair. It is a onestep solver; that is, in computing y(tn), it needs only the solution at the immediately preceding time point, y(tn1). In general, ode45 is the best solver to apply as a "first try" for most problems.  ode23  This is also based on an explicit RungeKutta proposed by Bogacki and Shampine. It may be more efficient than ode45 at crude tolerances and in the presence of mild stiffness. ode23 is a onestep solver.  ode113  This is a variable order AdamsBashforthMoulton PECE solver. It may be more efficient than ode45 at stringent tolerances. ode113 is a multistep solver; that is, it normally needs the solutions at several preceding time points to compute the current solution.  ode23s};  This is based on a modified Rosenbrock formula of order 2. Because it is a onestep solver, it may be more efficient than ode15s at crude tolerances. It can solve some kinds of stiff problems for which ode15s is not effective.  ode15s  This is a variable order solver based on the numerical differentiation formulas (NDFs). These are related to but are more efficient than the backward differentiation formulas, BDFs (also known as Gear's method). Like ode113, ode15s is a multistep method solver. If you suspect that a problem is stiff or if ode45 failed or was very inefficient, try ode15s.  myself  This is a standard fixed step 4th order Runge Kutta method．It is not recommended unless under particular circumstances. 
 Initial step : The initial step size that the ODE solver uses. When using one of the preinstalled solvers, since the step size is automatically adjusted, setting an initial step size unreasonably small (enter a value of the input field too big) will result in poor computation efficiency. When selecting the myself method, its step size is simply 1/(the value of the input field).
 Relative tolerance : measures the error relative to the size of each state. The relative tolerance represents a percentage of the state's value.
 Absolute tolerance : represents the acceptable error as the value of the measured state approaches zero*1.
Symmetrical property
 Symmetrical property：Definition of symmetrical properties of limit cycles and the detailed settings
 Not use this option： Out of consideration of a symmetrical property of limit cycles, i.e., normal action of FIX.
 Inversion：The case of limit cycles with an inversion symmetry. Thus, the limit cycle have the following property:
 Permutation：The case of limit cycles with a permutation symmetry and a time shift. The limit cycle is an invariant under the operation by the following matrix and by the time shift:
,
where indicates the following matrix:
,
 Arbitrary：The case of limit cycle with an arbitrary symmetrical property. User must define the arbitrary symmetrical matrix.
 divided number：When the symmetry operation is carried out to a limit cycle, you must also take into account the time shift([1])．This option is given as a divided number that need to shift how long time for the period of the limit cycle.
References
 [1] Y. Katsuta and H. Kawakami, Bifurcations of equilibria and periodic solutions in a nonlinear autonomous system with symmetry, Electronics and Communications in Japan, Scripta Technica, Inc., Part 3, vol.76, no.7, pp.114, 1993.
