generate simulation code based on a MapleSim/Multibody model module - MapleSim Help

 Calling Sequence BuildSimCode(MSimMultibodyModel, options)
 Parameters MSimMultibodyModel - Maple module of a system created using the BuildEQs command or retrieved from a MapleSim/Multibody store file (.lib) using the GetModel command. options - (optional) equations of the form keyword = value
 Description Important: The BuildSimCode command has been deprecated. Use the templates in the MapleSim templates folder to analyze your model instead. For more information, see Opening MapleSim Templates. The BuildSimCode command can be used to generate optimized procedures and functions in a variety of output languages. Unless the SilentMode option is set to true, BuildSimCode indicates the location and argument syntax for any generated functions. Matrix outputs are returned as column matrices in the column-major format.

Options

'TargetDir' = string

A string indicating the directory to which the BuildSimCode command exports the generated simulation code. If the string is assigned to a directory that does not exist or not assigned to a directory, BuildSimCode exports the simulation code to the current directory as defined by the currentdir command.  This option is ignored if the 'Language' option is set to MapleProc or MapleFunc.

'Exports' = set

The 'Exports' option indicates the functions for which BuildSimCode exports code. Two different syntaxes can be used depending on whether you want to use pre-defined or user-defined functions. The arguments of the generated functions are dependent on the selected 'Language' option.

- Pre-defined Functions: in this case, 'Exports'  is equated to a set of strings. Each string indicates a particular pre-defined function (that is, 'Exports' = {"Xdot", "Reactions"}). For more information about pre-defined functions, see the kinematics and dynamics help pages. The allowable strings are listed in the table below:

 Allowable Strings Description "Xdot" The first time derivatives of state variables, that is, $\frac{ⅆp}{ⅆ\phantom{\rule[-0.0ex]{0.2em}{0.0ex}}t}$ and $\frac{ⅆq}{ⅆ\phantom{\rule[-0.0ex]{0.2em}{0.0ex}}t}$ . "Reactions" The m reaction loads corresponding to either the vReactions or vLambda column matrix of the model, depending on the 'AugType' option used in BuildEQs when the model was generated. "PosCons" The m position-level constraint equations corresponding to the GetPosCons() function of the model. "VelCons" The m velocity-level constraint equations corresponding to the GetVelCons() function of the model. "JacPos" The m x n Jacobian matrix of the position constraint equations with respect to the generalized coordinates, corresponding to the model's xJacPos export.  The matrix is exported as an array using column-major ordering. "JacVel" The m x n Jacobian matrix of the velocity constraint equations with respect to the generalized speeds, corresponding to the model's xJacVel export.  The matrix is exported as an array using column-major ordering. "VelRHS" Right side of the velocity-level kinematic constraint equations, corresponding to the vVelRHS column matrix of the model. "AccRHS" Right side of the acceleration-level kinematic constraint equations, corresponding to the vAccRHS column matrix of the model. "C" The m x n coefficient matrix of constraint reactions, corresponding to the xC export of the model. "M" The n x n mass matrix, corresponding to the xM export of the model. "F" The n x 1 right side of the dynamic equations, corresponding to the vF column matrix of the model. "A" The first matrix returned by the export GetAYBSysODEs(). For more information, see the Dynamic Exports help page. "B" The third matrix returned by the export GetAYBSysODEs(). For more information, see the Dynamic Exports help page.

 • User-defined Vectors and Matrices:  if a user-defined Matrix or Vector has been created in a Maple worksheet, you can export it to any of the languages supported by BuildSimCode. You can use a list of lists instead of using strings.

'Exports' = {[["Func1Name", Func1VectorOrMatrix, Func1Subs], ["Func2Name", Func2VectorOrMatrix, Func2Subs], ...]}

The first element in each sublist is a string indicating the name of the function to be generated. The second element is a Matrix or Vector object that has been created in a Maple worksheet. The final element in each sublist is a list of equations defining symbols that are not known inputs, parameters, or state variables. In most cases, this entry is left as an empty list (that is, [ ]).

'Language'= string

