ecterrab

14085 Reputation

24 Badges

19 years, 121 days

MaplePrimes Activity


These are replies submitted by ecterrab

@nm 

A dot release is not a version ahead. You need to be clear about this. Only a few bug fixes enter the dot releases, fixes considered critical for one reason or another. The Physics Updates contain significantly more fixes than that, also developments (entirely absent from dot releases). So no, you do not have Maple 2 versions ahead.

About your comment regarding "... fixes for very old bugs", what is what you consider old? The fixes mentioned in this conversation are for bugs all of them newer than March 2024, this year.

Your statement about odetest hanging does not change what I said before. In that "odetest hangs" post of yours, first, you show it does not hang. Then you show that if you first interrupt computations using timelimit and assumptions, next a hang occurrs. Perhaps you are unaware that interrupting a computation in that way leaves fragments of results in caches. That does not usually generate problems, but sometimes it does. In particular, your problem in that post involves is, assume and coulditbe. As you probably remember, in some cases, is and coulditbe interrupt with an error the first time, then they do not. These nuances can produce problems as the one you mentioned, not a problem in odetest.

By the way, if you want to know what hangs, Maple offers the tracelast instruction - I already mentioned this to you before. That allows you to post the actual problem, which in your example is not a problem in odetest; if it were, then it would have hanged the first time. It is also not a problem in the Maplesoft Physics Updates, even if its presence uncovers a problem not noticed before - that sometimes happens, you know that too.

Then the way you wrote this post I am commenting, it does look like "why are the fixes found in the Physics Updates not present when I remove the Physics Updates?" For me, that indicates you are not clear about how this works. Maybe you were not clear of what a dot release is.

Anyway, I can see you expecting even more from the software; I don't share your expectation, the comments aim at explaining that.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

The Maplesoft Physics Updates is a project that aims at bringing to the current Maple release fixes that are present in the version under development, as well as new developments in the areas of Physics, Differential Equations and Mathematical Functions. Somehow, your posts these days, @nm , either indicate you are not clear about that or perhaps you are trying to say something else?

In this one above, you say that a version of the Maplesoft Physics Updates fixed a reported problem (and, you know, also other problems not reported in Mapleprimes). Well, if you remove the Maplesoft Physics Updates, the problems that it fixes are not fixed anymore. I find that natural—what were you expecting?

I also mentioned other times that installing the Physics Updates may reveal a problem that was not apparent without it. Yes, that can happen—report the case. Depending on the problem, it may get fixed in one of the next versions of the Updates, or you may need to wait until the next release because, sometimes, other problems are seen as a higher priority. All this is known.

In summary, you have two options: install this Maplesoft Physics Updates package and have all those fixes and developments "in advance", or wait for them until the next Maple dot or full release (for Maple 2024, that would be next year).

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Hullzie16 
Nice reply. Attached are a few comments intercalated in your worksheet that may be of interest. Physics is an environment with a myriad of simplifications for the input. They are not necessary but make entering things easier.

Tensor_Respond_(comments).mw

Best!
Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

This ODE system is solvable with the existing dsolve. But there is a problem in simplify subroutines. Before calling dsolve, input debug(`simplify/RootOf`, statements = false, depth = 1) and you see, at the end, the line that causes `simplify/RootOf` to interrupt with the Error message you saw. It shouldn't happen. I will see to include a fix for this in the next or so version of Maplesoft Physics Updates for Maple 2024, and write again here.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@mehran rajabi 
In that case you can use PDEtools:-dcoeffs. What confused me and I imagine others is that in your question you said "[Given the equation] ... I want to calculate separately, c1,c2,c3,c4,c5" (in a computer algebra context, that has the meaning of performing some calculation) while, actually, you do not want to calculate anything, only extract the coefficients.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

Hi @Carl Love, there are various, perhaps more like "many", things about the ODE and PDE code, the FunctionAdvisor, conversion network for math funcs, Physics ... that, while documented, the details are so many that end up buried, not noticed. The syntax you refer is explained in ?dsolve,system; scan for 'ranking', it appears almost on top of the page and also is explained in its own subsection, "The unknowns of the problem and rankings".

My impression of this thread is that the question was not properly formulated, or I may have understood it wrongly ... - I was in doubt. Mainly: my interpretation of the question is that this was about solving a system with f(x) and 5 c(x) "functions"whose derivatives are equal to 0, and the solution desired involved "f as a function of the c, and the c as functions of themselves at most", and that is what this ranking represents.

So, if I am correct, what was missing was a more precise description of what was intended. After that, the formulation followed naturally (for whoever knows the functionality available in 'ranking' details).

Best

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@nm 
Thank you. I added as a PS the worksheet containing the tracelast output, which is what allowed for reproducing the problem and fix it. And thanks to @aroche for his quick fix. The fix is already uploaded in the latest Maplesoft Physics Updates for Maple 2024, so to activate the fix, as usual, open Maple and input Physics:-Version(latest).

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@nm 

Generally speaking, the more extensive is your post, the most unreadable and time-consuming. I suggest you to try to make your reports more compact. By the way, I work with GUI. Always try GUI first. If I cannot reproduce I move into 1D Math input, then TTY.

Now on the issue: there is a bug in the subroutines of `simplify/common/factors`. I passed the problem to the person who wrote the routine, I suppose we will have a fix probably by the end of tomorrow or Friday. Regarding your comments about whether you display or not things ... I don't know. The bug - restricted to the input to the mentioned routine, however, is reproducible no matter what you display.

Best

Thanks Austin, for the so quick fix; it is uploaded now within the Maplesoft Physics Updates v.1793. To install: open Maple and enter Physics:-Version(latest)

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@C_R 
It is too time-consuming to answer about the details, but here are more comments. I think your approach, "based on observations", is the right approach, always. The page ?dsolve,setup is relevant only for those who want to have a glimpse on how dsolve is coded, its logic, and perhaps change the way it works. It is not a page of interest for many people. Its existence, perhaps, is more related to my academic origin/teaching: do not present black boxes, and never fall for "I know better". For the same rasons you have the PDEtools:-Library and its help page. Regarding the formulation known as d'Alembert's method, this method is present there, under `dsolve/methods`[1, high_degree]. The page on dsolve's setup could perhaps say "everything starting with 1 in eq.(1) of the page is about methods for 1st order ODEs".

Regarding PDEtools: I wrote this package before moving to Canada, then immersed myself into coding a new ODE solver and the symmetry methods for ODEs, followed by extending PDEtools with symmetry methods for PDEs, finally rewriting the DifferentialAlgebra package. I suppose that due to this historical sequence these two packages PDEtools and ODEtools (this one exists only on background, invisible) are strongly intertwined, and connected with the three differential elimination packages mentioned in PDEtools:-casesplit.

In the end, once you realize that from a solving-logic point of view all these things are differential equations of different orders and number of independent variables, the distinction between ODEs and PDEs (or even non-differential equations) is of little relevance: the majority of the tools for PDEs are of use for ODEs. And some commands are indeed present in both DEtools (the visible package for ODEs, not PDEs) and PDEtools, like casesplit, dpolyform, diff_table, dsubs, ... More important: the subroutines presented in ?PDEtools,Library, are used all around in the ODE code.

