Code covered by the BSD License

### Highlights fromDynamical Systems Toolbox

from Dynamical Systems Toolbox by Etienne Coetzee
Bifurcation analysis of dynamical systems. Integration of AUTO bifurcation software into MATLAB.

writeinterface.m
```Ndim=int8(34);      % Problem dimension
Ips=1;              % Problem type
Irs=0;              % Start solution label
Ilp=1;              % Fold detection
Nicp=1;             % Continuation parameters
Icp=[1];            % Array of continuation parameters
Ntst=10;            % Number of mesh intervals
Ncol=4;             % Number of collocation points
Isp=1;              % Bifurcation detection
Isw=0;              % Branch switching
Iplt=0;             % Select principal solution measure
Nbc=0;              % Number of boundary conditions
Nint=0;             % Number of integral conditions
Nmx=1000;           % Maximum number of steps
Rl0=-3.1416;        % Parameter interval RL0 = \$\$\lamda\$\$ = RL1
Rl1=3.1416;         %
A0=-1e6;            % Interval of principal solution measure A0 = |||| = A1
A1=1e6;             %
Npr=10;	           % Print and save restart data every NPR steps
Mxbf=10;            % Automatic branch switching for the fit MXBF bifurcation points
Iid=2;              % Control diagnostic output
Itmx=8;	           % Maximum number of iterations for locating special solutions/points
Itnw=5;             % Maximum number of correction steps
Nwtn=3;             % Corrector uses full Newton for NWTN steps
Jac=0;              % User defines derivatives
Epsl=1e-9;          % Convergence criterion for equation parameters
Epsu=1e-9;          % Convergence criterion for solution components
Epss=1e-8;          % Convergence criterion for special solutions/points
Ds=0.0175;          % Start step size
Dsmin=0.00175;      % Step size interval DSMIN =< DS =< DSMAX
Dsmax=0.175;        %
Nthl=[0];           % List of parameter and solution weights
Thl=[];
Nthu=[0];	       %
Thu=[];
Nuzr=[0];           % List of values for user defined output
Uzr=[];

C Assign NDIM values
NDIMptr1=mxGetProperty(PLHS(1),i,'Ndim')
C
IF( NDIMptr1 == 0 ) THEN
ENDIF
C
IF( mxIsInt8(NDIMptr1) == 0 ) THEN
CALL mexErrMsgTxt("Property 'NDIM' is not an integer")
ENDIF
C
NDIMptr2=mxGetPr(NDIMptr1)
CALL mxCopyPtrToInteger4(NDIMptr2,NDIMDST,n)
CALL mxDestroyArray(NDIMptr1)
C       ```