lemelinm

1425 Reputation

14 Badges

14 years, 285 days

 

 

--------------------------------------
Mario Lemelin
Maple 14.00 Win 7 64 bits
Maple 14.00 Ubuntu 10,04 64 bits
messagerie : mario.lemelin@cgocable.ca téléphone :  (819) 376-0987

MaplePrimes Activity


These are replies submitted by lemelinm

@acer 

As you will see in the attached document, I don't get the true numbers on the vertical axis but the unit is good. The horizontal axis is OK.  Then equation (4) print OK but why (6) is not?  Do you see what I am doing wrong?  I am using Maple 2020.1

Thank you in advance for your help.

Mario

Quark-Quark_interaction.mw

@acer 

Hello to you.  You have well understood the question.  The use of Units:-UseUnit(MeV) is new to me and I will explore it.  It's a nice trick.  And why I didn't think of this one: combine(pmin, units).

Thank you for your fast response.  I think that it would be nice to have this option in the Physics package: with(Physics, Particles) with the possibility to call a particle (for example Pi+) with the values of I, P, Q, S, C, and B.  Mass in MeV/c^2 and lifetime could be interesting too.  Anyway, this idea came to me while I am reading for the pleasure of the following book: Griffiths, Introduction to Elementary Particles, 2008.

Thank you, Acer.

 

@longrob 

Just to mention that 10 years later, I come accross that simple way to dol it in Maple 2020.  So I was able to do manually the calculation and check with Maple if I was right.  Thank you.

Since I will do 3 courses in Quantum Mechanics (MITx), I am wondering if it would be an added bonus to use the physics package?

 

Thank you again for the fast responses.

@ecterrab 

With you answer, I found out that what I can do is:

> DifferentialGeometry:-Library:-MetricSearch()

And I can find all the metrics as we can see in the book: Stephani - Exact Solutions of Einstein's Field Equations -Cambridge University

That is great.  Thank you!

@Christopher2222 

If I put t := 3 [s], I want to be able to transform it in [m] by using the contextual menu..  But for Maple, this is not possible.  I have to do it manually.  When I load the Physics package, I would like that units take into account that c = 1, etc).

Is this possible or too complex?

@tomleslie 

When you decide that c=1, because the equations are simpler, this make it problematic for an introduction on general relativity.  So I would like to be able to write an equation like this -t^2+x^2+y^2+z^2 and have the result in meter.  That help to check if your calculation is coherent.  This c=1 has an effect everywhere in the theory and I was albe to see in different instances that it complexify the comprehension for the beginners.

@Christopher2222 

Thank you for the reference.

Mario

@John Fredsted 

Here the real commands:

>Setup(metric = Matrix(3, 3, [[1/(-r^2+1), 0, 0], [0, r^2, 0], [0, 0, r^2*sin(theta)^2]], shape = symmetric))

Isn't that fun!  Having done a course on General Relativity, it help to understand what we are doing.

Mario

@John Fredsted 

Thank you again John for your help on this.  Even if I do it in 4D, I only need to put g[0,0]=1.  But your idea help me do that:

restart;
with(Physics); Setup(mathematicalnotation = true, dimension = 3); Coordinates(X = (r, theta, phi)); Setup(metric = Matrix(3, 3, [[Physics:-`^`(1-Physics:-`^`(r, 2), -1), 0, 0], [0, Physics:-`^`(r, 2), 0], [0, 0, Physics:-`*`(Physics:-`^`(r, 2), Physics:-`^`(sin(theta), 2))]], shape = symmetric));

g_[];
Define(A, B);
Setup(differentiationvariables = X);
seq(Christoffel[j, mu, nu, matrix], j = [`~1`, `~2`, `~3`]);
seq(Ricci[mu, nu, matrix], mu = 1 .. 3);
Ricci[scalar];
Einstein[mu, nu, matrix];

I am starting to manage well the physics package.

Mario

 

It would be great to have a book on general relativity with commands, templates for the most used metric and exeamples.  It should take into account that it would need two parts, one for introductory courses and one for advanced users.  I tried to find books like that but with no success.

You deserve all those thanks that you receive.  For me in particular, I saw on you a mentor for working with Maple.  You were always there when I needed.  So I really hope that your mate will be happy and yourself should feel happy too for all the helps you give to so numerous people who wanted to start with Maple.  And yes, clickable calculus is a big success of yours

Thank you again Robert,

Mario Lemelin

@John Fredsted 

First of all, a big thank you for all the work you put on this post.

Secondly, when I wrote the Riemann tensor, I was in an angry mood.  At that point, I had decided to do it by hand, since it took too much effort to do something that it was supposed to be simple in the Physics package.  It is my impression that the package, for the moment, is built to get faster calculation of known metrics.  I think that a function like linearize should be part of the package.  That would help us to test other kind of metric that are not listed in the metric list.

So I will stop working in this because I have other calculations to make.  So thank you again for your help.  Maybe later will I come back to this.

@John Fredsted 

I am sorry that the new linearize does not work with seconde derivative.  But to do what I need, we can use the definition that we get from

linearize(Christoffel[~1,2,3]);

Now, the calculation I want to do is this one, meaning the rhs of it

Can this be done.

Thanks again for your help

@John Fredsted 

First, thank you for this very well explain steps.  The first thing I noticed is that g_[mu,nu] must be equal to eta_[mu,nu]+h[mu,nu].  This is important because eta will be apply on the first derivative.  In fact, I should show you what I want to calculate.

Secondly, I just learn that semi-colon is the covariant derivative while the comma mean the standard partial derivative.  It seems to be somewhat different in books that I read.

And the linearize procedure is very interresting.  So I will play a bit with it and if I ran into problems, I will contact you in this post.

Thanks again,

Mario

1 2 3 4 5 6 Page 1 of 6