<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta>
      <journal-title-group>
        <journal-title>Kulyasov Nikita[</journal-title>
      </journal-title-group>
    </journal-meta>
    <article-meta>
      <title-group>
        <article-title>Graphical Method for the Knowledge Base Analysis of the Simulation Model*</article-title>
      </title-group>
      <contrib-group>
        <aff id="aff0">
          <label>0</label>
          <institution>Institute of Computational Modelling of the Siberian Branch of the Russian Academy of Sciences</institution>
          ,
          <addr-line>50/44 Akademgorodok, Krasnoyarsk, 660036</addr-line>
          ,
          <country country="RU">Russia</country>
        </aff>
      </contrib-group>
      <volume>000</volume>
      <fpage>0</fpage>
      <lpage>0001</lpage>
      <abstract>
        <p>This article presents a method of analysis of the simulation models which are peculiar due to their graphical structure and declarative knowledge bases containing the description of the state of the modeled object and rules of its functioning. The models consist of intellectual agents, each of those having interfaces of external influence perception and being able to influence other agents through the transmission of control commands via the specified switching links. We have formalized the model and criteria of the analysis, defining operations for the model elements in the form convenient for the software implementation. We suggest the data base structure, providing a combination of graphical methods. Software tools have been completed in the form of a web application of an infographics library. Using this method allows us to analyze the structural links of the model, calculate functional load, and reveal errors in the knowledge bases and elements with lacking or excessive parameters and rules. The results of the method operation are interactive graphical presentations and automatically generated tables of errors and recommendations.</p>
      </abstract>
      <kwd-group>
        <kwd>Simulation Modeling</kwd>
        <kwd>Spacecraft Onboard Equipment</kwd>
        <kwd>Knowledge Base</kwd>
        <kwd>Infographics</kwd>
        <kwd>Agent Modeling</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>-</title>
      <p>
        Modern scientific research has promoted agent-based simulation modeling which is a
presentation of a model in the form of a set of agents and environment of functioning
for the study of the behavior and interaction of complex objects. The science presents
subject-focused tools designed for solving industrial and R&amp;D tasks. For this purpose,
they are supplemented with abstract elements, language constructions and sets of
concepts taken directly from the subject area of research [
        <xref ref-type="bibr" rid="ref1">1</xref>
        ]. Agent-based modeling
does not have a solid set of standard methods of the model development.
* Copyright c 2020 for this paper by its authors. Use permitted under Creative Commons
License Attribution 4.0 International (CC BY 4.0).
      </p>
      <p>
        Widely used are intellectual simulation agents with the system of sensors for
external influence perception. They interpret the received data on the basis of the
embedded knowledge bases and reflect the events affecting the environment with the
help of effectors [
        <xref ref-type="bibr" rid="ref2">2</xref>
        ]. Within the tasks of technical system production support,
knowledge bases provide the accumulation and replication of the experience of highly
qualified specialists in the subject area [
        <xref ref-type="bibr" rid="ref3">3</xref>
        ]. Such an approach is non-deterministic and
allows one to define model functions as a result of the activity its agents. For
graphical construction of the models in different systems the UML language is used,
or a special set of graphical primitives allowing a high level of abstraction and
independence from the method of the model implementation [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. However, these
models have a complicated graphical architecture and a declarative format of the
simulation method presentation which requires the creation of special tools of the
model analysis [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ].
      </p>
      <p>
        Model evaluation methods are mostly about getting information on how well the
model describes real processes happening in the initial object, and how well it will
simulate the development of these processes. Monitoring is performed on the basis of
the statistical estimation of the simulation errors, however if there is no data of the
functional testing of the object, this approach is impossible. In this case, specialists
apply methods that use the qualified experience of experts in the subject area as well
as empirical data and knowledge of the related areas [
        <xref ref-type="bibr" rid="ref6">6</xref>
        ]. They verify the consistency
of knowledge on the basis of graph incidence matrices generalizing relations between
the targets of the knowledge bases, analysis of knowledge reference integrity and
search of the cycle dependences [
        <xref ref-type="bibr" rid="ref7">7</xref>
        ], etc.
      </p>
      <p>
        The purpose of our study is to create a graphical method allowing one to analyze
structural dependences in the knowledge bases of the simulation model in order to
control their completeness, adequacy and consistency. The model elements are
implemented in the form of intellectual agents simulating the logic of the onboard
equipment operation. The models have been designed by specialists of the Institute of
Computational Modeling within the software-and-hardware complex “Software-math
model of the spacecraft command-and-measuring system of onboard equipment” [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ].
This complex is currently a part of a space industry production process.
2
2.1
      </p>
    </sec>
    <sec id="sec-2">
      <title>Graphical Method of Knowledge Base Analysis</title>
      <sec id="sec-2-1">
        <title>Objectives of the Intellectual Simulation Model Design and Analysis</title>
        <p>In order to create a method of intellectual simulation model analysis we have done the
following: model formalization; data base creation; evaluation criteria definition;
graph construction; error analysis; recommendation.</p>
        <p>Model formalization implies the description of its elements and setting the criteria
for analysis, which determine operations for the model elements in the form
convenient for the software implementation. Intellectual simulation model is a set of
intellectual agents describing the condition of the onboard systems and their
functioning at each moment of time. The model elements are the following: B={Bi} –
a set of intellectual agents, I= { } – a set of switching interfaces, С={ } – a set of
links between the agents Bi and Bj. The connection is implemented by switching
interfaces of the agents. Let us set C ijnm=&lt;Iin, Ijm&gt;, where Iin is the interface of the
agent Bi , Ijm is the interface of the agent Bj. The agent contains a status block setting
the values of the onboard equipment characteristics and a behavior block determining
the strategy of its existence. The agent status is described by the following
parameters: X – a set of incoming impacts, K – a set of control commands, Y – a set of
output (monitored) parameters. The agent behavior is defined in a declarative form by
the “condition-action” rules in the knowledge base. The rule R: A→Z, where A = A1&amp;
… &amp;Ar is the logical condition, Z= Z1, …, Zm are the actions changing the model state.</p>
        <p>For the model analysis, we have designed a data base allowing one to store the
graphical model and rules of the agents functioning in universal structures. The data
base contains the following tables: ElementsOfModel – description of the model
elements; Timer – timers; Interface – switching interfaces; Variable – model
variables; Connection – links between the interfaces; Rule – knowledge base rules;
LogicItem – a list of all the functioning devices with the description of actions;
Timer_To_LogicItem, InterfaceToLogicItem, Variable_LogicItem – link tables of
timers, interfaces and variables with the rules. The model contains 1500 thousand of
items describing 13 blocks, 77 switching interfaces, 50 connections, 549 logical
elements in 145 rules of the knowledge base. The complexity of the algorithm of
preprocessing and analysis can be estimated by the expression (  + 5 + 3 + 4 +
3 +  + 2 ), where b is the number of blocks, l is number of conditions/actions in
the rules, i is the number of interfaces, c is the number of connections, r is the number
of rules, and t is the number of timers.</p>
        <p>On the basis of our formalization, and taking into account the given structure of
the data base, we suggest the criteria for analyzing the functional dependences and
control of compliance of the graphical structure of the model with the model
operation methods determined by the knowledge base rules.
2.2</p>
      </sec>
      <sec id="sec-2-2">
        <title>Method of Structural Link Analysis</title>
        <p>This method allows one to estimate the correctness of using the graphical model
elements in the knowledge base and to reveal structural errors. The criteria of
structural link evaluation are as follows:
5. For switching interface the only allowed connection is:
 Iin ! I jm | C ijnm =&lt;Iin, I jm&gt;  C.
6. Commutation connection is acceptable only through one-type interfaces: C ijnm
=&lt;Iin, I jm&gt; Tp(Iin)=Tp(Ijm), where Tp returns the type of the simulated interface.
7. Commutation connection is performed for multidirectional interfaces:  C ijnm
=&lt;Iin, I jm&gt; Rt(Iin)≠Rt(Ijm), where Rt returns the direction by which the data of the
simulated interface is transmitted, and takes the values of «Вх» (incoming) and
«Исх» (outcoming).
8. For the incoming interface, there must be a rule set for the data reception:  IinI
Rt(Iin)=«Вх» ⇒ Sel(Ai, I=Ii) ≠, where Sel is the function of selection from the set
of the elements with the given properties.
9. For the outcoming interface there must be a rule set for the data transmission:
 IinI Rt(Iin)=«Исх» ⇒ Sel(Zi, I=Ii) ≠.
10. For the interaction of the model elements specified by the rules, commutation
connections must be determined:  L(I in , I i+lm) ⇒СL={Cpq=&lt;Ip, Iq&gt;, p=i, …,
i+l-1, q=p+1}, where L(I in , I i+lm) is the path between the agents Bi , Bi+1, …, Bi+l
via the links set in СL.</p>
        <p>The criteria of the analysis are presented in a circle diagram of the structural links.
Fig.1 demonstrates a part of this diagram with the dependent interfaces.</p>
        <p>The diagram sections are the model elements which are the simulators of the onboard
devices, and the rays between them show the relations specified in the knowledge
base. The Figure shows the commutation connections between the model elements:
CCU 1, CCU 2 (simulations of the command-and-measuring system) and RECEIV
(simulator of the receiving device). The ray width shows the direction of the data
transmission.</p>
        <p>Graphical visualization helps to reveal lacking or excessive structures, as well as
model elements for which no rules have been specified in the knowledge base. This
allows finding inconsistency between the graphical presentation and the knowledge
base.
2.3</p>
      </sec>
      <sec id="sec-2-3">
        <title>Method of Functional Load Analysis</title>
        <p>Examination of the functional load is performed by the following criteria:
≤
|Sel(R, I=In)| ≤
1. For all the elements of the model, the methods of functioning must be set:  Bi </p>
        <p>Ri≠, where Ri denotes the rules of operation of Bi.
2. Equal functional load for one-type model elements must be provided:  BiB |Ri|
∑| | | |</p>
        <p>| |
load excess coefficient.
3. Equal functional load for one-type model interfaces must be provided:  InI
∑| | | ( , )|
+  , where |B| is the capacity of the set,  is the allowable functional
| |</p>
        <p>+  .</p>
        <p>The analysis is performed graphically. The nodes in the graph are the model
elements, and the curves are the ways of their interaction with other sub-models
(Fig. 2). The Figure shows the following model elements: CCU 1, CCU 2
(commandand-measuring system); OCS CU (onboard control complex); ODCC (onboard digital
computing complex); ODGS, ODGS 2 (onboard remote signaling equipment),
RECEIV 1, RECEIV 2 (receiving devices), TRANS 1, TRANS 2 (transmitting
devices), RECEIV.AERIAL, TRANS.AERIAL (antennas); GCC (ground segment).</p>
        <p>The size of the circles and curves demonstrates the amount of load calculated as the
capacity of the set of rules, in which the model elements participate. High load for the
model elements may become a reason for its revision or for the necessity of additional
equipment and commutation.
2.4</p>
      </sec>
      <sec id="sec-2-4">
        <title>Method of Complete Data Processing Analysis</title>
        <p>The control of consistency is performed for the operation parameters and commands
specified in the model description. There are the following criteria of the analysis:
4. There are rules for filling the model parameters: Pr(Ri, Xi)≠, where Ri is a set of
the rules Bi, Xi are the input parameters of Bi, Pr is the projection function
choosing all the elements from the set Xi, described in the rules.
5. There are rules of the model parameters interpretation: Pr(Ri, Yi)≠, where Yi are
the output parameters of Bi.
6. For each command k there is a model element performing its transmission:  kK
 Bi и  R1: A1→Z1, R1Sel(Ri, K= k) | kPr(Z1, K).
7. For each command k there is a model element performing its reception:  kK  Bj
и  R2: A2→Z2, R2Sel(Rj, K= k) | kPr(A2, K).</p>
        <p>The model implementation criteria are verified with the help of the coverage graph
(Fig. 3), which demonstrates the model elements containing the nodes with no rules,
or where no actions regarding the command reception or transmission have been
specified for the interfaces.</p>
        <p>The Figure shows the model elements listed in the load schedule and contains the
details on switching links. This interactive graphical tool allows one to analyze in
detail the implementation criteria with regard to the data for the model as a whole, as
well as for each of the model, or switching interfaces.</p>
        <p>The following notation is used: green – the logic of reception/transmission is
described for all switching interfaces of an element; red – all interfaces of an element
have no description of the communication logic; orange – one or more interfaces in
the element does not have any logic.
2.5</p>
      </sec>
      <sec id="sec-2-5">
        <title>Method of Model Coherence Analysis</title>
        <p>For the analysis of the functional dependences and in order to determine the model
working modes specified in the knowledge bases of the intellectual agents we suggest
the following criteria:
8. Each rule of the knowledge base must be included in one of the model functioning
modes:  RijRi  FModR | Sel(FMod, R=Rij) ≠, where FMod is a subset of the
logical chains of the logical Inference.
9. Between the model elements with the dependent rule chains, commutation links
must be set: Dep(Rchi, Rchj) ⇒  R1: A1→Z1, R2: A2→Z2, R1Rchi, R2Rchj, 
Cij=&lt;Ii1, Ij2&gt;, Ii1Pr(Z1, I), Ij2Pr(A2, I), where Dep defines the dependence
relation of the rule chains Rchi and Rchj.
10. For all the interfaces, the rules of data reception or transmission must be set: </p>
        <p>IinI Sel(Ri, I=Iin)≠.
11. If there is a commutation link between the model elements, there exists a way of
the data transmission via the interfaces included in this link:  Cijnm=&lt;Iin, I jm&gt; 
L(I in , I jm).</p>
        <p>Graphical interpretation of the specified criteria of the analysis is performed with
the help of a coherence graph, describing the interaction of the model elements
(Fig. 4). The graph shows the model elements for each switching interface, for
example, «CCU1 to ODGS» – interface CCU1 (command-and-measuring system),
which sends data to ODGS (onboard remote signaling equipment), «CCU1 from
ODGS» – interface CCU1, which receives data from ODGS. The graph is divided
into unrelated sections, for which the rules describe the events happening under
certain conditions independently from other processes. If the knowledge base
specifies multiple interactions between the model elements, the graph demonstrates
all possible ways of data transmission, regardless of the conditions of the rule
initiation. The graph gives specialist in the subject area a clear view of the simulation
model operation modes set in the knowledge base.</p>
        <p>
          In order to present the details the dependences realized during the logical output,
rule chains are built, where each rule is a separate node initiated under the specified
conditions. Such chains form precedents of simulation modeling and they are used for
the model analysis [
          <xref ref-type="bibr" rid="ref11">11</xref>
          ] through the comparison with the results of the tests
implemented on the simulated devices.
        </p>
      </sec>
    </sec>
    <sec id="sec-3">
      <title>The Application Results of Graphical Method of Analysis</title>
      <p>
        The criteria of the analysis have been developed in the form of a web-application on
the basis of the following software means and components: reception, pre-processing,
analysis and interpretation of data for visualization, which are implemented on the
server part of the application with the help of PHP7 scripts. At the core of the client’s
part, infographics libraries d3.js and sigma.js are used [
        <xref ref-type="bibr" rid="ref9">9</xref>
        ], providing interactive
graphical tools for visualization of the functional dependences [
        <xref ref-type="bibr" rid="ref10">10</xref>
        ].
      </p>
      <p>As a result of the analysis, tables of errors and recommendations are automatically
generated. The table of recommendations contains data on the excess of the average
level of the functional load for model elements (Fig 5).</p>
      <p>The table of errors contains a list of the model elements and interfaces with the
description of the errors found in the model structure or in the knowledge base. The
table of recommendations provides a list of the elements, switching interfaces and
connections for which the load is above the average calculated for all the model
elements. It also gives recommendations on additional devices or communication
lines. The model revision including the correction of errors and consideration of
recommendations increases the quality of simulation tests, forming the base for the
efficient task solution in the subject area.
4</p>
    </sec>
    <sec id="sec-4">
      <title>Conclusion</title>
      <p>The creation of an intellectual model is an important scientific task for the efficient
support of the design of complex technical systems. However, in order to effectively
use models and conduct simulation tests it is necessary to provide high quality and
adequacy of the constructed models. Our graphical method of the simulation model
analysis is designed to solve this problem. The method allows one to reveal the
dependences of the model elements, errors of the knowledge base, lacking or
excessive data and structures, for which no rules have been specified.</p>
      <p>High representativeness is an advantage of our method. It provides a specialist in
the subject area not only with automatic functions of the model control, but also with
the tools of infographics allowing the examination of the simulated devices and
checking of the compliance of the models with the technical tasks and design
documentation. Implementation of the web technology provides simultaneous work of
independent expert groups of different specialization. Owing to clarity, the acquired
graphs can be used for the transfer of unique knowledge and for education of
specialists via detailed immersion into the subject area.</p>
      <p>The future development of the graphical method implies the study of structures of
the model storage and its adaptation for other implementations. The distribution of
analytical tools to the related subject areas demands the creation of additional tools of
transformation of the existing formats of model storage into universal structures of the
developed data base being an intermediate representation between model and
software.</p>
      <p>Acknowledgments. The reported study was funded by RFBR and Government of
Krasnoyarsk Territory according to the research project № 18-47-242007 «The
technology of intellectual support of the spacecraft onboard systems design on the
basis of heterogeneous simulation models».</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          1.
          <string-name>
            <surname>Chung</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          : Simulation Modeling Handbook. London, CRC Press (
          <year>2004</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          2.
          <string-name>
            <surname>Russel</surname>
            ,
            <given-names>S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Norvig</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          : Artificial Intelligence:
          <string-name>
            <given-names>A Modern</given-names>
            <surname>Approach</surname>
          </string-name>
          . 3rd edn. Prentice Hall, New Jersey (
          <year>2010</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          3.
          <string-name>
            <surname>Tan</surname>
            ,
            <given-names>C.F.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Wahidin</surname>
            ,
            <given-names>L.S.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Khalil</surname>
            ,
            <given-names>S.N.</given-names>
          </string-name>
          and
          <string-name>
            <surname>Tamaldin</surname>
            ,
            <given-names>N.:</given-names>
          </string-name>
          <article-title>The application of expert system: a review of research and applications</article-title>
          .
          <source>ARPN Journal of Engineering and Applied Sciences</source>
          <volume>11</volume>
          (
          <issue>4</issue>
          ),
          <fpage>2448</fpage>
          -
          <lpage>2453</lpage>
          (
          <year>2016</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>4. AnyLogic 7.2 Released. AnyLogic Simulation Software, www.anylogic.com</mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          5.
          <string-name>
            <surname>Zhou</surname>
            ,
            <given-names>Y.</given-names>
          </string-name>
          :
          <article-title>An intelligent model validation method based on ECOC SVM</article-title>
          .
          <source>In Proceedings of the 10th International Conference on Computer Modeling and Simulation</source>
          , pp.
          <fpage>67</fpage>
          -
          <lpage>71</lpage>
          Association for Computing Machinery,
          <string-name>
            <surname>NY</surname>
          </string-name>
          , USA (
          <year>2018</year>
          ). doi:
          <volume>10</volume>
          .1145/3177457.3177487
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          6.
          <string-name>
            <surname>Min</surname>
            ,
            <given-names>F.</given-names>
          </string-name>
          :
          <article-title>Knowledge-based method for the validation of complex simulation models</article-title>
          .
          <source>Simulation modelling practice and theory</source>
          <volume>18</volume>
          (
          <issue>5</issue>
          ),
          <fpage>500</fpage>
          -
          <lpage>515</lpage>
          (
          <year>2010</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          7.
          <string-name>
            <surname>Zanon</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          :
          <article-title>The SimTG simulation modeling framework a domain specific language for space simulation</article-title>
          .
          <source>In. Proceedings of the Symposium on Theory of Modeling</source>
          &amp; Simulation
          <string-name>
            <surname>: DEVS Integrative M&amp;S Symposium</surname>
          </string-name>
          , pp.
          <fpage>16</fpage>
          -
          <lpage>23</lpage>
          (
          <year>2011</year>
          )
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          8.
          <string-name>
            <surname>Nozhenkova</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Isaeva</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Markov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Koldyrev</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Vogorovskiy</surname>
            ,
            <given-names>R.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Evsyukov</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Simulation infrastructure design on the basis of the space industry's international standards</article-title>
          .
          <source>Advances in Intelligent Systems Research</source>
          <volume>134</volume>
          ,
          <fpage>138</fpage>
          -
          <lpage>141</lpage>
          (
          <year>2017</year>
          ). doi:
          <volume>10</volume>
          .2991/caai-
          <fpage>17</fpage>
          .
          <year>2017</year>
          .28
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          9.
          <string-name>
            <surname>Bostock</surname>
            ,
            <given-names>M.</given-names>
          </string-name>
          :
          <string-name>
            <surname>Data-Driven</surname>
            <given-names>Documents</given-names>
          </string-name>
          , https://d3js.org/,
          <source>last accessed</source>
          <year>2019</year>
          /04/28
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          10.
          <string-name>
            <surname>Kulyasov</surname>
          </string-name>
          , N.:
          <article-title>Method of creation and verification of the spacecraft onboard equipment operation model</article-title>
          .
          <source>IOP Conference Series: Materials Science and Engineering</source>
          <volume>537</volume>
          ,
          <fpage>1</fpage>
          -
          <lpage>6</lpage>
          (
          <year>2019</year>
          ). doi:
          <volume>10</volume>
          .1088/
          <fpage>1757</fpage>
          -899X/537/2/022042
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          11.
          <string-name>
            <surname>Nozhenkova</surname>
            ,
            <given-names>L.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Isaeva</surname>
            ,
            <given-names>O.</given-names>
          </string-name>
          ,
          <string-name>
            <surname>Koldyrev</surname>
            ,
            <given-names>A.</given-names>
          </string-name>
          :
          <article-title>Creation of the base of a simulation model's precedents for analysis of the spacecraft onboard equipment testing results</article-title>
          .
          <source>Advances in Intelligent Systems Research</source>
          <volume>151</volume>
          ,
          <fpage>78</fpage>
          -
          <lpage>81</lpage>
          (
          <year>2018</year>
          ). doi:
          <volume>10</volume>
          .2991/cmsa-
          <fpage>18</fpage>
          .
          <year>2018</year>
          .18
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>