EnergyMomentum - Maple Help

Home : Support : Online Help : Physics : EnergyMomentum

Physics[EnergyMomentum] - The EnergyMomentum tensor

 Calling Sequence EnergyMomentum[$\mathrm{\mu },\mathrm{\nu }$] EnergyMomentum[keyword]

Parameters

 mu, nu - the indices, as names representing integer numbers between 0 and the spacetime dimension, they can also be the numbers themselves keyword - optional, it can be definition, array or matrix, nonzero

Description

 • The EnergyMomentum[mu, nu], displayed as $\mathrm{T__μ,ν}$, is a computational representation for the EnergyMomentum tensor. The components of this tensor are initially set equal to zero. To define these components in any particular way use the Define command passing to it an equation with EnergyMomentum[mu, nu] on the left-hand side and the desired definition, as a symmetric tensorial expression or a symmetric matrix on the right-hand side.
 • When the metric is set to represent any curved spacetime, the EnergyMomentum tensor is the source of the gravitational field, entering Einstein's equations

${G}_{\mathrm{\mu },\mathrm{\nu }}=8\mathrm{\pi }{\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$

 where $\mathrm{G__μ,ν}$ is the Einstein tensor, expressed in terms of the Ricci tensor $\mathrm{R__μ,ν}$ as

${G}_{\mathrm{\mu },\mathrm{\nu }}={R}_{\mathrm{\mu },\mathrm{\nu }}-\frac{1}{2}{g}_{\mathrm{\mu },\mathrm{\nu }}{R}_{\mathrm{\alpha }}^{\mathrm{\alpha }}$

 • The EnergyMomentum tensor also satisfies ${▿}_{\mathrm{\mu }}\left({\mathrm{Τ}}_{}^{\mathrm{\mu },\mathrm{\nu }}\right)=0$, where ${▿}_{\mathrm{\mu }}$ is the covariant derivative operator D_, and this 4D divergence is entered as D_[mu](EnergyMomentum[mu,nu]).
 • When the indices of EnergyMomentum assume integer values they are expected to be between 0 and the spacetime dimension, prefixed by ~ when they are contravariant, and the corresponding value of EnergyMomentum is returned. The values 0 and 4, or for the case any dimension set for the spacetime, represent the same object. When the indices have symbolic values EnergyMomentum returns unevaluated after normalizing its indices taking into account their symmetry property.
 • Computations performed with the Physics package commands take into account EnergyMomentum's sum rule for repeated indices - see . and Simplify. The distinction between covariant and contravariant indices in the input of tensors is done by prefixing contravariant ones with ~, say as in ~mu; in the output, contravariant indices are displayed as superscripts. For contracted indices, you can enter them one covariant and one contravariant. Note however that - provided that the spacetime metric is galilean (Euclidean or Minkowski), or the object is a tensor also in curvilinear coordinates - this distinction in the input is not relevant, and so contracted indices can be entered as both covariant or both contravariant, in which case they will be automatically rewritten as one covariant and one contravariant. Tensors can have spacetime and space indices at the same time. To change the type of letter used to represent spacetime or space indices see Setup.
 • Besides being indexed with two indices, EnergyMomentum accepts three keywords:
 – definition: returns the definition of the EnergyMomentum tensor in terms of the Ricci tensor.
 – matrix: (synonyms: Matrix, array, Array, or no indices whatsoever, as in EnergyMomentum[]) returns a Matrix that when indexed with numerical values from 1 to the dimension of spacetime returns the value of each of the components of EnergyMomentum. If this keyword is passed together with indices, that can be covariant or contravariant, the resulting matrix takes into account the character of the indices.
 – nonzero: returns a set of equations, with the left-hand-side as a sequence of two positive numbers identifying the element of $\mathrm{G__μ,ν}$ and the corresponding value on the right-hand-side. Note that this set is actually the output of the ArrayElems command when passing to it the Array obtained with the keyword array.
 • Some automatic checking and normalization are carried out each time you enter EnergyMomentum[...]. The checking is concerned with possible syntax errors. The automatic normalization takes into account the symmetry of EnergyMomentum[mu,nu] with respect to interchanging the positions of the indices mu and nu.
 • The %EnergyMomentum command is the inert form of EnergyMomentum, so it represents the same mathematical operation but without performing it. To perform the operation, use value.

Examples

 > $\mathrm{with}\left(\mathrm{Physics}\right):$
 > $\mathrm{Setup}\left(\mathrm{mathematicalnotation}=\mathrm{true}\right)$
 $\left[{\mathrm{mathematicalnotation}}{=}{\mathrm{true}}\right]$ (1)
 > $\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]$
 ${{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}$ (2)

When Physics is loaded, all the components of ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$ are set equal to 0

 > $\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$
 ${{\mathrm{EnergyMomentum}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{rrrr}0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right]\right)$ (3)

In order to set these components in any particular way use the Define command passing to it an equation with EnergyMomentum[mu, nu] on the left-hand side and the desired definition, as a symmetric tensorial expression or a symmetric matrix, on the right-hand side. For example, the general form of ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$ in terms of the energy density $W$, the flux density ${S}_{j}$ and the stress tensor ${\mathrm{\sigma }}_{j,k}$ of a system can be entered as

 > $\mathrm{Setup}\left(\mathrm{spaceindices}=\mathrm{lowercaselatin_is}\right)$
 $\left[{\mathrm{spaceindices}}{=}{\mathrm{lowercaselatin_is}}\right]$ (4)
 > $\mathrm{Define}\left(S\left[j\right],\mathrm{\sigma }\left[j,k\right],\mathrm{symmetric}\right)$
 $\mathrm{Defined objects with tensor properties}$
 $\left\{{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{S}}_{{j}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\gamma }}}_{{i}{,}{j}}{,}{{\mathrm{σ}}}_{{j}{,}{k}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (5)
 > $\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{Matrix}\left(4,\left(\mathrm{\mu },\mathrm{\nu }\right)↦\mathbf{if}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{\mu }=4\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{if}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{\nu }=4\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}W\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{else}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}S\left[\mathrm{\nu }\right]\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{end}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{if}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{elif}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{\nu }=4\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}S\left[\mathrm{\mu }\right]\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{else}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathrm{\sigma }\left[\mathrm{\mu },\mathrm{\nu }\right]\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{end}\phantom{\rule[-0.0ex]{0.3em}{0.0ex}}\mathbf{if}\right)$
 ${{\mathrm{EnergyMomentum}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{cccc}{\mathrm{σ}}_{1,1}& {\mathrm{σ}}_{1,2}& {\mathrm{σ}}_{1,3}& {S}_{1}\\ {\mathrm{σ}}_{2,1}& {\mathrm{σ}}_{2,2}& {\mathrm{σ}}_{2,3}& {S}_{2}\\ {\mathrm{σ}}_{3,1}& {\mathrm{σ}}_{3,2}& {\mathrm{σ}}_{3,3}& {S}_{3}\\ {S}_{1}& {S}_{2}& {S}_{3}& W\end{array}\right]\right)$ (6)