A string indicating the format of the exported simulation code. Allowable strings are listed in the table below:

 Allowable Strings Description "MapleProc" Generates a Maple procedure of the form PROCNAME(NumStates, t, vX, Output) This procedure is designed to be integrated using dsolve/numeric.   Tip: Use this option only when generating a procedure for the derivative of the model's state variables. All of the model parameters and inputs must be replaced by numeric values or functions of time to use this option. "MapleFunc" Creates a Maple function of the form  Output = FUNCNAME(t,vX,Params,Input)  where t is time, vX is the column matrix of state variables q and p, Params and Inputs are the parameters and inputs required to evaluate the procedure (use the model functions GetParams() and GetInputs() commands to identify them) and Output is the computed output. "Cproc" Creates a C-procedure of the form  PROCNAME(t,vX,Params,Inputs,Output) where t is time, vX is the column matrix of state variables q and p, Params and Inputs are the parameters and inputs required to evaluate the procedure (use the model functions GetParams() and GetInputs() commands to identify them) and Output is the computed output. "BB_SFunc" This option applies when the MapleSim Connector is installed. When selected, BuildSimCode creates an object in the Maple worksheet and saves the MATLAB® script required to automatically generate a Simulink® S-Function of the chosen procedure to the designated 'TargetDir'. The resulting Simulink® block will be an input or output block with no internal states. To facilitate this process, the state variables vX are moved from the State Variable Array to the Inputs Array. "BB_SFunc_WithIntegrator" This option applies when the MapleSim Connector is installed. When selected, BuildSimCode creates an object in the Maple worksheet and saves the MATLAB® script required to automatically generate a Simulink® S-Function of the chosen procedure to the designated 'TargetDir'.   The resulting Simulink® block will contain internal states (corresponding to the model's vX column matrix) and an integrator that will advance the state variables in time.

'Optimize'= boolean

A Boolean value indicating whether to apply optimization routines to the simulation code before the final code generation.  The default value is true.

'NumericSubs'= list

A list of equations of the form parameter=value. Any parameters or functions that are not included in this list are left as inputs to the generated simulation code.

'ICs'= list

A list of floats indicating the initial conditions for subsequent integration.  The ordering of this list is assumed to correspond to the order of the state variables in the model's vX export. This option is ignored unless the 'Language' option is set to "BB_SFunction_WithIntegrator".

'SilentMode'= boolean

Setting this option to true hides the output generated by the BuildSimCode command.  The default value is true.

Examples

Important: The BuildSimCode command has been deprecated. Use the templates in the MapleSim document folder to analyze your model instead.

