MaplePrimes Questions

R dJ(t)/dt+J(t)/C=f(t)

where f(t) is a driving electromotive force. Use the fourier transform to analyze this equation as follows.

 

 

Find the transfer function G(alpha)  then find g(t) .

 Thanks ....

Hi ;

I need your help to write the system contains all these equation:



# system 1
u[0,0]:=(1/2)*(u[1,0]+u[0,1]);
u[0,N+1]:=(1/2)*(u[0,N]+u[1,N+1]);
u[N+1,N+1]:=(1/2)*(u[N+1,N]+u[N,N+1]);
u[N+1,0]:=(1/2)*(u[N,0]+u[N+1,1]);
# system 2
for j from 1 to N do
u[0, j] := (1/4)*(u[1,j]+u[1,j]+u[0, j-1]+u[0,j+1]-f[0,j]*h^2);
end do;
# System 3
for j from 1 to N do
u[N+1, j] := (1/4)*(u[N, j]+u[N-1,j]+u[N+1, j-1]+u[N+1, j+1]-f[N+1,j]*h^2);
end do;
system4
 eqs := [ seq(seq(Stencil[1](h,i,j,u,f),i=1..N),j=1..N)];

How can collect these system 1 system2, system3 ans system4 in a set with one name.


sys:=[eqs,system3,system2,system2]; ( sys: here contains all the equation).

Thanks you.

 

Hi,

When I run two times my code, the results change.

Have you any idea, why???

 

 

Hi:

i will calculation of value function in every step by 0.001(x=0,x=0.001,x=0.002,...x=10):

my function:

f(x)=x^3+2*x+1

my domain:

0<x<10

Hi , everyone who love Maple and dsolve command, 

my ODE is :

sys_ode := diff(d11(m), m) = -(3*sin(m)^2-1)*d31(m)/a^(3/2)+(-3*cos(m)*sin(m)/a^(3/2))*d41(m), diff(d21(m), m) = (-3*cos(m)*sin(m)/a^(3/2))*d31(m)-(3*cos(m)^2-1)*d41(m)/a^(3/2), diff(d31(m), m) = -a^(3/2)*d11(m), diff(d41(m), m) = -a^(3/2)*d21(m)

using " dsolve([sys_ode]) " command could get the solution easily, and the solution contains "I" (imaginary domain).

However, when we substitute the solution into the ODE "sys_ode", find not correct !

we use the following command to check the solution :

 simplify(  -diff(d11(m), m) -(3*sin(m)^2-1)*d31(m)/a^(3/2)+(-3*cos(m)*sin(m)/a^(3/2))*d41(m)  )

the upper expression is supposed to be zero, but not ! Is it a bug in Maple dsolve ?

Hi:

how can i convert one nonlinear second order ode to first order in maple?

 

How can I plot multiple numeric plots, such that numeric solutions belong to different differential equations?

For solving problem sets, I have a pdf template I created for myself that has a header with a blank for the class name, TA, professor, date, etc. In addition to this header, I had a margin on the left to scribble questions I had and to holepunch. 

 

I used to print out the template and write on the template and turn in that as my pset.

 

I am nowthinking of doing everything on the computer. Writing out all of the problem set on the computer. Combining stuff from maple, combining handwritten stuff from the computer using a digitizer. However I want to write it all on top of the my template that I created, which is a pdf file. I can turn the pdf into an image file if need be. 

 

What would be the easiest way to do what I want? To open a program that automatically sets that pdf as the template and easily lets me handwrite stuff I want and paste in maple code? 

 

Right now if I tried my idea, I would basically be constantly copying and pasting stuff from maple and my digitizer drawn pictures/equations into one file and it would be very clumsy.

 

Basically there are problems that I do partially on maple and I just want to unify all my work into one easy, printable file. 

I'm doing a pset in maple.

 

I paste in a block of text in maple. I put an asterik at front and end, and hit enter and it spits out an error. Why? I thought two *'s around it comment it out. 

 

I also can't see the file I uploaded while I preview my post. 

 

 

Hai everyone.

I try to double integrate this generalized extreme distribution.

q[p] := 6.256: h := .8; t[c] := .45: S[di] := 0: k[v] := .32639: mu[v] := -0.1786e-1: sigma[v] := 2.1694: k[t] := .36132: mu[t] := .63543: sigma[t] := 3.1183:


int(exp(-(1+k[v]*(v-mu[v])/sigma[v])^(-1/k[v]))*(1+k[v]*(v-mu[v])/sigma[v])^(-1-1/k[v])*exp(-(1+k[t]*(t-mu[t])/sigma[t])^(-1/k[t]))*(1+k[t]*(t-mu[t])/sigma[t])^(-1-1/k[t])/(sigma[v]*sigma[t]), [v = q[p]/(2*h*t)+q[p]*t[c]/(2*h)+S[di] .. infinity, t = 0 .. infinity]);