You can now set these to be the components of ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$

 > $\mathrm{Define}\left(\right)$
 $\mathrm{Defined objects with tensor properties}$
 $\left\{{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{S}}_{{j}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\gamma }}}_{{i}{,}{j}}{,}{{\mathrm{σ}}}_{{j}{,}{k}}{,}{{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}\right\}$ (7)

After this definition, you can query about the definition, the nonzero components, any particular covariant or contravariant component via

 > $\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$
 ${{\mathrm{EnergyMomentum}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{cccc}{\mathrm{σ}}_{1,1}& {\mathrm{σ}}_{1,2}& {\mathrm{σ}}_{1,3}& {S}_{1}\\ {\mathrm{σ}}_{2,1}& {\mathrm{σ}}_{2,2}& {\mathrm{σ}}_{2,3}& {S}_{2}\\ {\mathrm{σ}}_{3,1}& {\mathrm{σ}}_{3,2}& {\mathrm{σ}}_{3,3}& {S}_{3}\\ {S}_{1}& {S}_{2}& {S}_{3}& W\end{array}\right]\right)$ (8)
 > $\mathrm{EnergyMomentum}\left[1,2\right]$
 ${{\mathrm{σ}}}_{{1}{,}{2}}$ (9)

In the definition above, all the components of ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$ are constant. To set part or all of them as depending on the coordinates, for instance in a generic coordinate system and using spherical coordinates,

 > $\mathrm{Setup}\left(\mathrm{coordinates}=\mathrm{spherical}\right)$
 ${}{}\mathrm{Systems of spacetime coordinates are:}{}{}{}\left\{X=\left(r{,}\mathrm{\theta }{,}\mathrm{\phi }{,}t\right)\right\}$
 $\mathrm{_______________________________________________________}$
 $\left[{\mathrm{coordinatesystems}}{=}\left\{{X}\right\}\right]$ (10)