The example below shows an embedded MapleSim model of a planar slider-crank mechanism. In this example, the embedded component name is Simulation0 and the model name is SliderCrank. When the code is generated, you can use the BuildSimCode command to export a Maple procedure to a worksheet.

 > $\mathrm{leArguments}≔\left['\mathrm{InputString}'=\mathrm{MapleSim}:-\mathrm{Multibody}:-\mathrm{GetMultibodyInput}\left("Simulation0"\right),'\mathrm{ModelName}'="SliderCrank",'\mathrm{KinSimpType}'="Simplify",'\mathrm{DynSimpType}'="Simplify",'\mathrm{AugType}'="Lagrange",'\mathrm{MaxSmallQOrder}'=1,'\mathrm{SaveToLib}'=\mathrm{true},'\mathrm{SilentMode}'=\mathrm{false}\right]:$$\mathrm{MapleSim}:-\mathrm{Multibody}:-\mathrm{BuildEQs}\left(\mathrm{op}\left(\mathrm{leArguments}\right)\right):$
 ${"Analyzing system..."}$
 ${"Performing constraint analysis..."}$
 ${"The system has 1 degree\left(s\right) of freedom. It is modeled using 3 generalized coordinate\left(s\right) coupled by 2 algebraic constraint\left(s\right)."}$
 ${"Peforming a dynamic analysis using an augmented Lagrange formulation - system variables shown below:"}$
 ${"vQ"}{,}\left[\begin{array}{c}{\mathrm{s_P2}}{}\left({t}\right)\\ {\mathrm{theta_R1}}{}\left({t}\right)\\ {\mathrm{theta_R5}}{}\left({t}\right)\end{array}\right]{,}{"vP"}{,}\left[\begin{array}{c}\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{s_P2}}{}\left({t}\right)\\ \frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{theta_R1}}{}\left({t}\right)\\ \frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{theta_R5}}{}\left({t}\right)\end{array}\right]{,}{"vPdot"}{,}\left[\begin{array}{c}\frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{\mathrm{s_P2}}{}\left({t}\right)\\ \frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{\mathrm{theta_R1}}{}\left({t}\right)\\ \frac{{{ⅆ}}^{{2}}}{{ⅆ}{{t}}^{{2}}}{}{\mathrm{theta_R5}}{}\left({t}\right)\end{array}\right]{,}{"vLambda"}{,}\left[\begin{array}{c}{\mathrm{λ1}}{}\left({t}\right)\\ {\mathrm{λ2}}{}\left({t}\right)\end{array}\right]$
 ${"Dynamic analysis complete."}$
 ${"Saving model for future use...:"}$
 ${"Model saved to: C:/SliderCrank.lib"}$
 ${"Use Model := MapleSim:-Multibody:-GetModel\left("C:/SliderCrank.lib"\right): to retrieve stored model."}$ (5.1)
 > $\mathrm{leNumericSubs}≔\left[G=9.81,\mathrm{m_1}=1,\mathrm{m_2}=1,\mathrm{m_3}=3,F\left(t\right)=0\right]:$
 > $\mathrm{leArguments}≔\left[\mathrm{SliderCrank},'\mathrm{Language}'="MapleProc",'\mathrm{Optimize}'=\mathrm{true},'\mathrm{NumericSubs}'=\mathrm{leNumericSubs},'\mathrm{Baumgarte}'=\left[2,2\right],'\mathrm{Exports}'=\left\{"Xdot"\right\},'\mathrm{SilentMode}'=\mathrm{false}\right]:$$\mathrm{MapleSim}:-\mathrm{Multibody}:-\mathrm{BuildSimCode}\left(\mathrm{op}\left(\mathrm{leArguments}\right)\right):$
 ${"Order of the State Variable Array:"}$
 $\left[\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{s_P2}}{}\left({t}\right){,}\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{theta_R1}}{}\left({t}\right){,}\frac{{ⅆ}}{{ⅆ}{t}}{}{\mathrm{theta_R5}}{}\left({t}\right){,}{\mathrm{s_P2}}{}\left({t}\right){,}{\mathrm{theta_R1}}{}\left({t}\right){,}{\mathrm{theta_R5}}{}\left({t}\right)\right]$
 ${"Order of the Constant Parameter Array:"}$
 $\left[{}\right]$
 ${"Order of the Inputs Array:"}$
 $\left[{}\right]$
 ${"--------------Processing Xdot-------------"}$
 ${"Procedure pXdot created in the Maple workspace."}$
 ${"pXdot\left(NumStates, t, States, StateDerivatives\right) #compatible with dsolve/numeric"}$
 ${4}{}{\mathrm{functions}}{+}{25}{}{\mathrm{subscripts}}{+}{96}{}{\mathrm{assignments}}{+}{83}{}{\mathrm{multiplications}}{+}{60}{}{\mathrm{additions}}{+}{3}{}{\mathrm{divisions}}$
 ${"----duration: .63e-1 seconds----"}$
 ${""}$ (5.2)
 > $\mathrm{eval}\left(\mathrm{pXdot}\right)$
 ${\mathbf{proc}}\left({N}{,}{t}{,}{\mathrm{vX_in}}{,}{\mathrm{XDOT}}\right)\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathbf{local}}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{:=}{\mathrm{array}}{}\left({1}{..}{59}\right){;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{1}{]}{:=}{\mathrm{sin}}{}\left({\mathrm{vX_in}}{[}{6}{]}\right){;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}{0.2500000000}{*}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{38}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{3}{]}{:=}{\mathrm{cos}}{}\left({\mathrm{vX_in}}{[}{6}{]}\right){;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{39}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{40}{]}{:=}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{4}{]}{:=}{\mathrm{sin}}{}\left({\mathrm{vX_in}}{[}{5}{]}\right){;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{5}{]}{:=}{\mathrm{cos}}{}\left({\mathrm{vX_in}}{[}{5}{]}\right){;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{6}{]}{:=}{−}{0.5000000000}{*}{\mathrm{zz}}{[}{3}{]}{*}{\mathrm{zz}}{[}{4}{]}{+}{0.5000000000}{*}{\mathrm{zz}}{[}{1}{]}{*}{\mathrm{zz}}{[}{5}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{41}{]}{:=}{\mathrm{zz}}{[}{6}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{7}{]}{:=}{0.5000000000}{*}{\mathrm{zz}}{[}{1}{]}{*}{\mathrm{zz}}{[}{4}{]}{+}{0.5000000000}{*}{\mathrm{zz}}{[}{3}{]}{*}{\mathrm{zz}}{[}{5}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{42}{]}{:=}{\mathrm{zz}}{[}{7}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{43}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{6}{]}{+}{\mathrm{zz}}{[}{1}{]}{*}{\mathrm{vX_in}}{[}{4}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{44}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{8}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{7}{]}{+}{\mathrm{zz}}{[}{3}{]}{*}{\mathrm{vX_in}}{[}{4}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{45}{]}{:=}{\mathrm{zz}}{[}{8}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{46}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{47}{]}{:=}{\mathrm{zz}}{[}{6}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{48}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{49}{]}{:=}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{50}{]}{:=}{\mathrm{zz}}{[}{7}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{51}{]}{:=}{\mathrm{zz}}{[}{8}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}{\mathrm{vX_in}}{[}{3}{]}{^}{2}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{52}{]}{:=}{−}{0.2500000000}{*}{\mathrm{zz}}{[}{2}{]}{*}{\mathrm{zz}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{53}{]}{:=}{−}{0.2943000000}{*}{\mathrm{zz}}{[}{5}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{54}{]}{:=}{2.452500000}{*}{\mathrm{zz}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{55}{]}{:=}{\mathrm{vX_in}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{56}{]}{:=}{\mathrm{vX_in}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{57}{]}{:=}{\mathrm{vX_in}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{6}{]}{:=}{2.}{*}{\mathrm{zz}}{[}{5}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{7}{]}{:=}{0.5000000000}{*}{\mathrm{vX_in}}{[}{2}{]}{^}{2}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{8}{]}{:=}{2.}{*}{\mathrm{zz}}{[}{4}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{4}{]}{:=}{\mathrm{zz}}{[}{7}{]}{*}{\mathrm{zz}}{[}{4}{]}{-}{1.}{*}{\mathrm{zz}}{[}{6}{]}{*}{\mathrm{vX_in}}{[}{2}{]}{-}{0.5000000000}{*}{\mathrm{zz}}{[}{2}{]}{*}{\mathrm{zz}}{[}{4}{]}{+}\left({\mathrm{zz}}{[}{6}{]}{-}{4.}{*}{\mathrm{vX_in}}{[}{4}{]}\right){*}{\mathrm{vX_in}}{[}{3}{]}{-}{1.}{*}{\mathrm{zz}}{[}{8}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{9}{]}{:=}{4.}{*}{\mathrm{vX_in}}{[}{4}{]}{+}{4.}{*}{\mathrm{vX_in}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}\left({\mathrm{vX_in}}{[}{4}{]}{-}{0.5000000000}{*}{\mathrm{zz}}{[}{5}{]}\right){*}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{5}{]}{:=}{\mathrm{zz}}{[}{7}{]}{*}{\mathrm{zz}}{[}{5}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{7}{]}{:=}{\mathrm{zz}}{[}{8}{]}{*}{\mathrm{vX_in}}{[}{3}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{8}{]}{:=}{\mathrm{zz}}{[}{8}{]}{*}{\mathrm{vX_in}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{58}{]}{:=}{\mathrm{zz}}{[}{1}{]}{*}{\mathrm{zz}}{[}{4}{]}{+}\left({\mathrm{zz}}{[}{5}{]}{+}{\mathrm{zz}}{[}{8}{]}{+}{\mathrm{zz}}{[}{2}{]}{-}{1.}{*}{\mathrm{zz}}{[}{7}{]}{+}{\mathrm{zz}}{[}{9}{]}{-}{1.}{*}{\mathrm{zz}}{[}{6}{]}\right){*}{\mathrm{zz}}{[}{3}{]}{-}{4.}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{59}{]}{:=}\left({−}{1.}{*}{\mathrm{zz}}{[}{5}{]}{-}{1.}{*}{\mathrm{zz}}{[}{8}{]}{-}{1.}{*}{\mathrm{zz}}{[}{2}{]}{+}{\mathrm{zz}}{[}{7}{]}{-}{1.}{*}{\mathrm{zz}}{[}{9}{]}{+}{\mathrm{zz}}{[}{6}{]}\right){*}{\mathrm{zz}}{[}{1}{]}{+}{\mathrm{zz}}{[}{3}{]}{*}{\mathrm{zz}}{[}{4}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{21}{]}{:=}{0.4000000000}{*}{\mathrm{zz}}{[}{38}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{22}{]}{:=}{0.4000000000}{*}{\mathrm{zz}}{[}{39}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{23}{]}{:=}{0.4000000000}{*}{\mathrm{zz}}{[}{40}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{26}{]}{:=}{0.1290000000}{-}{1.}{*}{\mathrm{zz}}{[}{43}{]}{*}{\mathrm{zz}}{[}{21}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{27}{]}{:=}{\mathrm{zz}}{[}{44}{]}{-}{1.}{*}{\mathrm{zz}}{[}{43}{]}{*}{\mathrm{zz}}{[}{22}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{28}{]}{:=}{\mathrm{zz}}{[}{45}{]}{-}{1.}{*}{\mathrm{zz}}{[}{43}{]}{*}{\mathrm{zz}}{[}{23}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{29}{]}{:=}{\mathrm{zz}}{[}{48}{]}{-}{1.}{*}{\mathrm{zz}}{[}{46}{]}{*}{\mathrm{zz}}{[}{21}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{30}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{46}{]}{*}{\mathrm{zz}}{[}{22}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{31}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{46}{]}{*}{\mathrm{zz}}{[}{23}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{32}{]}{:=}{\mathrm{zz}}{[}{51}{]}{-}{1.}{*}{\mathrm{zz}}{[}{49}{]}{*}{\mathrm{zz}}{[}{21}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{33}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{49}{]}{*}{\mathrm{zz}}{[}{22}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{34}{]}{:=}{−}{1.}{*}{\mathrm{zz}}{[}{49}{]}{*}{\mathrm{zz}}{[}{23}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{24}{]}{:=}{133.3333333}{*}{\mathrm{zz}}{[}{41}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{25}{]}{:=}{133.3333333}{*}{\mathrm{zz}}{[}{42}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{30}{]}{:=}{\mathrm{zz}}{[}{30}{]}{-}{1.}{*}{\mathrm{zz}}{[}{47}{]}{*}{\mathrm{zz}}{[}{24}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{31}{]}{:=}{\mathrm{zz}}{[}{31}{]}{-}{1.}{*}{\mathrm{zz}}{[}{47}{]}{*}{\mathrm{zz}}{[}{25}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{33}{]}{:=}{\mathrm{zz}}{[}{33}{]}{-}{1.}{*}{\mathrm{zz}}{[}{50}{]}{*}{\mathrm{zz}}{[}{24}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{34}{]}{:=}{\mathrm{zz}}{[}{34}{]}{-}{1.}{*}{\mathrm{zz}}{[}{50}{]}{*}{\mathrm{zz}}{[}{25}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{1}{]}{:=}{1}{/}{\mathrm{zz}}{[}{26}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{27}{]}{:=}{\mathrm{zz}}{[}{27}{]}{*}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{28}{]}{:=}{\mathrm{zz}}{[}{28}{]}{*}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{30}{]}{:=}{\mathrm{zz}}{[}{30}{]}{-}{1.}{*}{\mathrm{zz}}{[}{29}{]}{*}{\mathrm{zz}}{[}{27}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{31}{]}{:=}{\mathrm{zz}}{[}{31}{]}{-}{1.}{*}{\mathrm{zz}}{[}{29}{]}{*}{\mathrm{zz}}{[}{28}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{33}{]}{:=}{\mathrm{zz}}{[}{33}{]}{-}{1.}{*}{\mathrm{zz}}{[}{32}{]}{*}{\mathrm{zz}}{[}{27}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{34}{]}{:=}{\mathrm{zz}}{[}{34}{]}{-}{1.}{*}{\mathrm{zz}}{[}{32}{]}{*}{\mathrm{zz}}{[}{28}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}{1}{/}{\mathrm{zz}}{[}{30}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{31}{]}{:=}{\mathrm{zz}}{[}{31}{]}{*}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{34}{]}{:=}{\mathrm{zz}}{[}{34}{]}{-}{1.}{*}{\mathrm{zz}}{[}{33}{]}{*}{\mathrm{zz}}{[}{31}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{13}{]}{:=}{0.4000000000}{*}{\mathrm{zz}}{[}{52}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{10}{]}{:=}{\mathrm{zz}}{[}{54}{]}{-}{1.}{*}{\mathrm{zz}}{[}{13}{]}{*}{\mathrm{zz}}{[}{43}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{11}{]}{:=}{\mathrm{zz}}{[}{58}{]}{-}{1.}{*}{\mathrm{zz}}{[}{13}{]}{*}{\mathrm{zz}}{[}{46}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{12}{]}{:=}{\mathrm{zz}}{[}{59}{]}{-}{1.}{*}{\mathrm{zz}}{[}{13}{]}{*}{\mathrm{zz}}{[}{49}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{14}{]}{:=}{133.3333333}{*}{\mathrm{zz}}{[}{53}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{11}{]}{:=}{\mathrm{zz}}{[}{11}{]}{-}{1.}{*}{\mathrm{zz}}{[}{14}{]}{*}{\mathrm{zz}}{[}{47}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{12}{]}{:=}{\mathrm{zz}}{[}{12}{]}{-}{1.}{*}{\mathrm{zz}}{[}{14}{]}{*}{\mathrm{zz}}{[}{50}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{15}{]}{:=}{\mathrm{zz}}{[}{10}{]}{*}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{11}{]}{:=}{\mathrm{zz}}{[}{11}{]}{-}{1.}{*}{\mathrm{zz}}{[}{15}{]}{*}{\mathrm{zz}}{[}{29}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{12}{]}{:=}{\mathrm{zz}}{[}{12}{]}{-}{1.}{*}{\mathrm{zz}}{[}{15}{]}{*}{\mathrm{zz}}{[}{32}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{16}{]}{:=}{\mathrm{zz}}{[}{55}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{17}{]}{:=}{\mathrm{zz}}{[}{56}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{18}{]}{:=}{\mathrm{zz}}{[}{57}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{19}{]}{:=}{\mathrm{zz}}{[}{11}{]}{*}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{12}{]}{:=}{\mathrm{zz}}{[}{12}{]}{-}{1.}{*}{\mathrm{zz}}{[}{19}{]}{*}{\mathrm{zz}}{[}{33}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{20}{]}{:=}{\mathrm{zz}}{[}{12}{]}{/}{\mathrm{zz}}{[}{34}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{37}{]}{:=}{\mathrm{zz}}{[}{19}{]}{-}{1.}{*}{\mathrm{zz}}{[}{20}{]}{*}{\mathrm{zz}}{[}{31}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{1}{]}{:=}{\mathrm{zz}}{[}{15}{]}{-}{1.}{*}{\mathrm{zz}}{[}{20}{]}{*}{\mathrm{zz}}{[}{28}{]}{-}{1.}{*}{\mathrm{zz}}{[}{37}{]}{*}{\mathrm{zz}}{[}{27}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{36}{]}{:=}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{2}{]}{:=}{\mathrm{zz}}{[}{14}{]}{-}{1.}{*}{\mathrm{zz}}{[}{20}{]}{*}{\mathrm{zz}}{[}{25}{]}{-}{1.}{*}{\mathrm{zz}}{[}{37}{]}{*}{\mathrm{zz}}{[}{24}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{zz}}{[}{35}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{1}{]}{:=}{\mathrm{zz}}{[}{13}{]}{-}{1.}{*}{\mathrm{zz}}{[}{20}{]}{*}{\mathrm{zz}}{[}{23}{]}{-}{1.}{*}{\mathrm{zz}}{[}{37}{]}{*}{\mathrm{zz}}{[}{22}{]}{-}{1.}{*}{\mathrm{zz}}{[}{36}{]}{*}{\mathrm{zz}}{[}{21}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{2}{]}{:=}{\mathrm{zz}}{[}{2}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{3}{]}{:=}{\mathrm{zz}}{[}{1}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{4}{]}{:=}{\mathrm{zz}}{[}{16}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{5}{]}{:=}{\mathrm{zz}}{[}{17}{]}{;}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{XDOT}}{[}{6}{]}{:=}{\mathrm{zz}}{[}{18}{]}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{end pr}{oc}$ (5.3)