however, I got an error, as follow:

Error, (in assuming) when calling '`root/fraction`'. Received: 'numeric exception: overflow'

Any recomendations and tips to solve this integration? or this integration may cannot solve?

Thank you.

 

I have an expression like this:

Since it is linear I want Maple to rewrite it into this:

(with the benefit that Maple then can solve it at least up to a point). i have tried to conceive a rule to do that but got stuck relatively quickly. Does anybody have a way to do this (in some genrality)?

Thanks,

Mac Dude.

 

Hi, I want to define functions recursively... 

I don't know how to do it with the for loop in Maple.

I have a[0](t)=0, a[1](t)=t, and then some recursion connecting a[j+1](t) = f( a[j](t), a[j-1](t)) for some explicit function f. 

Then I want to plot the graph of a[N](t) as a function of t. 

Thanks!!

Hi, please a Need a help.
I have an error in my code: Error, (in collect) cannot collect 0
Here, is the code:
Question8X.mw

Many thinks.

 

Hello,

I am a student and using Maple to type homework assignments because of the math symbols available.  I was using Maple 17 but upgraded to 18, and its not as easy to use as 17.  For right now I am using it in text mode, because that is what I need.  But I can't figure out an easy way to do subscripts and superscripts in text mode.  Also, the "element symbol," where is it?  I feel like a lot of stuff is missing.  For instance, the arrows, the sideways triangle, lots of symbols I used before the upgrade I cant find.  Can someone help please???

After using Simplify the indices are are arranged in the tensor.  I am using the April 14th update from the Physics R&D page.

 


restart

with(Physics):

Setup(mathematicalnotation = true, coordinatesystems = X)

[coordinatesystems = {X}, mathematicalnotation = true]

(1)

Define(l[mu], eta[mu, nu] = -rhs(g_[Minkowski]))

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

(2)

declare(l(X))

l(x1, x2, x3, x4)*`will now be displayed as`*l

(3)

InitialMetric := g_[mu, nu] = eta[mu, nu]+Physics:-`*`(l[mu](X), l[nu](X)); 1; Define(G[mu, nu] = rhs(InitialMetric))

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

(4)

Setup(metric = rhs(G[]))

[metric = {(1, 1) = 1+l[1](X)^2, (1, 2) = l[1](X)*l[2](X), (1, 3) = l[1](X)*l[3](X), (1, 4) = l[1](X)*l[4](X), (2, 2) = 1+l[2](X)^2, (2, 3) = l[2](X)*l[3](X), (2, 4) = l[2](X)*l[4](X), (3, 3) = 1+l[3](X)^2, (3, 4) = l[3](X)*l[4](X), (4, 4) = -1+l[4](X)^2}]

(5)

NULL

We first define the Christoffel symbol in terms of the metric,   `g__&mu;,&nu;`.

``

Christoffel[`~rho`, mu, nu] = convert(Christoffel[`~rho`, mu, nu], g_)

Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-g_[`~alpha`, `~rho`]*(Physics:-d_[nu](Physics:-g_[alpha, mu], [X])+Physics:-d_[mu](Physics:-g_[alpha, nu], [X])-Physics:-d_[alpha](Physics:-g_[mu, nu], [X]))

(6)

SubstituteTensor(g_[mu, nu] = eta[mu, nu]+Physics:-`*`(l[mu](X), l[nu](X)), Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-g_[`~alpha`, `~rho`]*(Physics:-d_[nu](Physics:-g_[alpha, mu], [X])+Physics:-d_[mu](Physics:-g_[alpha, nu], [X])-Physics:-d_[alpha](Physics:-g_[mu, nu], [X])), evaluateexpression)

Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*(eta[`~alpha`, `~rho`]+l[`~alpha`](X)*l[`~rho`](X))*(Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])+Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-Physics:-d_[alpha](l[mu](X), [X])*l[nu](X)-l[mu](X)*Physics:-d_[alpha](l[nu](X), [X]))

(7)

Simplify(SubstituteTensor(Physics:-`*`(l[`~alpha`](X), l[`~rho`](X)) = -Physics:-`*`(l[`~alpha`](X), l[`~rho`](X)), Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*(eta[`~alpha`, `~rho`]+l[`~alpha`](X)*l[`~rho`](X))*(Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])+Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-Physics:-d_[alpha](l[mu](X), [X])*l[nu](X)-l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])), evaluateexpression))

Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])

(8)

SubstituteTensor(Physics:-`*`(l[alpha](X), eta[`~alpha`, `~rho`]) = l[`~rho`](X), Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X]))

Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~rho`](X)*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~rho`](X)*Physics:-d_[mu](l[nu](X), [X])+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])