you can indicate the functionality in the definition. For example, set $W$ to be constant (i.e. no functionality) but the flux density ${S}_{j}$ and stress ${\mathrm{sigma}}_{j,k}$ tensors depending on $X$

 > $\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{Matrix}\left(4,\left(\mathrm{μ},\mathrm{ν}\right)→\mathbf{if}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathrm{μ}=4\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{if}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathrm{ν}=4\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}W\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{else}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{S}_{\mathrm{ν}}\left(X\right)\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{end if}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{elif}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathrm{ν}=4\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{then}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{S}_{\mathrm{μ}}\left(X\right)\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{else}\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}{\mathrm{σ}}_{\mathrm{μ},\mathrm{ν}}\left(X\right)\phantom{\rule[-0.0ex]{0.5em}{0.0ex}}\mathbf{end if}\right)$
 ${{\mathrm{EnergyMomentum}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{cccc}{\mathrm{σ}}_{1,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\\ {\mathrm{σ}}_{2,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{2,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{2,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\\ {\mathrm{σ}}_{3,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{3,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {\mathrm{σ}}_{3,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\\ {S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& {S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)& W\end{array}\right]\right)$ (11)
 > $\mathrm{CompactDisplay}\left(\right)$
 ${S}{}\left({X}\right){}{\mathrm{will now be displayed as}}{}{S}$
 ${\mathrm{σ}}{}\left({X}\right){}{\mathrm{will now be displayed as}}{}{\mathrm{σ}}$ (12)

Define now ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$ with these components

 > $\mathrm{Define}\left(\right)$
 $\mathrm{Defined objects with tensor properties}$
 $\left\{{{\mathrm{\gamma }}}_{{\mathrm{\mu }}}{,}{{\mathrm{\sigma }}}_{{\mathrm{\mu }}}{,}{{S}}_{{j}}{,}{{\partial }}_{{\mathrm{\mu }}}{,}{{g}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\gamma }}}_{{i}{,}{j}}{,}{{\mathrm{σ}}}_{{j}{,}{k}}{,}{{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{\mathrm{\epsilon }}}_{{\mathrm{\alpha }}{,}{\mathrm{\beta }}{,}{\mathrm{\mu }}{,}{\mathrm{\nu }}}{,}{{X}}_{{\mathrm{\mu }}}\right\}$ (13)

To see the continuity equations for the components of ${\mathrm{Τ}}_{\mathrm{\mu },\mathrm{\nu }}$, use for instance the inert version of the covariant derivative operator D_ and the TensorArray command

 > $\left(\mathrm{%D_}\left[\mathrm{\mu }\right]=\mathrm{D_}\left[\mathrm{\mu }\right]\right)\left(\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]\right)$
 ${{\mathrm{%D_}}}_{{\mathrm{μ}}}{}\left({{\mathrm{EnergyMomentum}}}_{{\mathrm{~mu}}{,}{\mathrm{ν}}}\right){=}{0}$ (14)
 > $\mathrm{Vector}\left[\mathrm{column}\right]\left(\mathrm{TensorArray}\left(\right)\right)$
 $\left[\begin{array}{c}{\mathrm{%d_}}_{1}{}\left(-{\mathrm{σ}}_{1,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{2}{}\left(-{\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{3}{}\left(-{\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{4}{}\left({S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)=0\\ {\mathrm{%d_}}_{1}{}\left(-{\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{2}{}\left(-{\mathrm{σ}}_{2,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{3}{}\left(-{\mathrm{σ}}_{2,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{4}{}\left({S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)=0\\ {\mathrm{%d_}}_{1}{}\left(-{\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{2}{}\left(-{\mathrm{σ}}_{2,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{3}{}\left(-{\mathrm{σ}}_{3,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{4}{}\left({S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)=0\\ {\mathrm{%d_}}_{1}{}\left(-{S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{2}{}\left(-{S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{3}{}\left(-{S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right),\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)+{\mathrm{%d_}}_{4}{}\left(W,\left[r,\mathrm{θ},\mathrm{φ},t\right]\right)=0\end{array}\right]$ (15)
 > $\mathrm{value}\left(\right)$
 $\left[\begin{array}{c}-\left(\frac{\partial }{\partial r}{}{\mathrm{σ}}_{1,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{θ}}{}{\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{φ}}{}{\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)+\frac{\partial }{\partial t}{}{S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0\\ -\left(\frac{\partial }{\partial r}{}{\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{θ}}{}{\mathrm{σ}}_{2,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{φ}}{}{\mathrm{σ}}_{2,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)+\frac{\partial }{\partial t}{}{S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0\\ -\left(\frac{\partial }{\partial r}{}{\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{θ}}{}{\mathrm{σ}}_{2,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{φ}}{}{\mathrm{σ}}_{3,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)+\frac{\partial }{\partial t}{}{S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0\\ -\left(\frac{\partial }{\partial r}{}{S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{θ}}{}{S}_{2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)-\left(\frac{\partial }{\partial \mathrm{φ}}{}{S}_{3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)\right)=0\end{array}\right]$ (16)

In curved spaces, EnergyMomentum enters Einstein's equations as the source of the gravitational field

 > $\mathrm{g_}\left[\mathrm{sc}\right]$
 ${}{}\mathrm{_______________________________________________________}$
 ${}{}\mathrm{The Schwarzschild metric in coordinates}{}{}{}\left[r{,}\mathrm{\theta }{,}\mathrm{\phi }{,}t\right]$
 $\mathrm{Parameters:}{}\left[m\right]$
 $\mathrm{Signature:}{}\left(\mathrm{- - - +}\right)$
 ${}{}\mathrm{_______________________________________________________}$
 ${{\mathrm{g_}}}_{{\mathrm{μ}}{,}{\mathrm{ν}}}{=}\left(\left[\begin{array}{cccc}\frac{r}{2{}m-r}& 0& 0& 0\\ 0& -{r}^{2}& 0& 0\\ 0& 0& -{r}^{2}{}{\mathrm{sin}{}\left(\mathrm{θ}\right)}^{2}& 0\\ 0& 0& 0& \frac{r-2{}m}{r}\end{array}\right]\right)$ (17)
 > $\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$
 ${{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}\frac{{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}}{{8}{}{\mathrm{\pi }}}{,}{{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}\left[\begin{array}{cccc}{{\mathrm{σ}}}_{{1}{,}{1}}{}\left({X}\right)& {{\mathrm{σ}}}_{{1}{,}{2}}{}\left({X}\right)& {{\mathrm{σ}}}_{{1}{,}{3}}{}\left({X}\right)& {{S}}_{{1}}{}\left({X}\right)\\ {{\mathrm{σ}}}_{{2}{,}{1}}{}\left({X}\right)& {{\mathrm{σ}}}_{{2}{,}{2}}{}\left({X}\right)& {{\mathrm{σ}}}_{{2}{,}{3}}{}\left({X}\right)& {{S}}_{{2}}{}\left({X}\right)\\ {{\mathrm{σ}}}_{{3}{,}{1}}{}\left({X}\right)& {{\mathrm{σ}}}_{{3}{,}{2}}{}\left({X}\right)& {{\mathrm{σ}}}_{{3}{,}{3}}{}\left({X}\right)& {{S}}_{{3}}{}\left({X}\right)\\ {{S}}_{{1}}{}\left({X}\right)& {{S}}_{{2}}{}\left({X}\right)& {{S}}_{{3}}{}\left({X}\right)& {W}\end{array}\right]$ (18)

Take the first of these two equations and compute a tensor array for it

 > $\left[1\right]$
 ${{\mathrm{Τ}}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}{=}\frac{{{G}}_{{\mathrm{\mu }}{,}{\mathrm{\nu }}}}{{8}{}{\mathrm{\pi }}}$ (19)
 > $\mathrm{TensorArray}\left(\left[1\right]\right)$
 $\left[\begin{array}{cccc}{\mathrm{σ}}_{1,1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0& {\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0& {\mathrm{σ}}_{1,3}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0& {S}_{1}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0\\ {\mathrm{σ}}_{1,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0& {\mathrm{σ}}_{2,2}{}\left(r,\mathrm{θ},\mathrm{φ},t\right)=0& {\mathrm{σ}}_{2,3}{}\left(r\right)\end{array}\right]$