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

where $\mathrm{G\_\_\μ,\ν}$ is the Einstein tensor, expressed in terms of the Ricci tensor $\mathrm{R\_\_\μ,\ν}$ as

The EnergyMomentum tensor also satisfies ${▿}_{\mathrm{\mu }}\left({\mathrm{{\rm T}}}_{}^{\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.

$\mathrm{with}\left(\mathrm{Physics}\right)\:$

$\mathrm{Setup}\left(\mathrm{mathematicalnotation}=\mathrm{true}\right)$

$\left[\mathrm{mathematicalnotation}=\mathrm{true}\right]$

$\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}$

$\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\left[\begin{array}{cccc}0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& 0& 0\end{array}\right]$

$\mathrm{Setup}\left(\mathrm{spaceindices}=\mathrm{lowercaselatin\_is}\right)$

$\left[\mathrm{spaceindices}=\mathrm{lowercaselatin\_is}\right]$

$\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\}$

$\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{Matrix}\left(4\,\left(\mathrm{\mu }\,\mathrm{\nu }\right)\mapsto \mathbf{if}\mathrm{\mu }=4\mathbf{then}\mathbf{if}\mathrm{\nu }=4\mathbf{then}W\mathbf{else}S\left[\mathrm{\nu }\right]\mathbf{end}\mathbf{if}\mathbf{elif}\mathrm{\nu }=4\mathbf{then}S\left[\mathrm{\mu }\right]\mathbf{else}\mathrm{\sigma }\left[\mathrm{\mu },\mathrm{\nu }\right]\mathbf{end}\mathbf{if}\right)$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\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]$

$\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{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\right\}$

$\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\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]$

$\mathrm{EnergyMomentum}\left[1,2\right]$

${\mathrm{\σ}}_{1,2}$

$\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]$

$\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{Matrix}\left(4\,\left(\mathrm{\mu }\,\mathrm{\nu }\right)\mapsto \mathbf{if}\mathrm{\mu }=4\mathbf{then}\mathbf{if}\mathrm{\nu }=4\mathbf{then}W\mathbf{else}S\left[\mathrm{\nu }\right]\left(X\right)\mathbf{end}\mathbf{if}\mathbf{elif}\mathrm{\nu }=4\mathbf{then}S\left[\mathrm{\mu }\right]\left(X\right)\mathbf{else}\mathrm{\sigma }\left[\mathrm{\mu },\mathrm{\nu }\right]\left(X\right)\mathbf{end}\mathbf{if}\right)$

${\mathrm{{\rm T}}}_{\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]$

$\mathrm{CompactDisplay}\left(\right)$

$S\left(r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right)\mathrm{will\; now\; be\; displayed\; as}S$

$\mathrm{\σ}\left(r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right)\mathrm{will\; now\; be\; displayed\; as}\mathrm{\σ}$

$\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{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\,X_{\mathrm{\mu }}\right\}$

$\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{\mu }}\left({\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}\right)=0$

$\mathrm{Vector}\left[\mathrm{column}\right]\left(\mathrm{TensorArray}\left(\right)\right)$

$\left[\begin{array}{c}{\partial }_1\left(-{\mathrm{\σ}}_{1,1}\right)+{\partial }_2\left(-{\mathrm{\σ}}_{1,2}\right)+{\partial }_3\left(-{\mathrm{\σ}}_{1,3}\right)+{\partial }_4\left(S_1\right)=0\\ {\partial }_1\left(-{\mathrm{\σ}}_{1,2}\right)+{\partial }_2\left(-{\mathrm{\σ}}_{2,2}\right)+{\partial }_3\left(-{\mathrm{\σ}}_{2,3}\right)+{\partial }_4\left(S_2\right)=0\\ {\partial }_1\left(-{\mathrm{\σ}}_{1,3}\right)+{\partial }_2\left(-{\mathrm{\σ}}_{2,3}\right)+{\partial }_3\left(-{\mathrm{\σ}}_{3,3}\right)+{\partial }_4\left(S_3\right)=0\\ {\partial }_1\left(-S_1\right)+{\partial }_2\left(-S_2\right)+{\partial }_3\left(-S_3\right)+{\partial }_4\left(W\right)=0\end{array}\right]$

$\mathrm{value}\left(\right)$

$\left[\begin{array}{c}-{\left({\mathrm{\σ}}_{1,1}\right)}_r-{\left({\mathrm{\σ}}_{1,2}\right)}_{\mathrm{\theta }}-{\left({\mathrm{\σ}}_{1,3}\right)}_{\mathrm{\phi }}+{\stackrel{\mathbf{.}}{S}}_1=0\\ -{\left({\mathrm{\σ}}_{1,2}\right)}_r-{\left({\mathrm{\σ}}_{2,2}\right)}_{\mathrm{\theta }}-{\left({\mathrm{\σ}}_{2,3}\right)}_{\mathrm{\phi }}+{\stackrel{\mathbf{.}}{S}}_2=0\\ -{\left({\mathrm{\σ}}_{1,3}\right)}_r-{\left({\mathrm{\σ}}_{2,3}\right)}_{\mathrm{\theta }}-{\left({\mathrm{\σ}}_{3,3}\right)}_{\mathrm{\phi }}+{\stackrel{\mathbf{.}}{S}}_3=0\\ -{\left(S_1\right)}_r-{\left(S_2\right)}_{\mathrm{\theta }}-{\left(S_3\right)}_{\mathrm{\phi }}=0\end{array}\right]$

