Note: Claims are shown in the official language in which they were submitted.
-13-
What is claimed is:
1. A computational model of the heart comprising:
a space defining lattice having a set of nodes;
each of said nodes positioned in said lattice such that each node is
adjacent to neighboring nodes;
whereby said lattice and node together define a three dimensional
representation of at least a portion of the heart;
each node of said set of nodes representing a unit of the myocardium;
each node having associated with it a set of state defining node
equations for computing an action potential at the location of said node
and for computing a node potential;
each node having associated with it a set of coupling relationships
related to the anatomic structure of the heart for computing the
contribution to said node potential contributed by neighboring nodes;
whereby the total set of voltages at each of said nodes represents the
global depolarization state of the heart.
2. A computational model of the heart comprising:
a cubic lattice having a set N nodes;
each of said nodes positioned at a vertex of said lattice;
each node representing a biophysical computational unit of the heart;
each node having associated with it a set of state defining equations
expressing state variables sufficient to compute node response at time
t+deltat given node response at time t, said equations including at least
one equation selected from biophysical processes selected from the set:
i) equations for time-varying, voltage-dependent transmembrane
conductances permeable to sodium (Na), potassium (K), calcium (Ca),
and/or chloride (CI) ions, the temporal evolution of which are modeled
using ordinary differential equations;
-14-
ii) equations for time-varying transmembrane ionic pump and
exchanger currents, properties of which are modeled using algebraic
equations;
iii) equations for time-varying Ca uptake, sequestration, and release
currents modeling the regulation of intracellular Ca levels by cellular
organelles, the properties of which are modeled using coupled systems of
ordinary differential equations;
iv) equations for total transmembrane flux of each ionic species to
which the cell membrane is permeable;
each node connection to neighboring nodes represented by a coupling
conductance, the value of which is defined by the anatomic relationship of
the two nodes;
each node having a coupling current associated with it determined by
the total current entering the node through all said coupling conductances
defined at that node;
whereby said coupling currents and node state equations may be solved
for the temporal evolution of all state variables at all nodes.
3. A process for computing a model of a heart comprising the steps of:
a) defining a set of lattice nodes;
each lattice node representing a biophysical subunit of the heart;
each node having a set of biophysical equations for computing a total
node current associated with Na, Ca and K ion transport in said biophysical
subunit;
each node associated with at least five adjacent nodes;
each node coupled to neighboring nodes;
said coupling modeled as a resistance;
b) solving said set of biophysical equations to determine said node
currents;
c) summing said node current to determine the voltage present at each
node resulting from the depolarization of each neighboring node;
-15-
d) displaying a representation of said computed voltage.
4. A method for modeling a heart with a computer comprising:
a) defining a set of nodes;
each node representing a biophysical subunit of the heart;
each node having a membrane, defining an intracellular and an
extracellular space;
each node having a set of biophysical equations associated with it for
computing a coupling current associated with ion transport in said
intracellular space;
each node having a set of biophysical equations associated with it for
computing a coupling current associated with ion transport in said
extracellular space;
each node having a set of biophysical equations associated with it for
computing a transmembrane current between the intracellular and
extracellular spaces;
b) defining a lattice including each of said nodes;
each node associated with a physical location in heart tissue;
each node connected with at least five adjacent nodes;
each node exhibiting an anisotropic coupling with neighboring nodes
reflecting the anatomic relationship between such nodes;
said coupling modeled as a coupling resistance
c) solving said set of biophysical equations to determine said node
currents associated with extracellular and intracellular currents;
d) summing said node current to compute the voltage present at each
node resulting from the depolarization of each node and the propagation of a
depolarization wave front through each node;
e) displaying a representation of said computed voltage.
-16-
5. A model of a biological organ system comprising:
a plurality of nodes;
each node having associated therewith a set of state defining equations
which take their variables biophysical data, which are substantially local to
the node itself;
each node having at least one state variable which takes as its argument
parameters communicated solely from near neighbors;
a network containing all such nodes;
said network including coupling relationships defined between each
node and its near neighbor node wherein said coupling relationships reflect
the anatomical structure of the biological organ, and are expressed as set of
state defining equations;
whereby said state defining variables can be solved interatively and
state variable values communicated to the network model, whereby said
network equations can be computed for each node.
6. A method of defining a model comprising the steps of:
applying nuclear magnetic resonance imaging to a biological organ to
generate a data file defining tensors associated with the distribution of
water
molecules within the physiologic organ;
applying said tensor dataset to a network to determine the magnitude of
coupling conductances between nodes of the network such that the tensor
information is reflected by the internodal conduction values of the network,
thus embedding the anatomic characteristics of the organ in the model.
7. A method for using a multiple processor computer to solve for the
temporal evolution of state variables defining biophysical and biochemical
properties of nodes within a network model of a biological system
comprising:
-17-
partitioning the lattice of nodes making up the network into N number
of lattice subsets, such that each node within the network is included in at
least one of said subsets;
associating state equations and state variables describing the
biophysical and biochemical properties of each node in said subsets with
particular processing units in the said parallel computer such that each
subset
of nodes is associated with a distinct processing unit, and that all
processing
units are associated with at least one subset of nodes;
computing the temporal evolution of all state variables within each of
the node subsets associated with distinct processing units concurrently across
each of the processing units;
concurrently storing said state variables in computer memory;
displaying a sequence of said state variables.
8. A method of modeling a cardiac disease states in a model of the heart as
set forth in claim 3 comprising the additional steps of:
a) adjusting the lattice defining the anatomic structure of the heart
such that anatomical structural changes known to occur in the heart during
said disease state;
b) using said adjusted tensor dataset to determine the magnitude of
coupling conductances between nodes of the network such that the altered
structure of the heart in the disease state, as represented by the adjusted
tensor
dataset, is embedded in the network model;
c) adjusting parameters of the set of state defining equations,
specifying properties of biochemical cellular processes defined at each node
of the lattice, to reproduce changes in these parameters known to occur in
said
disease state;
d) adjusting initial values of state variables, the temporal evolution of
which is determined by the set of equations specifying properties of
biochemical cellular processes defined at each node of the lattice, in such a
-18-
way as to reproduce changes in these state variables known to occur in said
disease state;
e) solving for the temporal evolution of state variables, defined at
each node of the lattice;
f) storing said state variables in a computer memory for graphical
display and analysis.
9. A computational model of the heart comprising:
a set of n nodes arranged in a cubic lattice with each node having at
least five neighbors;
each node having a set of state variable equations associated with it;
said equations describing biophysical and biochemical reactions at
each node which may be solved to find the currents present at each node;
the lattice having a set of coupling equations associated therewith
defining the coupling between each node which may be solved to determine
the total voltage present at each node in the lattice.
10. A finite-difference computations model of a heart comprising:
a) a set of N nodes arranged in a cubic lattice with each node having at
least five neighbors;
b) a subroutine procedure which uses experimentally measured data on
fiber orientation within the heart (obtained in a number of ways, including
direct anatomical measurements or use of magnetic resonance imaging
methods) for computing the connection currents between lattice nodes;
c) a multi-dimensional data array S(NI,N2,N3,S1,S2,S3, ... SNeq) in
which indices N1, N2, and N3 specify position of a node within the cubic
lattice, and variables S1.... SNeq are values of state variables defined at
node
(NI,NZ,N3) at time t;
d) a multi-dimensional data array F(N1,N2,N3,S1- S2-S3- .. SNeq-) in
which the variable Sx~ denotes the time rate of change, dSx(t)/dt, of state
variable x at time t;
-19-
e) a subroutine procedure which specifies an algorithm for computing
elements of the multi-dimensional data array F(), with this subroutine
procedure itself being composed of a number of subroutine procedures
selected from a library of procedures corresponding to biophysical processes
modeled at the sub-cellular level;
f) a subroutine procedure specifying a numerical integration algorithm
for the computation of values of the array So at time t+At given values of So
and Fo at time t.
11. A computational model of a heart comprising:
a) a set of n nodes arranged in a cubic lattice with each node having at
least five neighbors;
b) each node having a set of state variable equations associated with
it;
c) said equations describing biophysical and biochemical reactions at
each node which may be solved to find the currents present at each node
based upon local parameters, and the voltages present at each node based
upon local and near neighbor parameters;
d) the lattice having a set of coupling equations associated therewith
defining the coupling between each node, which may be solved to determine
the total node voltage present at each node in the lattice.
12. A method of computing the model of claim 8 comprising the steps of:
i) partitioning the lattice into m nodes and assigning state variable
computations to one of (p) processing units;
ii) concurrently computing state variables with the p processing
units;
iii) computing the temporal evolution of all state variables and
storing or displaying the results.
-20-
13. A method of computing the temporal evolution of state variables in the
model of claim 3 on multiprocessor computers consisting of a set of P
processors each with local memory of size M, a communications network
supporting data exchange between processors, and a global shared memory,
with the method comprising the steps of:
a) partitioning the lattice into Q sets of nodes, referred to as node sets,
in such a way that all state variable values at time t, parameters, and
required
temporary storage locations (work space) for any node set fits into the local
memory available to any processing unit;
b) assignment of each node set to a specific processing unit and it's
local memory;
c) inter-processor communication of any data required to compute new
values of state variables for nodes positioned at the borders between
different
node sets;
d) concurrent execution on all P processors of the computations to
compute new values of local state variables at time t+Dt in each of the node
sets that have been assigned to processing units;
e) iteration of steps b) and d) until computation of state variable
values at time t+Dt, where Dt is small, for all Q node sets are completed;
f) storage of the computed set of state variables on external storage
devices for subsequent graphical display;
g) iteration of steps b) - f) until t+Dt equals some end time T.
14. The system of claim 4 with the inclusion of chemical or compound data:
a) a means for modifying one or more of said state variables as a
result of imposing one of said chemical or compound data;
chemical or compound data includes mechanism of action data; kinetic data;
biochemical, mechanical;
b) a means for solving and updating the complete set of integral
network algorithms and chemical/composition data in time sets;
-21-
c) a means to determine the currents present at each node based on
changes in the local and near neighbor state variable parameters as a result
of
chemical/compound data;
d) a means to determine the action potential each node in the lattice;
e) a means for modifying state variables as a result of
chemical/compound data;
f) a means to use parallel computers to solve for the temporal
evolution of state variables and chemical/compound data defined at each node
of the lattice;
g) a means to graphically display the changes at each node in the
lattice as a result of said modification of the system of Claim 4 with
chemical/compound data;
h) a means to store said state variables in computer memory and or
external data storage devices for graphical display and analysis.
15. The system of claim 4 with electrical input data comprising:
a) a means for modifying one or more of said state variables as a
result of imposing one of said electrical input data;
electrical input data includes electric waveform stimulation protocols, etc.
b) a means for solving and updating the complete set of integral
network algorithms and electrical input data in time sets;
c) a means to determine the currents present at each node based on
changes in the local and near neighbor state variable parameters as a result
of
electrical input data;
d) a means to localize the origin the origin of the electrical stimulus;
e) a means to determine the ecg at each node in the lattice;
f) a means for modifying state variables as a result of electrical input
data;
g) a means to use parallel computers to solve for the temporal
evolution of state variables and electrical input data defined at each node of
the lattice;
-22-
h) a means to graphically display the changes at each node in the
lattice as a result of said modification of the system of Claim 4 with
electrical
input data;
i) a means to store said state variables in computer memory and or
external data storage devices for graphical display and analysis.
16. The model of claim 4 further comprising:
a) a means for modifying one or more of said state variables as a
result of imposing chemical compound data in said state equations, wherein
said chemical compound data includes mechanism of action data; kinetic data;
biochemical, or mechanical data;
b) a means for solving and updating the complete set of integral
network algorithms and chemical/composition data in time sets;
c) a means to determine the currents present at each node based on
changes in the local and near neighbor state variable parameters as a result
of
chemical/compound data;
d) a means to determine the ecg at each node in the lattice;
e) a means for modifying state variables as a result of
chemical/compound data;
f) a means to use parallel computers to solve for the temporal
evolution of state variables and chemical/compound data defined at each node
in the lattice;
g) a means to graphically display the changes at each node in the
lattice as a result of said modification of the system with chemical compound
data;
h) a means to store said state variables in computer memory and or
external data storage devices for graphical display and analysis.
1. A computational model of the heart comprising:
a space defining lattice having a set of nodes;
each of said nodes positioned in said lattice such that each node is adjacent
to
neighboring nodes;
whereby said lattice and node together define a three dimensional
representation of at least a portion of the heart;
each node of said set of nodes representing a biophysical subunit of the
myocardium;
each node having associated with it a set of state defining node equations for
computing an action potential at the location of said node based upon a
biophysical model that explicitly relies on the computation of individual
ionic
current that collectively give rise to said action potential, and for
computing a
node potential at the node from the action potentials within the biophysical
subunit derived directly from the computed action potentials;
each node having associated with it a set of coupling relationships related to
the anatomic structure of the heart for computing the contribution to said
node
potential contributed by the voltages and currents exchanged with the
neighboring nodes;
whereby the total set of voltages at each of said nodes represents the global
depolarization state of the heart.
2. A computational model of the heart comprising:
a cubic lattice having a set N nodes;
each of said nodes positioned at a vertex of said lattice;
each node representing a biophysical computational unit of the heart;
-24-
each node having associated with it a set of state defining equations
expressing state variables sufficient to compute node response at time
t+deltat
given node response at time t, said equations including at least one equation
selected from biophysical processes selected from the set:
(i) equations for time-varying, voltage-dependent transmembrane
conductance permeable to sodium (Na), potassium (K), calcium (Ca),
and/or chloride (CI) ions, the temporal evolution of which are molded
using ordinary differential equations;
(ii) equations for time-varying transmembrane ionic pump and exchanger
currents, properties of which are modeled using algebraic equations;
(iii) equations for time-varying Ca uptake, sequestration, and release
currents modeling the regulation of intracellular Ca levels by cellular
organelles, the properties of which are modeled using coupled systems
of ordinary differential equations;
(iv) equations for total transmembrane flux of each ionic species to which
the cell membrane is permeable;
each node connection to neighboring nodes represented by a coupling
conductance, the value of which is defined by the anatomic relationship of the
two nodes;
each node having a coupling current associated with it determined by
the total current entering the node through all said coupling conductance
defined at that node;
whereby said coupling currents and node state equations may be
solved for the temporal evolution of all state variables at all nodes.
3. A process for computing a model of a heart comprising the steps of:
(a) defining a set of lattice nodes;
-25-
each lattice node representing a biophysical subunit of the heart;
each node having a set of biophysical equations for computing a total
node current associated with Na, Ca and K ion transport across cellular
membranes in said biophysical subunit;
each node associated with at least five adjacent nodes;
each node coupled to neighboring nodes;
said node coupling modeled as a resistance between adjacent nodes;
(b) solving said set of biophysical equations to determine said node currents;
(c) summing said node current to determine the voltage present at each node
resulting from the depolarization of each node in conjunction with each
neighboring node;
(d) displaying a representation of said computed voltage to show the time
course
of changes in voltages at nodes.
4. A method for modeling a heart with a computer comprising:
(a) defining a set of nodes;
each node representing a biophysical subunit of the heart;
each node having a membrane, defining an intracellular and an extracellular
space;
each node having a set of biophysical equations associated with it for
computing a coupling current associated with ion transport in said
intracellular
space;
each node having a set of biophysical equations associated with it for
computing a coupling current associated with ion transport in said
extracellular
space;
(b) defining a lattice including each of said nodes;
-26-
each node associated with a physical location in heart tissue;
each node connected with at least five adjacent nodes;
each node exhibiting an anisotropic coupling with neighboring nodes
reflecting the anatomic relationship between such nodes;
said coupling modeled as a coupling resistance
(c) solving said set of biophysical equations to determine said node currents
associated with extracellular and intracellular currents;
(d) summing said node current to compute the voltage present at each node
resulting from the depolarization of each node and the propagation of a
depolarization wave front through each node;
(e) displaying a representation of said computed voltage.
5. A model of a biological organ system comprising:
A plurality of nodes;
each node having associated therewith a set of state defining equations which
take their variables biophysical data, which are substantially local to the
node itself,
and which include explicit computations for ionic flux across cellular
membranes for
giving rise to action potentials at the node;
each node having at least one state variable which takes as its argument
parameters communicated solely from near neighbors which is either current,
voltage;
a network containing all such nodes;
said network including coupling relationships defined between each node where
each coupling relationship is modeled as a conductance which permits voltage
and
current summations at each node, and its near neighbor node wherein said
coupling
-27-
relationships reflect the anatomical structure of the biological organ, and
are
expressed as set of state defining equations;
whereby said state defining variables can be solved interactively and state
variable values communicated to the network model, whereby said network
equations can be computed for each node.
6. A method of defining a model comprising the steps of:
applying nuclear magnetic resonance imaging to a biological organ to
generate a data file defining tensors associated with the distribution of
water
molecules within the physiologic organ;
applying said tensor dataset to a network to determine the magnitude of
coupling conductance between nodes of the network such that the tensor
information is reflected by the internodal conduction values between the nodes
of the network, thus embedding the anatomic characteristics of the organ in
the
model by representing the coupling relationships as conductanes between
neighboring nodes.
7. A method for using a multiple processor computer to solve for the temporal
evolution of state variables defining biophysical and biochemical properties
of
nodes within a network model of a biological system comprising:
partitioning the lattice of nodes making up the network into N number of
lattice subsets, such that each node within the network is included in at
least on of
said subsets;
associating state equations and state variables describing the biophysical and
biochemical properties of each node in said subsets with particular processing
units
in the said parallel computer such that each subset of nodes is associated
with a
-28-
distinct processing unit, and that all processing units are associated with at
least one
subset of nodes;
computing the temporal evolution of all state variables within each of the
node subsets associated with distinct processing units concurrently across
each of
the processing units;
concurrently storing said variables in computer memory;
displaying a sequence of said state variables.
8. A method of modeling a cardiac disease states in a model of the heart as
set forth
in claim 3 comprising the additional steps of:
(a) adjusting the lattice defining the anatomic structure of the heart such
that
anatomical structural changes known to occur in the heart during said disease
state;
(b) using said adjusted tensor dataset to determine the magnitude of coupling
conductance between nodes of the network such that the altered structure of
the heart in the disease state, as represented by the adjusted tensor dataset,
is
embedded in the network model;
(c) adjusting parameters of the set of state defining equations, specifying
properties of biochemical cellular processes defined at each node of the
lattice, to reproduce changes in these parameters known to occur in said
disease state;
(d) adjusting initial values of state variables, the temporal evolution of
which is
determined by the set of equations specifying properties of biochemical
cellular ionic current exchange processes defined at each node of the lattice,
in
such a way as to reproduce changes in these state variables known to occur in
said disease state;
-29-
(e) solving for the temporal evolution of state variables, defined at each
node of
the lattice;
(f) storing said state variables in a computer memory for graphical display
and
analysis.
9. A computational model of the heart comprising:
a set of N nodes arranged in a cubic lattice with each node having at least
five
neighbors;
each node having a set of state variables equations associated with it;
said equations describing biophysical and biochemical reactions at each node
which may be solved to find the currents present at each node based upon
computed
ionic currents across cellular membranes;
the lattice having a set of coupling equations associated therewith defining
the coupling conductance between each node which may be solved to determine
the
total voltage present at each node in the lattice derived from voltages
generated
within the node and communicated to the node from near neighbors.
10. A finite-difference computations model of the heart comprising:
(a) a set of N nodes arranged in a cubic lattice with each node having at
least five
neighbors;
(b) a subroutine procedure which uses experimentally measured data on fiber
orientation within the heart (obtained in a number of ways, including direct
anatomical measurements or use of magnetic resonance imaging methods) for
computing the connection currents between lattice nodes;
(c) a multi-dimensional data array S(N1,N2,N3,S1,S2,S3, ... Sneq) in which
indices N1, N2, and N3 specify position of a node within the cubic lattice,
and
-30-
variables S1 .... Sneq are values of state variables defined at node
(N1,N2,N3)
at time t;
(d) a multi-dimensional data array F(N1,N2,N3,S1-S2-S3-.Sneq-) in which the
variable Sx' denotes the time rate of change, dSx(t)/dt, of state variable x
at
time t;
(e) a subroutine procedure which specifies an algorithm for computing elements
of the multi-dimensional data array F~, with the subroutine procedure itself
being composed of a number of subroutine procedures selected from a library
of procedures corresponding to biophysical processes modeled at the sub-
cellular level;
(f) a subroutine procedure specifying a numerical integration algorithm for
the
computation of values of the array So at time t+At given values of So and Fo
at time t.
11. A computational model of a heart comprising:
(a) a set of n nodes arranged in a cubic lattice with each node having at
least five
neighbors;
(b) each node having a set of state variable equations associated with it;
(c) said equations describing biophysical and biochemical reactions at each
node
which may be solved to find the currents present at each node based upon
local parameters, and the voltages present at each node based upon local and
near neighbor parameters;
(d) the lattice having a set of coupling equations associated therewith
defining
the coupling between each node, which may be Solved to determine the total
node voltage present at each node in the lattice.
-31-
12. A method of computing the model of claim 8 comprising the steps of:
(i) partitioning the lattice into m nodes and assigning state variable
computations to one of (p) processing units;
(ii) concurrently computing state variables with the p processing units;
(iii) computing the temporal evolution of all state variables and storing or
displaying the results.
13. A method of computing the temporal evolution of state variables in the
model of
claim 3 on multiprocessor computers consisting of a set of P processors each
with
local memory of size M, a communications network supporting data exchange
between processors, and a global shared memory, with the method comprising
the steps of:
(a) partitioning the lattice into Q sets of nodes, referred to as node sets,
in such a
way that all state variable values at time t, parameters, and required
temporary storage locations (work space) for any node set fits into the local
memory available to any processing unit;
(b) assignment of each node set to a specific processing unit and it's local
memory;
(c) inter-processor communication of any data required to compute new values
of state variables for nodes positioned at the borders between different node
sets;
(d) concurrent execution on all P processors of the computations to compute
new
values of local state variables at time t+Dt in each of the node sets that
have
been assigned to processing units;
(e) iteration of steps b) and d) until computation of state variable values at
time
t+Dt, where Dt is small, for all Q node sets are completed;
-32-
(f) storage of the computed set of state variables on external storage devices
for
subsequent graphical display;
(g) iteration of steps b) - f) until t-Dt equals some end time T.
14. The system of claim 4 with the inclusion of chemical or compound data
which
has an impact on the ionic currents computed for a node:
(a) a means for modifying one or more or said state variables as a result of
imposing one of said chemical or compound data;
chemical or compound data includes mechanism of action data; kinetic data;
biochemical, mechanical;
(b) a means for solving and updating the complete set of integral network
algorithms and chemical/composition data in time sets;
(c) a means to determine the currents present at each node based on changes in
the local and near neighbor state variable parameters as a result of
chemical/compound data;
(d) a means to determine the action potential each node in the lattice;
(e) a means for modifying state variables as a results of chemical/compound
data;
(f) a means to use parallel computers to solve for the temporal evolution of
state
variables and chemical/compound data defined at each node of the lattice;
(g) a means to graphically display the changes at each node in the lattice as
a
result of said modification of the system of Claim 4 with chemical/compound
data;
(h) a means to store said state variables in computer memory and or external
data
storage devices for graphical display and analysis.
-33-
15. The system for claim 4 with electrical input data comprising:
(a) a means for modifying one or more of said state variables as a results of
imposing one of said electrical input data;
electrical input data includes electric waveform stimulation protocols, etc.
(b) a means for solving and updating the complete set of integral network
algorithms and electrical input data in time sets;
(c) a means to determine the currents present at each node based on changes in
the local and near neighbor state variable parameters as a result of
electrical
input data;
(d) a means to localize the origin the origin of the electrical stimulus;
(e) a means to determine the ecg at each node in the lattice;
(f) a means for modifying state variables as a result of electrical input
data;
(g) a means to use parallel computers to solve for the temporal evolution of
state
variables and electrical input data defined at each node of the lattice;
(h) a means to graphically display the changes at each node in the lattice as
a
result of said modification of the system of Claim 4 with electrical input
data;
(i) a means to store said state variables in computer memory and or external
data
storage devices for graphical display and analysis.
16. The model of claim 4 further comprising:
(a) a means for modifying one or more of said state variables as a result of
imposing chemical compound data in said state equations, wherein said
chemical compound data includes mechanism of action data; kinetic data;
biochemical, or mechanical data;
(b) a means for solving and updating the complete set of integral network
algorithms and chemical/composition data in time sets;
-34-
(c) a means to determine the currents present at each node based on changes in
the local and near neighbor state variable parameters as a result of
chemical/compound data;
(d) a means to determine the ecg at each node in the lattice;
(e) a means for modifying state variables as a result of chemical/compound
data;
(f) a means to use parallel computer to solve for the temporal evolution of
state
variables and chemical/compound data defined at each node in the lattice;
(g) a means to graphically display the changes at each node in the lattice as
a
result of said modification of the system with chemical compound data;
(h) a means to store said state variables in computer memory and or external
data
storage devices for graphical display and analysis.