Regarding your last question, there is a help page, ?dsolve,education, that presents most things, though indeed not an example as the one discussed in this post (branches and related stuff). It may be good to add such an example. Still, the documentation for differential equations is already so large and detailed that adding more may help in some way but, on the other hand, adds clutter. In part, for that reason, I wrote the odeadvisor and its help pages together with dsolve - trying to provide a more structured way to peek at what dsolve does. Anyway, I will think about one more example.

Finally, I share your impression that "studying ODEs with a textbook" does not make sense. In my case, I learned about symmetry methods to then discover the standard classifications shown in textbooks (mainly Kamke's and Zwillinger's) and note that they were only a minimal subset of particular cases of symmetries or integrating factors.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@vv 
Yes. It is adjusted in v.1786 of the Physics Updates.

@vv 
You are right. And although the FunctionAdvisor cannot list all possible formulae, these ones, related to Mittag-Leffler's theorem, actually not only for tan, are relevant; Wolinski is also right. I missed all these when coding the FunctionAdvisor. For tan, the issue is resolved installing the Maplesoft Physics Updates v.1875 or newer, just uploaded. As usual, to install, open Maple and input Physics:-Version(latest);

The change includes a change in convert(tan(z), Sum), and this is a pic of the FunctionAdvisor output (btw, originally called FunctionWizard .. I liked that name more :)

 

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@dharr 
This Warning appears in the updated Maple 2024.1, as per this image, all in TTY with no initialization and only the libraries of the main installation:

I know from where is this "Warning" coming (the DifferentialThomas elimination package); will take a closer look later today.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@KM4@Josecherukara

 

This worksheet is the one I used to answer before, cutting after introducing the Wikipedia formula for the linearized form of the Ricci tensor (see Wikipedia for Linearized gravity), and introducing a new subsection showing a couple of steps deriving the linearized forms with the help of a new Physics:-Library:-Linearize command. Note that to reproduce this material you need to install the latest Maplesoft Physics updates for your copy of Maple 2024

Physics:-Version()

`The "Physics Updates" version in the MapleCloud is 1782 and is the same as the version installed in this computer, created 2024, August 8, 16:5 hours Pacific Time.`

(1)

 

restart

with(Physics)

 

Set cartesian coordinates (you could use spherical, cylindrical, or define your own ) and automatic simplification

Setup(coordinates = cartesian, simplification = true)

`* Partial match of '`*simplification*`' against keyword '`*automaticsimplification*`' `

 

`Systems of spacetime coordinates are:`*{X = (x, y, z, t)}

 

_______________________________________________________

 

[automaticsimplification = true, coordinatesystems = {X}]

(2)

Suppose you want to define a perturbation around the Minkowski metric (the steps are the same for any other metric)

g_[]

Physics:-g_[mu, nu] = Matrix(%id = 36893488152329516324)

(3)

Define a perturbation h[mu, nu], that in the general case depends on the coordinates and s not diagonal; the only requirement is that it is symmetric (to have it diagonal, change symmetric by diagonal)

h[mu, nu] = Matrix(4, proc (i, j) options operator, arrow; delta[i, j](X) end proc, shape = symmetric)

h[mu, nu] = Matrix(%id = 36893488152329494764)

(4)

"Define(?)"

`Defined objects with tensor properties`

 

{Physics:-Dgamma[mu], Physics:-Psigma[mu], Physics:-d_[mu], Physics:-g_[mu, nu], h[mu, nu], Physics:-LeviCivita[alpha, beta, mu, nu], Physics:-SpaceTimeVector[mu](X)}

(5)

If you are interested in working the details, you can define now the "new metric" as a perturbation over the existing metric uncommenting these two lines, but to follow the Wikipedia page do not redefine the metric.


Define a tensor representing the Ricci tensor according to the Wikipedia page (it is possible to derive this formula starting from the definition of Ricci in terms of the metric)

%Ricci[mu, nu] = 1/2*(d_[sigma](d_[mu](h[`~sigma`, nu]))+d_[sigma](d_[nu](h[`~sigma`, mu]))-d_[mu](d_[nu](h[`~sigma`, sigma]))-dAlembertian(h[mu, nu]))

%Ricci[mu, nu] = (1/2)*Physics:-d_[mu](Physics:-d_[sigma](h[`~sigma`, nu], [X]), [X])+(1/2)*Physics:-d_[nu](Physics:-d_[sigma](h[`~sigma`, mu], [X]), [X])-(1/2)*Physics:-d_[mu](Physics:-d_[nu](h[`~sigma`, sigma], [X]), [X])-(1/2)*Physics:-dAlembertian(h[mu, nu], [X])

(6)

Deriving the linearized form of the Ricci and Einstein tensors

 

Here I change the approach of the previous worksheet and this time redefine the metric in order to derive the linearized form of the Ricci or any other tensor. To make things more transparent, let's define first a tensor equal to the Minkowski metric

"eta[mu,nu] = rhs(?)"

eta[mu, nu] = Matrix(%id = 36893488152306604380)

(7)

"Define(?)"

`Defined objects with tensor properties`

 

{Physics:-Dgamma[mu], Physics:-Psigma[mu], Physics:-d_[mu], eta[mu, nu], Physics:-g_[mu, nu], h[mu, nu], Physics:-LeviCivita[alpha, beta, mu, nu], Physics:-SpaceTimeVector[mu](X)}

(8)

We will assume as usual that `h__μ,ν` is a small perturbation of `η__μ,ν` so that the metric is

g_[mu, nu] = eta[mu, nu]+h[mu, nu]

Physics:-g_[mu, nu] = eta[mu, nu]+h[mu, nu]

(9)

And this time we define the metric using this equation above

Define(Physics[g_][mu, nu] = eta[mu, nu]+h[mu, nu])

_______________________________________________________

 

`Coordinates: `[x, y, z, t]*`. Signature: `(`- - - +`)

 

_______________________________________________________

 

Physics:-g_[mu, nu] = Matrix(%id = 36893488152536018204)

 

_______________________________________________________

 

`Setting `*lowercaselatin_is*` letters to represent `*space*` indices`

 

`Defined objects with tensor properties`

 

{Physics:-D_[mu], Physics:-Dgamma[mu], Physics:-Psigma[mu], Physics:-Ricci[mu, nu], Physics:-Riemann[mu, nu, alpha, beta], Physics:-Weyl[mu, nu, alpha, beta], Physics:-d_[mu], eta[mu, nu], Physics:-g_[mu, nu], Physics:-gamma_[i, j], h[mu, nu], Physics:-Christoffel[mu, nu, alpha], Physics:-Einstein[mu, nu], Physics:-LeviCivita[alpha, beta, mu, nu], Physics:-SpaceTimeVector[mu](X)}

(10)

Your first question was about deriving this linearized form of the Ricci tensor shown in Wikipedia, entered above as (6). The following steps apply to any tensor that can be expressed in terms of the perturbation `h__μ,ν`. For that, you need to express everything in terms of the metric `g__μ,ν` and its derivatives, then Substitute "=g[mu,nu]=eta[mu,nu]+h[mu,nu]"., as you attempted, Josechecukara.  

 

Start rewriting the Ricci tensor in terms of the metric. You could do this in two ways:

Ricci[definition]

Physics:-Ricci[mu, nu] = Physics:-d_[alpha](Physics:-Christoffel[`~alpha`, mu, nu], [X])-Physics:-d_[nu](Physics:-Christoffel[`~alpha`, mu, alpha], [X])+Physics:-Christoffel[`~beta`, mu, nu]*Physics:-Christoffel[`~alpha`, beta, alpha]-Physics:-Christoffel[`~beta`, mu, alpha]*Physics:-Christoffel[`~alpha`, nu, beta]