$\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{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left[\begin{array}{cccc}\frac{r}{2m-r}& 0& 0& 0\\ 0& -r^2& 0& 0\\ 0& 0& -r^2{\sin \left(\mathrm{\theta }\right)}^2& 0\\ 0& 0& 0& \frac{r-2m}{r}\end{array}\right]$

$\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\frac{G_{\mathrm{\mu },\mathrm{\nu }}}{8\mathrm{\pi }},{\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\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]$

$\left[1\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\frac{G_{\mathrm{\mu },\mathrm{\nu }}}{8\mathrm{\pi }}$

$\mathrm{TensorArray}\left(\left[1\right]\right)$

$\left[\begin{array}{cccc}{\mathrm{\σ}}_{1,1}=0& {\mathrm{\σ}}_{1,2}=0& {\mathrm{\σ}}_{1,3}=0& S_1=0\\ {\mathrm{\σ}}_{1,2}=0& {\mathrm{\σ}}_{2,2}=0& {\mathrm{\σ}}_{2,3}=0& S_2=0\\ {\mathrm{\σ}}_{1,3}=0& {\mathrm{\σ}}_{2,3}=0& {\mathrm{\σ}}_{3,3}=0& S_3=0\\ S_1=0& S_2=0& S_3=0& W=0\end{array}\right]$

$\mathrm{g\_}\left[\mathrm{tol}\right]$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$\mathrm{The\; Tolman\; metric\; in\; coordinates}\left[r\,\mathrm{\theta }\,\mathrm{\phi }\,t\right]$

$\mathrm{Parameters:}\left[R\left(t\,r\right)\,E\left(r\right)\right]$

$\mathrm{Signature:}\left(\mathrm{-\; -\; -\; +}\right)$

$\mathrm{\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_\_}$

$g_{\mathrm{\mu },\mathrm{\nu }}=\left[\begin{array}{cccc}-\frac{{R_r}^2}{1+2E\left(r\right)}& 0& 0& 0\\ 0& -{R\left(t\,r\right)}^2& 0& 0\\ 0& 0& -{R\left(t\,r\right)}^2{\sin \left(\mathrm{\theta }\right)}^2& 0\\ 0& 0& 0& 1\end{array}\right]$

$\mathrm{CompactDisplay}\left(\right)$

$E\left(r\right)\mathrm{will\; now\; be\; displayed\; as}E$

$R\left(t\,r\right)\mathrm{will\; now\; be\; displayed\; as}R$

$\mathrm{TensorArray}\left(\left[1\right]\right)$

$\left[\begin{array}{cccc}{\mathrm{\σ}}_{1,1}=-\frac{{R_r}^2\left(2R\stackrel{\mathbf{..}}{R}+{{\stackrel{\mathbf{.}}{R}}^{}}^2-2E\right)}{8R^2\left(1+2E\right)\mathrm{\pi }}& {\mathrm{\σ}}_{1,2}=0& {\mathrm{\σ}}_{1,3}=0& S_1=0\\ {\mathrm{\σ}}_{1,2}=0& {\mathrm{\σ}}_{2,2}=-\frac{R\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& {\mathrm{\σ}}_{2,3}=0& S_2=0\\ {\mathrm{\σ}}_{1,3}=0& {\mathrm{\σ}}_{2,3}=0& {\mathrm{\σ}}_{3,3}=-\frac{R{\sin \left(\mathrm{\theta }\right)}^2\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& S_3=0\\ S_1=0& S_2=0& S_3=0& W=\frac{{{\stackrel{\mathbf{.}}{R}}^{}}^2{R}_r+2\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)R-2E_{r}R-2R_{r}E}{8R_r{R}^2\mathrm{\pi }}\end{array}\right]$

$\frac{\mathrm{Einstein}\left[i,j\right]}{8\mathrm{\pi }}$

$\frac{G_{i,j}}{8\mathrm{\pi }}$

$\mathrm{\sigma }\left[j,k\right]=\mathrm{TensorArray}\left(\right)$

${\mathrm{\σ}}_{j,k}=\left[\begin{array}{ccc}-\frac{{R_r}^2\left(2R\stackrel{\mathbf{..}}{R}+{{\stackrel{\mathbf{.}}{R}}^{}}^2-2E\right)}{8R^2\left(1+2E\right)\mathrm{\pi }}& 0& 0\\ 0& -\frac{R\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& 0\\ 0& 0& -\frac{R{\sin \left(\mathrm{\theta }\right)}^2\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}\end{array}\right]$

$\mathrm{Define}\left(\right)$

$\mathrm{Defined\; objects\; with\; tensor\; properties}$

$\left\{{▿}_{\mathrm{\mu }}\,{\mathrm{\gamma }}_{\mathrm{\mu }}\,{\mathrm{\sigma }}_{\mathrm{\mu }}\,R_{\mathrm{\mu },\mathrm{\nu }}\,R_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,S_j\,C_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\partial }_{\mathrm{\mu }}\,g_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\gamma }}_{i,j}\,{\mathrm{\σ}}_{j,k}\,{\mathrm{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,G_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\,X_{\mathrm{\mu }}\right\}$

$\mathrm{\sigma }\left[\mathrm{definition}\right]$

${\mathrm{\σ}}_{j,k}=\left[\begin{array}{ccc}-\frac{{R_r}^2\left(2R\stackrel{\mathbf{..}}{R}+{{\stackrel{\mathbf{.}}{R}}^{}}^2-2E\right)}{8R^2\left(1+2E\right)\mathrm{\pi }}& 0& 0\\ 0& -\frac{R\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& 0\\ 0& 0& -\frac{R{\sin \left(\mathrm{\theta }\right)}^2\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}\end{array}\right]$

$\mathrm{\sigma }\left[\mathrm{`~`}\right]$

${\mathrm{\σ}}_{\phantom{}\phantom{i,j}}^{\phantom{}i,j}=\left[\begin{array}{ccc}-\frac{\left(2R\stackrel{\mathbf{..}}{R}+{{\stackrel{\mathbf{.}}{R}}^{}}^2-2E\right)\left(1+2E\right)}{8R^2\mathrm{\pi }{R_r}^2}& 0& 0\\ 0& \frac{-\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R-\stackrel{\mathbf{..}}{R}{R}_r-\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)+E_r}{8R^3{R}_r\mathrm{\pi }}& 0\\ 0& 0& -\frac{{\csc \left(\mathrm{\theta }\right)}^2\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R^3{R}_r\mathrm{\pi }}\end{array}\right]$

$\mathrm{Define}\left(\mathrm{EnergyMomentum}\left[\mathrm{\mu },\mathrm{\nu }\right]=\mathrm{Matrix}\left(4\,\left(\mathrm{\mu }\,\mathrm{\nu }\right)\mapsto \mathbf{if}\mathrm{\mu }=4\mathbf{then}\mathbf{if}\mathrm{\nu }=4\mathbf{then}W\mathbf{else}S\left[\mathrm{\nu }\right]\left(X\right)\mathbf{end}\mathbf{if}\mathbf{elif}\mathrm{\nu }=4\mathbf{then}S\left[\mathrm{\mu }\right]\left(X\right)\mathbf{else}\mathrm{\sigma }\left[\mathrm{\mu },\mathrm{\nu }\right]\mathbf{end}\mathbf{if}\right)\right)$

$\mathrm{Defined\; objects\; with\; tensor\; properties}$

$\left\{{▿}_{\mathrm{\mu }}\,{\mathrm{\gamma }}_{\mathrm{\mu }}\,{\mathrm{\sigma }}_{\mathrm{\mu }}\,R_{\mathrm{\mu },\mathrm{\nu }}\,R_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,S_j\,C_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha },\mathrm{\beta }}\,{\partial }_{\mathrm{\mu }}\,g_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\gamma }}_{i,j}\,{\mathrm{\σ}}_{j,k}\,{\mathrm{\Gamma }}_{\mathrm{\mu },\mathrm{\nu },\mathrm{\alpha }}\,G_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}\,{\mathrm{\epsilon }}_{\mathrm{\alpha },\mathrm{\beta },\mathrm{\mu },\mathrm{\nu }}\,X_{\mathrm{\mu }}\right\}$