(9)

SubstituteTensor(Physics:-`*`(l[`~alpha`](X), l[alpha](X)) = 0, Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~rho`](X)*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~rho`](X)*Physics:-d_[mu](l[nu](X), [X])+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[nu](l[mu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[alpha](X)*Physics:-d_[mu](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X]))

Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~rho`](X)*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~rho`](X)*Physics:-d_[mu](l[nu](X), [X])+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])

(10)

NULL

NULL

Now we can substitute into the null condition for the Ricci tensor, `R__&mu;&nu;`*`#mi("l")`^mu*l^nu = 0.

convert(Physics:-`*`(Physics:-`*`(Ricci[mu, nu], l[`~mu`](X)), l[`~nu`](X)), Christoffel)

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

(11)

NULL

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

(Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])-Physics:-d_[nu](Physics:-Christoffel[`~rho`, mu, rho], [X])+Physics:-Christoffel[`~beta`, mu, nu]*Physics:-Christoffel[`~rho`, beta, rho]-Physics:-Christoffel[`~beta`, mu, rho]*Physics:-Christoffel[`~rho`, beta, nu])*l[`~mu`](X)*l[`~nu`](X)

(12)

  Do the first term

 

expand(Physics:-`*`(Physics:-`*`(d_[rho](Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~rho`](X)*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~rho`](X)*Physics:-d_[mu](l[nu](X), [X])+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])), l[`~mu`](X)), l[`~nu`](X)))

Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = (1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[mu](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[mu](l[nu](X), [X]), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[nu](l[mu](X), [X]), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[mu](X), [X]), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[mu](l[alpha](X), [X]), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*Physics:-d_[rho](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[nu](X), [X]), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[nu](l[alpha](X), [X]), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*Physics:-d_[rho](l[mu](X), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](Physics:-d_[alpha](l[mu](X), [X]), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[alpha](l[mu](X), [X])*Physics:-d_[rho](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[nu](X)*Physics:-d_[rho](Physics:-d_[mu](l[alpha](X), [X]), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[mu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[nu](X), [X]), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[nu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[nu](l[alpha](X), [X]), [X])*eta[`~alpha`, `~rho`]

(13)

NULL

SubstituteTensor(Physics:-`*`(l[`~nu`](X), l[nu](X)) = 0, Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = (1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[mu](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[mu](l[nu](X), [X]), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[nu](l[mu](X), [X]), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[alpha](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[mu](X), [X]), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[mu](l[alpha](X), [X]), [X])*l[nu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*Physics:-d_[rho](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[nu](X), [X]), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~alpha`](X), [X])*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[nu](l[alpha](X), [X]), [X])*l[mu](X)-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*Physics:-d_[rho](l[mu](X), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](Physics:-d_[alpha](l[mu](X), [X]), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[alpha](l[mu](X), [X])*Physics:-d_[rho](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[nu](X)*Physics:-d_[rho](Physics:-d_[mu](l[alpha](X), [X]), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[mu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[alpha](l[nu](X), [X]), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[nu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[mu](X)*Physics:-d_[rho](Physics:-d_[nu](l[alpha](X), [X]), [X])*eta[`~alpha`, `~rho`])

Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = (1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[alpha](l[mu](X), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*Physics:-d_[rho](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*Physics:-d_[rho](l[mu](X), [X])-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[alpha](l[mu](X), [X])*Physics:-d_[rho](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[nu](X), [X])*Physics:-d_[mu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[mu](X), [X])*Physics:-d_[nu](l[alpha](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[mu](l[nu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[mu](l[nu](X), [X]), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*Physics:-d_[rho](l[`~rho`](X), [X])*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~mu`](X)*l[`~nu`](X)*l[`~rho`](X)*Physics:-d_[rho](Physics:-d_[nu](l[mu](X), [X]), [X])

(14)

 

 

Do same thing with the first term but use the Simplify command

 

Simplify(Physics:-`*`(Physics:-`*`(d_[rho](Physics:-Christoffel[`~rho`, mu, nu] = (1/2)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])*eta[`~alpha`, `~rho`]+(1/2)*l[`~rho`](X)*Physics:-d_[nu](l[mu](X), [X])+(1/2)*l[`~rho`](X)*Physics:-d_[mu](l[nu](X), [X])+(1/2)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)*eta[`~alpha`, `~rho`]-(1/2)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])*eta[`~alpha`, `~rho`]-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[nu](l[alpha](X), [X])*l[mu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[mu](X)*Physics:-d_[alpha](l[nu](X), [X])-(1/2)*l[`~alpha`](X)*l[`~rho`](X)*Physics:-d_[mu](l[alpha](X), [X])*l[nu](X)+(1/2)*l[`~alpha`](X)*l[`~rho`](X)*l[nu](X)*Physics:-d_[alpha](l[mu](X), [X])), l[`~mu`](X)), l[`~nu`](X)))

Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = -Physics:-d_[alpha2](l[`~alpha3`](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha5`](X)*l[`~rho`](X)*l[rho](X)+Physics:-d_[alpha4](l[`~alpha5`](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha3`](X)*l[`~alpha4`](X)-Physics:-d_[alpha](l[alpha3](X), [X])*eta[`~alpha`, `~alpha1`]*Physics:-d_[alpha1](l[alpha2](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+eta[`~alpha`, `~alpha5`]*Physics:-d_[alpha2](l[alpha](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)-eta[`~alpha5`, `~alpha3`]*Physics:-d_[alpha3](Physics:-d_[alpha5](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha6`](X)+eta[`~alpha6`, `~alpha5`]*Physics:-d_[alpha3](Physics:-d_[alpha5](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](l[alpha3](X), [X])*Physics:-d_[rho](l[`~rho`](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](Physics:-d_[alpha3](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)*l[`~alpha6`](X)

(15)

SubstituteTensor(Physics:-`*`(l[`~nu`](X), l[nu](X)) = 0, Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = -Physics:-d_[alpha2](l[`~alpha3`](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha5`](X)*l[`~rho`](X)*l[rho](X)+Physics:-d_[alpha4](l[`~alpha5`](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha3`](X)*l[`~alpha4`](X)-Physics:-d_[alpha](l[alpha3](X), [X])*eta[`~alpha`, `~alpha1`]*Physics:-d_[alpha1](l[alpha2](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+eta[`~alpha`, `~alpha5`]*Physics:-d_[alpha2](l[alpha](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)-eta[`~alpha5`, `~alpha3`]*Physics:-d_[alpha3](Physics:-d_[alpha5](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha6`](X)+eta[`~alpha6`, `~alpha5`]*Physics:-d_[alpha3](Physics:-d_[alpha5](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[alpha2](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](l[alpha3](X), [X])*Physics:-d_[rho](l[`~rho`](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](Physics:-d_[alpha3](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)*l[`~alpha6`](X))

Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = -Physics:-d_[alpha](l[alpha3](X), [X])*eta[`~alpha`, `~alpha1`]*Physics:-d_[alpha1](l[alpha2](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+eta[`~alpha`, `~alpha5`]*Physics:-d_[alpha2](l[alpha](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](l[alpha3](X), [X])*Physics:-d_[rho](l[`~rho`](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](Physics:-d_[alpha3](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)*l[`~alpha6`](X)

(16)

 

Simplify command does make the algebra easier, but the indices are not the same. Now, equation 16 should correspond to equation 14, but there is no combination of alphas that is consistent.  The variables alpha1 and alpha5 must be rho. One term is always wrong when I try to change the other indices.

alpha2 and alpha3 must be either mu or nu based on the first term.  But alpha6 should also be either mu or nu based on the last term, however that will make alpha2 and alpha3 either (nu and rho) or (mu and rho).  Neither combination makes all of the terms consistent with (14).  Very frustrating.

 

SubstituteTensorIndices({alpha1 = rho, alpha5 = rho}, Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = -Physics:-d_[alpha](l[alpha3](X), [X])*eta[`~alpha`, `~alpha1`]*Physics:-d_[alpha1](l[alpha2](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+eta[`~alpha`, `~alpha5`]*Physics:-d_[alpha2](l[alpha](X), [X])*Physics:-d_[alpha5](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](l[alpha3](X), [X])*Physics:-d_[rho](l[`~rho`](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](Physics:-d_[alpha3](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)*l[`~alpha6`](X))

Physics:-d_[rho](Physics:-Christoffel[`~rho`, mu, nu], [X])*l[`~mu`](X)*l[`~nu`](X) = -Physics:-d_[alpha](l[alpha3](X), [X])*eta[`~alpha`, `~rho`]*Physics:-d_[rho](l[alpha2](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+eta[`~alpha`, `~rho`]*Physics:-d_[alpha2](l[alpha](X), [X])*Physics:-d_[rho](l[alpha3](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](l[alpha3](X), [X])*Physics:-d_[rho](l[`~rho`](X), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)+Physics:-d_[alpha2](Physics:-d_[alpha3](l[alpha6](X), [X]), [X])*l[`~alpha2`](X)*l[`~alpha3`](X)*l[`~alpha6`](X)

(17)

``

 

``

``

``

``

``


Download Vacuum_Solutions_(Kerr-Schild)_3.mw

First 1457 1458 1459 1460 1461 1462 1463 Last Page 1459 of 2456