(11)

convert(Physics[Ricci][mu, nu] = Physics[d_][alpha](Physics[Christoffel][`~alpha`, mu, nu], [X])-Physics[d_][nu](Physics[Christoffel][`~alpha`, mu, alpha], [X])+Physics[Christoffel][`~beta`, mu, nu]*Physics[Christoffel][`~alpha`, beta, alpha]-Physics[Christoffel][`~beta`, mu, alpha]*Physics[Christoffel][`~alpha`, nu, beta], g_, only = Christoffel)

Physics:-Ricci[mu, nu] = (1/2)*Physics:-d_[alpha](Physics:-g_[`~alpha`, `~beta`], [X])*(Physics:-d_[nu](Physics:-g_[beta, mu], [X])+Physics:-d_[mu](Physics:-g_[beta, nu], [X])-Physics:-d_[beta](Physics:-g_[mu, nu], [X]))+(1/2)*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[alpha](Physics:-d_[nu](Physics:-g_[beta, mu], [X]), [X])+Physics:-d_[alpha](Physics:-d_[mu](Physics:-g_[beta, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[beta](Physics:-g_[mu, nu], [X]), [X]))-(1/2)*Physics:-d_[nu](Physics:-g_[`~alpha`, `~beta`], [X])*Physics:-d_[mu](Physics:-g_[alpha, beta], [X])-(1/2)*Physics:-g_[`~alpha`, `~beta`]*Physics:-d_[mu](Physics:-d_[nu](Physics:-g_[alpha, beta], [X]), [X])+(1/4)*Physics:-g_[`~beta`, `~lambda`]*(Physics:-d_[nu](Physics:-g_[lambda, mu], [X])+Physics:-d_[mu](Physics:-g_[lambda, nu], [X])-Physics:-d_[lambda](Physics:-g_[mu, nu], [X]))*Physics:-g_[`~alpha`, `~kappa`]*Physics:-d_[beta](Physics:-g_[alpha, kappa], [X])-(1/4)*Physics:-g_[`~beta`, `~lambda`]*(Physics:-d_[mu](Physics:-g_[alpha, lambda], [X])+Physics:-d_[alpha](Physics:-g_[lambda, mu], [X])-Physics:-d_[lambda](Physics:-g_[alpha, mu], [X]))*Physics:-g_[`~alpha`, `~kappa`]*(Physics:-d_[nu](Physics:-g_[beta, kappa], [X])+Physics:-d_[beta](Physics:-g_[kappa, nu], [X])-Physics:-d_[kappa](Physics:-g_[beta, nu], [X]))

(12)

Or directly

Ricci[mu, nu] = convert(Ricci[mu, nu], g_)

Physics:-Ricci[mu, nu] = (1/2)*Physics:-d_[alpha](Physics:-g_[`~alpha`, `~beta`], [X])*(Physics:-d_[nu](Physics:-g_[beta, mu], [X])+Physics:-d_[mu](Physics:-g_[beta, nu], [X])-Physics:-d_[beta](Physics:-g_[mu, nu], [X]))+(1/2)*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[alpha](Physics:-d_[nu](Physics:-g_[beta, mu], [X]), [X])+Physics:-d_[alpha](Physics:-d_[mu](Physics:-g_[beta, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[beta](Physics:-g_[mu, nu], [X]), [X]))-(1/2)*Physics:-d_[nu](Physics:-g_[`~alpha`, `~beta`], [X])*Physics:-d_[mu](Physics:-g_[alpha, beta], [X])-(1/2)*Physics:-g_[`~alpha`, `~beta`]*Physics:-d_[mu](Physics:-d_[nu](Physics:-g_[alpha, beta], [X]), [X])+(1/4)*Physics:-g_[`~beta`, `~lambda`]*(Physics:-d_[nu](Physics:-g_[lambda, mu], [X])+Physics:-d_[mu](Physics:-g_[lambda, nu], [X])-Physics:-d_[lambda](Physics:-g_[mu, nu], [X]))*Physics:-g_[`~alpha`, `~kappa`]*Physics:-d_[beta](Physics:-g_[alpha, kappa], [X])-(1/4)*Physics:-g_[`~beta`, `~lambda`]*(Physics:-d_[mu](Physics:-g_[alpha, lambda], [X])+Physics:-d_[alpha](Physics:-g_[lambda, mu], [X])-Physics:-d_[lambda](Physics:-g_[alpha, mu], [X]))*Physics:-g_[`~alpha`, `~kappa`]*(Physics:-d_[nu](Physics:-g_[beta, kappa], [X])+Physics:-d_[beta](Physics:-g_[kappa, nu], [X])-Physics:-d_[kappa](Physics:-g_[beta, nu], [X]))

(13)

Introduce the perturbation using the metric's definition (or directly )

g_[definition]

Physics:-g_[mu, nu] = eta[mu, nu]+h[mu, nu]

(14)

Substitute(g_[definition], Physics[Ricci][mu, nu] = (1/2)*Physics[d_][alpha](Physics[g_][`~alpha`, `~beta`], [X])*(Physics[d_][nu](Physics[g_][beta, mu], [X])+Physics[d_][mu](Physics[g_][beta, nu], [X])-Physics[d_][beta](Physics[g_][mu, nu], [X]))+(1/2)*Physics[g_][`~alpha`, `~beta`]*(Physics[d_][alpha](Physics[d_][nu](Physics[g_][beta, mu], [X]), [X])+Physics[d_][alpha](Physics[d_][mu](Physics[g_][beta, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][beta](Physics[g_][mu, nu], [X]), [X]))-(1/2)*Physics[d_][nu](Physics[g_][`~alpha`, `~beta`], [X])*Physics[d_][mu](Physics[g_][alpha, beta], [X])-(1/2)*Physics[g_][`~alpha`, `~beta`]*Physics[d_][mu](Physics[d_][nu](Physics[g_][alpha, beta], [X]), [X])+(1/4)*Physics[g_][`~beta`, `~lambda`]*(Physics[d_][nu](Physics[g_][lambda, mu], [X])+Physics[d_][mu](Physics[g_][lambda, nu], [X])-Physics[d_][lambda](Physics[g_][mu, nu], [X]))*Physics[g_][`~alpha`, `~kappa`]*Physics[d_][beta](Physics[g_][alpha, kappa], [X])-(1/4)*Physics[g_][`~beta`, `~lambda`]*(Physics[d_][mu](Physics[g_][alpha, lambda], [X])+Physics[d_][alpha](Physics[g_][lambda, mu], [X])-Physics[d_][lambda](Physics[g_][alpha, mu], [X]))*Physics[g_][`~alpha`, `~kappa`]*(Physics[d_][nu](Physics[g_][beta, kappa], [X])+Physics[d_][beta](Physics[g_][kappa, nu], [X])-Physics[d_][kappa](Physics[g_][beta, nu], [X])))

Physics:-Ricci[mu, nu] = (1/2)*Physics:-d_[alpha](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*(Physics:-d_[nu](eta[beta, mu]+h[beta, mu], [X])+Physics:-d_[mu](eta[beta, nu]+h[beta, nu], [X])-Physics:-d_[beta](eta[mu, nu]+h[mu, nu], [X]))+(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics:-d_[alpha](Physics:-d_[nu](eta[beta, mu]+h[beta, mu], [X]), [X])+Physics:-d_[alpha](Physics:-d_[mu](eta[beta, nu]+h[beta, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[beta](eta[mu, nu]+h[mu, nu], [X]), [X]))-(1/2)*Physics:-d_[nu](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*Physics:-d_[mu](eta[alpha, beta]+h[alpha, beta], [X])-(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics:-d_[mu](Physics:-d_[nu](eta[alpha, beta]+h[alpha, beta], [X]), [X])+(1/4)*(eta[`~beta`, `~lambda`]+h[`~beta`, `~lambda`])*(Physics:-d_[nu](eta[lambda, mu]+h[lambda, mu], [X])+Physics:-d_[mu](eta[lambda, nu]+h[lambda, nu], [X])-Physics:-d_[lambda](eta[mu, nu]+h[mu, nu], [X]))*(eta[`~alpha`, `~kappa`]+h[`~alpha`, `~kappa`])*Physics:-d_[beta](eta[alpha, kappa]+h[alpha, kappa], [X])-(1/4)*(eta[`~beta`, `~lambda`]+h[`~beta`, `~lambda`])*(Physics:-d_[mu](eta[alpha, lambda]+h[alpha, lambda], [X])+Physics:-d_[alpha](eta[lambda, mu]+h[lambda, mu], [X])-Physics:-d_[lambda](eta[alpha, mu]+h[alpha, mu], [X]))*(eta[`~alpha`, `~kappa`]+h[`~alpha`, `~kappa`])*(Physics:-d_[nu](eta[beta, kappa]+h[beta, kappa], [X])+Physics:-d_[beta](eta[kappa, nu]+h[kappa, nu], [X])-Physics:-d_[kappa](eta[beta, nu]+h[beta, nu], [X]))

(15)

Linearize with respect to h[mu, nu] and its derivatives (note: this Library:-Linearize command is present only in the latest Maplesoft Physics Updates; to install before using this worksheet, input
Physics:-Version(latest)

Library:-Linearize(Physics[Ricci][mu, nu] = (1/2)*Physics[d_][alpha](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*(Physics[d_][nu](eta[beta, mu]+h[beta, mu], [X])+Physics[d_][mu](eta[beta, nu]+h[beta, nu], [X])-Physics[d_][beta](eta[mu, nu]+h[mu, nu], [X]))+(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics[d_][alpha](Physics[d_][nu](eta[beta, mu]+h[beta, mu], [X]), [X])+Physics[d_][alpha](Physics[d_][mu](eta[beta, nu]+h[beta, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][beta](eta[mu, nu]+h[mu, nu], [X]), [X]))-(1/2)*Physics[d_][nu](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*Physics[d_][mu](eta[alpha, beta]+h[alpha, beta], [X])-(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics[d_][mu](Physics[d_][nu](eta[alpha, beta]+h[alpha, beta], [X]), [X])+(1/4)*(eta[`~beta`, `~lambda`]+h[`~beta`, `~lambda`])*(Physics[d_][nu](eta[lambda, mu]+h[lambda, mu], [X])+Physics[d_][mu](eta[lambda, nu]+h[lambda, nu], [X])-Physics[d_][lambda](eta[mu, nu]+h[mu, nu], [X]))*(eta[`~alpha`, `~kappa`]+h[`~alpha`, `~kappa`])*Physics[d_][beta](eta[alpha, kappa]+h[alpha, kappa], [X])-(1/4)*(eta[`~beta`, `~lambda`]+h[`~beta`, `~lambda`])*(Physics[d_][mu](eta[alpha, lambda]+h[alpha, lambda], [X])+Physics[d_][alpha](eta[lambda, mu]+h[lambda, mu], [X])-Physics[d_][lambda](eta[alpha, mu]+h[alpha, mu], [X]))*(eta[`~alpha`, `~kappa`]+h[`~alpha`, `~kappa`])*(Physics[d_][nu](eta[beta, kappa]+h[beta, kappa], [X])+Physics[d_][beta](eta[kappa, nu]+h[kappa, nu], [X])-Physics[d_][kappa](eta[beta, nu]+h[beta, nu], [X])), h)

Physics:-Ricci[mu, nu] = -(1/2)*eta[`~alpha`, `~beta`]*(Physics:-d_[mu](Physics:-d_[nu](h[alpha, beta], [X]), [X])+Physics:-d_[alpha](Physics:-d_[beta](h[mu, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[nu](h[beta, mu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[mu](h[beta, nu], [X]), [X]))

(16)

At this point, we want to replace "eta[]^(alpha,beta)" by the Minkowski metric. The most direct way is to restate the metric as a flat Minkowski, but then Ricci automatically evaluates to 0, so first introduce its inert version %Ricci,

subs(Ricci = %Ricci, Physics[Ricci][mu, nu] = -(1/2)*eta[`~alpha`, `~beta`]*(Physics[d_][mu](Physics[d_][nu](h[alpha, beta], [X]), [X])+Physics[d_][alpha](Physics[d_][beta](h[mu, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][nu](h[beta, mu], [X]), [X])-Physics[d_][alpha](Physics[d_][mu](h[beta, nu], [X]), [X])))

%Ricci[mu, nu] = -(1/2)*eta[`~alpha`, `~beta`]*(Physics:-d_[mu](Physics:-d_[nu](h[alpha, beta], [X]), [X])+Physics:-d_[alpha](Physics:-d_[beta](h[mu, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[nu](h[beta, mu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[mu](h[beta, nu], [X]), [X]))

(17)

Set the metric to be flat Minkowski

g_[min]

_______________________________________________________

 

`The Minkowski metric in coordinates `*[x, y, z, t]

 

`Signature: `(`- - - +`)

 

_______________________________________________________

 

Physics:-g_[mu, nu] = Matrix(%id = 36893488152329528252)

(18)

Replace eta[mu, nu] by the metric g[mu, nu], then Simplify

subs(eta = g_, %Ricci[mu, nu] = -(1/2)*eta[`~alpha`, `~beta`]*(Physics[d_][mu](Physics[d_][nu](h[alpha, beta], [X]), [X])+Physics[d_][alpha](Physics[d_][beta](h[mu, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][nu](h[beta, mu], [X]), [X])-Physics[d_][alpha](Physics[d_][mu](h[beta, nu], [X]), [X])))

%Ricci[mu, nu] = -(1/2)*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[mu](Physics:-d_[nu](h[alpha, beta], [X]), [X])+Physics:-d_[alpha](Physics:-d_[beta](h[mu, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[nu](h[beta, mu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[mu](h[beta, nu], [X]), [X]))

(19)

Simplify(%Ricci[mu, nu] = -(1/2)*Physics[g_][`~alpha`, `~beta`]*(Physics[d_][mu](Physics[d_][nu](h[alpha, beta], [X]), [X])+Physics[d_][alpha](Physics[d_][beta](h[mu, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][nu](h[beta, mu], [X]), [X])-Physics[d_][alpha](Physics[d_][mu](h[beta, nu], [X]), [X])))

%Ricci[mu, nu] = (1/2)*Physics:-d_[beta](Physics:-d_[nu](h[mu, `~beta`], [X]), [X])-(1/2)*Physics:-dAlembertian(h[mu, nu], [X])-(1/2)*Physics:-d_[mu](Physics:-d_[nu](h[beta, `~beta`], [X]), [X])+(1/2)*Physics:-d_[beta](Physics:-d_[mu](h[nu, `~beta`], [X]), [X])

(20)

This expression is the same one shown in the Wikipedia, entered lines above as (6)

%Ricci[mu, nu] = (1/2)*Physics[d_][mu](Physics[d_][sigma](h[`~sigma`, nu], [X]), [X])+(1/2)*Physics[d_][nu](Physics[d_][sigma](h[`~sigma`, mu], [X]), [X])-(1/2)*Physics[d_][mu](Physics[d_][nu](h[`~sigma`, sigma], [X]), [X])-(1/2)*Physics[dAlembertian](h[mu, nu], [X])

%Ricci[mu, nu] = (1/2)*Physics:-d_[mu](Physics:-d_[sigma](h[`~sigma`, nu], [X]), [X])+(1/2)*Physics:-d_[nu](Physics:-d_[sigma](h[`~sigma`, mu], [X]), [X])-(1/2)*Physics:-d_[mu](Physics:-d_[nu](h[`~sigma`, sigma], [X]), [X])-(1/2)*Physics:-dAlembertian(h[mu, nu], [X])

(21)

Simplify((%Ricci[mu, nu] = (1/2)*Physics[d_][beta](Physics[d_][nu](h[mu, `~beta`], [X]), [X])-(1/2)*Physics[dAlembertian](h[mu, nu], [X])-(1/2)*Physics[d_][mu](Physics[d_][nu](h[beta, `~beta`], [X]), [X])+(1/2)*Physics[d_][beta](Physics[d_][mu](h[nu, `~beta`], [X]), [X]))-(%Ricci[mu, nu] = (1/2)*Physics[d_][mu](Physics[d_][sigma](h[`~sigma`, nu], [X]), [X])+(1/2)*Physics[d_][nu](Physics[d_][sigma](h[`~sigma`, mu], [X]), [X])-(1/2)*Physics[d_][mu](Physics[d_][nu](h[`~sigma`, sigma], [X]), [X])-(1/2)*Physics[dAlembertian](h[mu, nu], [X])))

0 = 0

(22)

The question by Josecherukara  is the same, just about the linearized form of Einstein's tensor, and the steps are the same. Again to avoid the automatic evaluation to 0 of  `Einstein__μ,ν`when in a flat space, redefine again the metric as in (9)

Physics[g_][mu, nu] = eta[mu, nu]+h[mu, nu]

Physics:-g_[mu, nu] = eta[mu, nu]+h[mu, nu]

(23)

Define(Physics[g_][mu, nu] = eta[mu, nu]+h[mu, nu])

_______________________________________________________

 

`Coordinates: `[x, y, z, t]*`. Signature: `(`- - - +`)

 

_______________________________________________________

 

Physics:-g_[mu, nu] = Matrix(%id = 36893488152323722524)

 

_______________________________________________________

 

`Defined objects with tensor properties`

 

{Physics:-D_[mu], Physics:-Dgamma[mu], Physics:-Psigma[mu], Physics:-Ricci[mu, nu], Physics:-Riemann[mu, nu, alpha, beta], Physics:-Weyl[mu, nu, alpha, beta], Physics:-d_[mu], eta[mu, nu], Physics:-g_[mu, nu], Physics:-gamma_[i, j], h[mu, nu], Physics:-Christoffel[mu, nu, alpha], Physics:-Einstein[mu, nu], Physics:-LeviCivita[alpha, beta, mu, nu], Physics:-SpaceTimeVector[mu](X)}

(24)

%Einstein[mu, nu] = convert(Einstein[mu, nu], g_)

%Einstein[mu, nu] = (1/2)*Physics:-d_[beta](Physics:-g_[`~alpha`, `~beta`], [X])*(Physics:-d_[nu](Physics:-g_[alpha, mu], [X])+Physics:-d_[mu](Physics:-g_[alpha, nu], [X])-Physics:-d_[alpha](Physics:-g_[mu, nu], [X]))+(1/2)*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[beta](Physics:-d_[nu](Physics:-g_[alpha, mu], [X]), [X])+Physics:-d_[beta](Physics:-d_[mu](Physics:-g_[alpha, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[beta](Physics:-g_[mu, nu], [X]), [X]))-(1/2)*Physics:-d_[nu](Physics:-g_[`~alpha`, `~beta`], [X])*Physics:-d_[mu](Physics:-g_[alpha, beta], [X])-(1/2)*Physics:-g_[`~alpha`, `~beta`]*Physics:-d_[mu](Physics:-d_[nu](Physics:-g_[alpha, beta], [X]), [X])+(1/4)*Physics:-g_[`~kappa`, `~lambda`]*(Physics:-d_[nu](Physics:-g_[lambda, mu], [X])+Physics:-d_[mu](Physics:-g_[lambda, nu], [X])-Physics:-d_[lambda](Physics:-g_[mu, nu], [X]))*Physics:-g_[`~alpha`, `~beta`]*Physics:-d_[kappa](Physics:-g_[alpha, beta], [X])-(1/4)*Physics:-g_[`~kappa`, `~lambda`]*(Physics:-d_[mu](Physics:-g_[beta, lambda], [X])+Physics:-d_[beta](Physics:-g_[lambda, mu], [X])-Physics:-d_[lambda](Physics:-g_[beta, mu], [X]))*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[nu](Physics:-g_[alpha, kappa], [X])+Physics:-d_[kappa](Physics:-g_[alpha, nu], [X])-Physics:-d_[alpha](Physics:-g_[kappa, nu], [X]))-(1/2)*Physics:-g_[mu, nu]*Physics:-g_[`~alpha`, `~tau`]*((1/2)*Physics:-d_[kappa](Physics:-g_[`~kappa`, `~sigma`], [X])*(Physics:-d_[tau](Physics:-g_[alpha, sigma], [X])+Physics:-d_[alpha](Physics:-g_[sigma, tau], [X])-Physics:-d_[sigma](Physics:-g_[alpha, tau], [X]))+(1/2)*Physics:-g_[`~kappa`, `~sigma`]*(Physics:-d_[kappa](Physics:-d_[tau](Physics:-g_[alpha, sigma], [X]), [X])+Physics:-d_[alpha](Physics:-d_[kappa](Physics:-g_[sigma, tau], [X]), [X])-Physics:-d_[kappa](Physics:-d_[sigma](Physics:-g_[alpha, tau], [X]), [X]))-(1/2)*Physics:-d_[tau](Physics:-g_[`~beta`, `~kappa`], [X])*Physics:-d_[alpha](Physics:-g_[beta, kappa], [X])-(1/2)*Physics:-g_[`~beta`, `~kappa`]*Physics:-d_[alpha](Physics:-d_[tau](Physics:-g_[beta, kappa], [X]), [X])+(1/4)*Physics:-g_[`~lambda`, `~omega`]*(Physics:-d_[tau](Physics:-g_[alpha, omega], [X])+Physics:-d_[alpha](Physics:-g_[omega, tau], [X])-Physics:-d_[omega](Physics:-g_[alpha, tau], [X]))*Physics:-g_[`~kappa`, `~upsilon`]*Physics:-d_[lambda](Physics:-g_[kappa, upsilon], [X])-(1/4)*Physics:-g_[`~lambda`, `~psi`]*(Physics:-d_[kappa](Physics:-g_[alpha, psi], [X])+Physics:-d_[alpha](Physics:-g_[kappa, psi], [X])-Physics:-d_[psi](Physics:-g_[alpha, kappa], [X]))*Physics:-g_[`~chi`, `~kappa`]*(Physics:-d_[tau](Physics:-g_[chi, lambda], [X])+Physics:-d_[lambda](Physics:-g_[chi, tau], [X])-Physics:-d_[chi](Physics:-g_[lambda, tau], [X])))

(25)

Substitute(Physics[g_][mu, nu] = eta[mu, nu]+h[mu, nu], %Einstein[mu, nu] = (1/2)*Physics[d_][beta](Physics[g_][`~alpha`, `~beta`], [X])*(Physics[d_][nu](Physics[g_][alpha, mu], [X])+Physics[d_][mu](Physics[g_][alpha, nu], [X])-Physics[d_][alpha](Physics[g_][mu, nu], [X]))+(1/2)*Physics[g_][`~alpha`, `~beta`]*(Physics[d_][beta](Physics[d_][nu](Physics[g_][alpha, mu], [X]), [X])+Physics[d_][beta](Physics[d_][mu](Physics[g_][alpha, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][beta](Physics[g_][mu, nu], [X]), [X]))-(1/2)*Physics[d_][nu](Physics[g_][`~alpha`, `~beta`], [X])*Physics[d_][mu](Physics[g_][alpha, beta], [X])-(1/2)*Physics[g_][`~alpha`, `~beta`]*Physics[d_][mu](Physics[d_][nu](Physics[g_][alpha, beta], [X]), [X])+(1/4)*Physics[g_][`~kappa`, `~lambda`]*(Physics[d_][nu](Physics[g_][lambda, mu], [X])+Physics[d_][mu](Physics[g_][lambda, nu], [X])-Physics[d_][lambda](Physics[g_][mu, nu], [X]))*Physics[g_][`~alpha`, `~beta`]*Physics[d_][kappa](Physics[g_][alpha, beta], [X])-(1/4)*Physics[g_][`~kappa`, `~lambda`]*(Physics[d_][mu](Physics[g_][beta, lambda], [X])+Physics[d_][beta](Physics[g_][lambda, mu], [X])-Physics[d_][lambda](Physics[g_][beta, mu], [X]))*Physics[g_][`~alpha`, `~beta`]*(Physics[d_][nu](Physics[g_][alpha, kappa], [X])+Physics[d_][kappa](Physics[g_][alpha, nu], [X])-Physics[d_][alpha](Physics[g_][kappa, nu], [X]))-(1/2)*Physics[g_][mu, nu]*Physics[g_][`~alpha`, `~tau`]*((1/2)*Physics[d_][kappa](Physics[g_][`~kappa`, `~sigma`], [X])*(Physics[d_][tau](Physics[g_][alpha, sigma], [X])+Physics[d_][alpha](Physics[g_][sigma, tau], [X])-Physics[d_][sigma](Physics[g_][alpha, tau], [X]))+(1/2)*Physics[g_][`~kappa`, `~sigma`]*(Physics[d_][kappa](Physics[d_][tau](Physics[g_][alpha, sigma], [X]), [X])+Physics[d_][alpha](Physics[d_][kappa](Physics[g_][sigma, tau], [X]), [X])-Physics[d_][kappa](Physics[d_][sigma](Physics[g_][alpha, tau], [X]), [X]))-(1/2)*Physics[d_][tau](Physics[g_][`~beta`, `~kappa`], [X])*Physics[d_][alpha](Physics[g_][beta, kappa], [X])-(1/2)*Physics[g_][`~beta`, `~kappa`]*Physics[d_][alpha](Physics[d_][tau](Physics[g_][beta, kappa], [X]), [X])+(1/4)*Physics[g_][`~lambda`, `~omega`]*(Physics[d_][tau](Physics[g_][alpha, omega], [X])+Physics[d_][alpha](Physics[g_][omega, tau], [X])-Physics[d_][omega](Physics[g_][alpha, tau], [X]))*Physics[g_][`~kappa`, `~upsilon`]*Physics[d_][lambda](Physics[g_][kappa, upsilon], [X])-(1/4)*Physics[g_][`~lambda`, `~psi`]*(Physics[d_][kappa](Physics[g_][alpha, psi], [X])+Physics[d_][alpha](Physics[g_][kappa, psi], [X])-Physics[d_][psi](Physics[g_][alpha, kappa], [X]))*Physics[g_][`~chi`, `~kappa`]*(Physics[d_][tau](Physics[g_][chi, lambda], [X])+Physics[d_][lambda](Physics[g_][chi, tau], [X])-Physics[d_][chi](Physics[g_][lambda, tau], [X]))))

%Einstein[mu, nu] = (1/2)*Physics:-d_[beta](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*(Physics:-d_[nu](eta[alpha, mu]+h[alpha, mu], [X])+Physics:-d_[mu](eta[alpha, nu]+h[alpha, nu], [X])-Physics:-d_[alpha](eta[mu, nu]+h[mu, nu], [X]))+(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics:-d_[beta](Physics:-d_[nu](eta[alpha, mu]+h[alpha, mu], [X]), [X])+Physics:-d_[beta](Physics:-d_[mu](eta[alpha, nu]+h[alpha, nu], [X]), [X])-Physics:-d_[alpha](Physics:-d_[beta](eta[mu, nu]+h[mu, nu], [X]), [X]))-(1/2)*Physics:-d_[nu](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*Physics:-d_[mu](eta[alpha, beta]+h[alpha, beta], [X])-(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics:-d_[mu](Physics:-d_[nu](eta[alpha, beta]+h[alpha, beta], [X]), [X])+(1/4)*(eta[`~kappa`, `~lambda`]+h[`~kappa`, `~lambda`])*(Physics:-d_[nu](eta[lambda, mu]+h[lambda, mu], [X])+Physics:-d_[mu](eta[lambda, nu]+h[lambda, nu], [X])-Physics:-d_[lambda](eta[mu, nu]+h[mu, nu], [X]))*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics:-d_[kappa](eta[alpha, beta]+h[alpha, beta], [X])-(1/4)*(eta[`~kappa`, `~lambda`]+h[`~kappa`, `~lambda`])*(Physics:-d_[mu](eta[beta, lambda]+h[beta, lambda], [X])+Physics:-d_[beta](eta[lambda, mu]+h[lambda, mu], [X])-Physics:-d_[lambda](eta[beta, mu]+h[beta, mu], [X]))*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics:-d_[nu](eta[alpha, kappa]+h[alpha, kappa], [X])+Physics:-d_[kappa](eta[alpha, nu]+h[alpha, nu], [X])-Physics:-d_[alpha](eta[kappa, nu]+h[kappa, nu], [X]))+(1/8)*((-2*h[`~kappa`, `~sigma`]-2*eta[`~kappa`, `~sigma`])*Physics:-d_[alpha](Physics:-d_[kappa](eta[sigma, tau]+h[sigma, tau], [X]), [X])+(2*h[`~beta`, `~kappa`]+2*eta[`~beta`, `~kappa`])*Physics:-d_[alpha](Physics:-d_[tau](eta[beta, kappa]+h[beta, kappa], [X]), [X])+(2*h[`~kappa`, `~sigma`]+2*eta[`~kappa`, `~sigma`])*Physics:-d_[kappa](Physics:-d_[sigma](eta[alpha, tau]+h[alpha, tau], [X]), [X])+(-2*h[`~kappa`, `~sigma`]-2*eta[`~kappa`, `~sigma`])*Physics:-d_[kappa](Physics:-d_[tau](eta[alpha, sigma]+h[alpha, sigma], [X]), [X])+(Physics:-d_[tau](eta[chi, lambda]+h[chi, lambda], [X])+Physics:-d_[lambda](eta[chi, tau]+h[chi, tau], [X])-Physics:-d_[chi](eta[lambda, tau]+h[lambda, tau], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics:-d_[alpha](eta[kappa, psi]+h[kappa, psi], [X])-Physics:-d_[lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics:-d_[alpha](eta[omega, tau]+h[omega, tau], [X])-2*Physics:-d_[kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics:-d_[alpha](eta[sigma, tau]+h[sigma, tau], [X])-(Physics:-d_[kappa](eta[alpha, psi]+h[alpha, psi], [X])-Physics:-d_[psi](eta[alpha, kappa]+h[alpha, kappa], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics:-d_[chi](eta[lambda, tau]+h[lambda, tau], [X])+(Physics:-d_[lambda](eta[chi, tau]+h[chi, tau], [X])+Physics:-d_[tau](eta[chi, lambda]+h[chi, lambda], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics:-d_[kappa](eta[alpha, psi]+h[alpha, psi], [X])+Physics:-d_[lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics:-d_[omega](eta[alpha, tau]+h[alpha, tau], [X])+2*Physics:-d_[kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics:-d_[sigma](eta[alpha, tau]+h[alpha, tau], [X])-Physics:-d_[psi](eta[alpha, kappa]+h[alpha, kappa], [X])*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics:-d_[tau](eta[chi, lambda]+h[chi, lambda], [X])+2*Physics:-d_[tau](eta[`~beta`, `~kappa`]+h[`~beta`, `~kappa`], [X])*Physics:-d_[alpha](eta[beta, kappa]+h[beta, kappa], [X])-Physics:-d_[lambda](eta[chi, tau]+h[chi, tau], [X])*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics:-d_[psi](eta[alpha, kappa]+h[alpha, kappa], [X])-Physics:-d_[lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics:-d_[tau](eta[alpha, omega]+h[alpha, omega], [X])-2*Physics:-d_[kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics:-d_[tau](eta[alpha, sigma]+h[alpha, sigma], [X]))*(eta[mu, nu]+h[mu, nu])*(eta[`~alpha`, `~tau`]+h[`~alpha`, `~tau`])

(26)

Library:-Linearize(%Einstein[mu, nu] = (1/2)*Physics[d_][beta](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*(Physics[d_][nu](eta[alpha, mu]+h[alpha, mu], [X])+Physics[d_][mu](eta[alpha, nu]+h[alpha, nu], [X])-Physics[d_][alpha](eta[mu, nu]+h[mu, nu], [X]))+(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics[d_][beta](Physics[d_][nu](eta[alpha, mu]+h[alpha, mu], [X]), [X])+Physics[d_][beta](Physics[d_][mu](eta[alpha, nu]+h[alpha, nu], [X]), [X])-Physics[d_][alpha](Physics[d_][beta](eta[mu, nu]+h[mu, nu], [X]), [X]))-(1/2)*Physics[d_][nu](eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`], [X])*Physics[d_][mu](eta[alpha, beta]+h[alpha, beta], [X])-(1/2)*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics[d_][mu](Physics[d_][nu](eta[alpha, beta]+h[alpha, beta], [X]), [X])+(1/4)*(eta[`~kappa`, `~lambda`]+h[`~kappa`, `~lambda`])*(Physics[d_][nu](eta[lambda, mu]+h[lambda, mu], [X])+Physics[d_][mu](eta[lambda, nu]+h[lambda, nu], [X])-Physics[d_][lambda](eta[mu, nu]+h[mu, nu], [X]))*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*Physics[d_][kappa](eta[alpha, beta]+h[alpha, beta], [X])-(1/4)*(eta[`~kappa`, `~lambda`]+h[`~kappa`, `~lambda`])*(Physics[d_][mu](eta[beta, lambda]+h[beta, lambda], [X])+Physics[d_][beta](eta[lambda, mu]+h[lambda, mu], [X])-Physics[d_][lambda](eta[beta, mu]+h[beta, mu], [X]))*(eta[`~alpha`, `~beta`]+h[`~alpha`, `~beta`])*(Physics[d_][nu](eta[alpha, kappa]+h[alpha, kappa], [X])+Physics[d_][kappa](eta[alpha, nu]+h[alpha, nu], [X])-Physics[d_][alpha](eta[kappa, nu]+h[kappa, nu], [X]))+(1/8)*((-2*h[`~kappa`, `~sigma`]-2*eta[`~kappa`, `~sigma`])*Physics[d_][alpha](Physics[d_][kappa](eta[sigma, tau]+h[sigma, tau], [X]), [X])+(2*h[`~beta`, `~kappa`]+2*eta[`~beta`, `~kappa`])*Physics[d_][alpha](Physics[d_][tau](eta[beta, kappa]+h[beta, kappa], [X]), [X])+(2*h[`~kappa`, `~sigma`]+2*eta[`~kappa`, `~sigma`])*Physics[d_][kappa](Physics[d_][sigma](eta[alpha, tau]+h[alpha, tau], [X]), [X])+(-2*h[`~kappa`, `~sigma`]-2*eta[`~kappa`, `~sigma`])*Physics[d_][kappa](Physics[d_][tau](eta[alpha, sigma]+h[alpha, sigma], [X]), [X])+(Physics[d_][tau](eta[chi, lambda]+h[chi, lambda], [X])+Physics[d_][lambda](eta[chi, tau]+h[chi, tau], [X])-Physics[d_][chi](eta[lambda, tau]+h[lambda, tau], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics[d_][alpha](eta[kappa, psi]+h[kappa, psi], [X])-Physics[d_][lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics[d_][alpha](eta[omega, tau]+h[omega, tau], [X])-2*Physics[d_][kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics[d_][alpha](eta[sigma, tau]+h[sigma, tau], [X])-(Physics[d_][kappa](eta[alpha, psi]+h[alpha, psi], [X])-Physics[d_][psi](eta[alpha, kappa]+h[alpha, kappa], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics[d_][chi](eta[lambda, tau]+h[lambda, tau], [X])+(Physics[d_][lambda](eta[chi, tau]+h[chi, tau], [X])+Physics[d_][tau](eta[chi, lambda]+h[chi, lambda], [X]))*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics[d_][kappa](eta[alpha, psi]+h[alpha, psi], [X])+Physics[d_][lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics[d_][omega](eta[alpha, tau]+h[alpha, tau], [X])+2*Physics[d_][kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics[d_][sigma](eta[alpha, tau]+h[alpha, tau], [X])-Physics[d_][psi](eta[alpha, kappa]+h[alpha, kappa], [X])*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics[d_][tau](eta[chi, lambda]+h[chi, lambda], [X])+2*Physics[d_][tau](eta[`~beta`, `~kappa`]+h[`~beta`, `~kappa`], [X])*Physics[d_][alpha](eta[beta, kappa]+h[beta, kappa], [X])-Physics[d_][lambda](eta[chi, tau]+h[chi, tau], [X])*(eta[`~lambda`, `~psi`]+h[`~lambda`, `~psi`])*(eta[`~chi`, `~kappa`]+h[`~chi`, `~kappa`])*Physics[d_][psi](eta[alpha, kappa]+h[alpha, kappa], [X])-Physics[d_][lambda](eta[kappa, upsilon]+h[kappa, upsilon], [X])*(eta[`~kappa`, `~upsilon`]+h[`~kappa`, `~upsilon`])*(eta[`~lambda`, `~omega`]+h[`~lambda`, `~omega`])*Physics[d_][tau](eta[alpha, omega]+h[alpha, omega], [X])-2*Physics[d_][kappa](eta[`~kappa`, `~sigma`]+h[`~kappa`, `~sigma`], [X])*Physics[d_][tau](eta[alpha, sigma]+h[alpha, sigma], [X]))*(eta[mu, nu]+h[mu, nu])*(eta[`~alpha`, `~tau`]+h[`~alpha`, `~tau`]), h)

%Einstein[mu, nu] = -(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics:-d_[alpha](Physics:-d_[kappa](h[sigma, tau], [X]), [X])+(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics:-d_[kappa](Physics:-d_[sigma](h[alpha, tau], [X]), [X])-(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics:-d_[kappa](Physics:-d_[tau](h[alpha, sigma], [X]), [X])+(1/4)*eta[`~beta`, `~kappa`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics:-d_[alpha](Physics:-d_[tau](h[beta, kappa], [X]), [X])-(1/2)*eta[`~alpha`, `~beta`]*(Physics:-d_[alpha](Physics:-d_[beta](h[mu, nu], [X]), [X])-Physics:-d_[beta](Physics:-d_[mu](h[alpha, nu], [X]), [X])-Physics:-d_[beta](Physics:-d_[nu](h[alpha, mu], [X]), [X])+Physics:-d_[mu](Physics:-d_[nu](h[alpha, beta], [X]), [X]))

(27)

Reintroduce the Minkowski metric:

g_[min]

_______________________________________________________

 

`The Minkowski metric in coordinates `*[x, y, z, t]

 

`Signature: `(`- - - +`)

 

_______________________________________________________

 

Physics:-g_[mu, nu] = Matrix(%id = 36893488152325079036)

(28)

Replace occurrences of eta[mu, nu] by g[mu, nu] to proceed simplifying

subs(eta = g_, %Einstein[mu, nu] = -(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics[d_][alpha](Physics[d_][kappa](h[sigma, tau], [X]), [X])+(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics[d_][kappa](Physics[d_][sigma](h[alpha, tau], [X]), [X])-(1/4)*eta[`~kappa`, `~sigma`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics[d_][kappa](Physics[d_][tau](h[alpha, sigma], [X]), [X])+(1/4)*eta[`~beta`, `~kappa`]*eta[mu, nu]*eta[`~alpha`, `~tau`]*Physics[d_][alpha](Physics[d_][tau](h[beta, kappa], [X]), [X])-(1/2)*eta[`~alpha`, `~beta`]*(Physics[d_][alpha](Physics[d_][beta](h[mu, nu], [X]), [X])-Physics[d_][beta](Physics[d_][mu](h[alpha, nu], [X]), [X])-Physics[d_][beta](Physics[d_][nu](h[alpha, mu], [X]), [X])+Physics[d_][mu](Physics[d_][nu](h[alpha, beta], [X]), [X])))

%Einstein[mu, nu] = -(1/4)*Physics:-g_[`~kappa`, `~sigma`]*Physics:-g_[mu, nu]*Physics:-g_[`~alpha`, `~tau`]*Physics:-d_[alpha](Physics:-d_[kappa](h[sigma, tau], [X]), [X])+(1/4)*Physics:-g_[`~kappa`, `~sigma`]*Physics:-g_[mu, nu]*Physics:-g_[`~alpha`, `~tau`]*Physics:-d_[kappa](Physics:-d_[sigma](h[alpha, tau], [X]), [X])-(1/4)*Physics:-g_[`~kappa`, `~sigma`]*Physics:-g_[mu, nu]*Physics:-g_[`~alpha`, `~tau`]*Physics:-d_[kappa](Physics:-d_[tau](h[alpha, sigma], [X]), [X])+(1/4)*Physics:-g_[`~beta`, `~kappa`]*Physics:-g_[mu, nu]*Physics:-g_[`~alpha`, `~tau`]*Physics:-d_[alpha](Physics:-d_[tau](h[beta, kappa], [X]), [X])-(1/2)*Physics:-g_[`~alpha`, `~beta`]*(Physics:-d_[alpha](Physics:-d_[beta](h[mu, nu], [X]), [X])-Physics:-d_[beta](Physics:-d_[mu](h[alpha, nu], [X]), [X])-Physics:-d_[beta](Physics:-d_[nu](h[alpha, mu], [X]), [X])+Physics:-d_[mu](Physics:-d_[nu](h[alpha, beta], [X]), [X]))

(29)

And this is the linearized form of Einstein's tensor

Simplify(%Einstein[mu, nu] = -(1/4)*Physics[g_][`~kappa`, `~sigma`]*Physics[g_][mu, nu]*Physics[g_][`~alpha`, `~tau`]*Physics[d_][alpha](Physics[d_][kappa](h[sigma, tau], [X]), [X])+(1/4)*Physics[g_][`~kappa`, `~sigma`]*Physics[g_][mu, nu]*Physics[g_][`~alpha`, `~tau`]*Physics[d_][kappa](Physics[d_][sigma](h[alpha, tau], [X]), [X])-(1/4)*Physics[g_][`~kappa`, `~sigma`]*Physics[g_][mu, nu]*Physics[g_][`~alpha`, `~tau`]*Physics[d_][kappa](Physics[d_][tau](h[alpha, sigma], [X]), [X])+(1/4)*Physics[g_][`~beta`, `~kappa`]*Physics[g_][mu, nu]*Physics[g_][`~alpha`, `~tau`]*Physics[d_][alpha](Physics[d_][tau](h[beta, kappa], [X]), [X])-(1/2)*Physics[g_][`~alpha`, `~beta`]*(Physics[d_][alpha](Physics[d_][beta](h[mu, nu], [X]), [X])-Physics[d_][beta](Physics[d_][mu](h[alpha, nu], [X]), [X])-Physics[d_][beta](Physics[d_][nu](h[alpha, mu], [X]), [X])+Physics[d_][mu](Physics[d_][nu](h[alpha, beta], [X]), [X])))

%Einstein[mu, nu] = (1/2)*(-Physics:-d_[beta](Physics:-d_[sigma](h[`~beta`, `~sigma`], [X]), [X])+Physics:-dAlembertian(h[beta, `~beta`], [X]))*Physics:-g_[mu, nu]-(1/2)*Physics:-dAlembertian(h[mu, nu], [X])+(1/2)*Physics:-d_[beta](Physics:-d_[mu](h[nu, `~beta`], [X]), [X])+(1/2)*Physics:-d_[beta](Physics:-d_[nu](h[mu, `~beta`], [X]), [X])-(1/2)*Physics:-d_[mu](Physics:-d_[nu](h[beta, `~beta`], [X]), [X])

(30)

 


 

Download perturbed_metric_(general_case_reviewed).mw

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft​​​​​​​

@KM4 , @Josecherukara
I will return to this today, or tomorrow, and prepare a worksheet showing how.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

 

1 2 3 4 5 6 7 Last Page 1 of 62