$\mathrm{EnergyMomentum}\left[\mathrm{definition}\right]$

${\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\frac{G_{\mathrm{\mu },\mathrm{\nu }}}{8\mathrm{\pi }},{\mathrm{{\rm T}}}_{\mathrm{\mu },\mathrm{\nu }}=\left[\begin{array}{cccc}-\frac{{R_r}^2\left(2R\stackrel{\mathbf{..}}{R}+{{\stackrel{\mathbf{.}}{R}}^{}}^2-2E\right)}{8R^2\left(1+2E\right)\mathrm{\pi }}& 0& 0& S_1\\ 0& -\frac{R\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& 0& S_2\\ 0& 0& -\frac{R{\sin \left(\mathrm{\theta }\right)}^2\left(\left(\frac{{\partial }^{3}R}{\partial t^2\partial r}\right)R+\stackrel{\mathbf{..}}{R}{R}_r+\stackrel{\mathbf{.}}{R}\left(\frac{{\partial }^{2}R}{\partial t\partial r}\right)-E_r\right)}{8R_r\mathrm{\pi }}& S_3\\ S_1& S_2& S_3& W\end{array}\right]$

[1] Landau, L.D., and Lifshitz, E.M. The Classical Theory of Fields, Course of Theoretical Physics Volume 2, fourth revised English edition. Elsevier, 1975.