Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi,

I am new to maple. I want to write the following codes in Maple in a file and execute it. I am writing the code in matlab. thanks

 

for i =1:1:5

for j=1:1:5

M(i,j)=i+j;

end

end

 

thanks

I am trying to take a list of prime numbers that I already generated with the following:

L4:=NULL:

for i from 2 to 50 do

if nops(ifactor(i))=1 then L4:=L4,i end if:

end do;

[L4]

This gives me the primes bewteen 2-50, [2,3,5,7,11,13,17,23,29,31,37,43,47].  I want to write something that takes the products of my list like this

2*3=6

2*3*5=30

2*3*5*7=210

and so on through the whole list.  Any ideas?  I figured out how to take the product of 2*3, 2*5, 2*7 and so on, but I cant figure how to make it do 2*3, 2*3*7 and so on.

I am doing a presentation on Euclids infinite prime proof and thought it would be cool to write a maple program that will show his proof for a certain amount of prime numbers.

 

Thanks,  Matt

 

 

 

Hi,

When I show the context menu on maple 18 (CTRL click on mac or right click with a mouse) on an equation the system doesn't show the tutor or series options as written in the example I use.   (after having loaded the multivariate calculus as asked)

Is that an issue in my installation ?  Where do I have to place the cursor ?  Is the multivariate calculus module well loaded ?

Thanks for your help

 

Stéphane

Dear friends,

I would like to know the agorithm of isovle for solving 2 order equations for more than 3 variables.

I see that not all the solutions are shown for example for:

x^2+y^2=z^2+t^2

 Could you please help me with this?

 

hi,

I have bought maple18 student edition. I want to learn GPU programming through Maple. Please suggest how to do this. I have a notebook with i7 processor and NVIDIA geforce 750m graphics. I want to solve system of algebraic equations, integral equations etc in parallel using GPUs.

thanks

I am trying to simplify sums of a few LaguerreL polinomials of different n using the identities in the function advisor such as recurrsion relations. How does one go about in using the FunctionAdvisor identities when trying to simplify expressions containing orthogonal polynomials? 

 

 

I want to generate all 4x4 matrices of rank-1 over the field with two elements, F_2. I don't want to use actual arrays though, I just want to use symbols representing the unit matrices. So for example Eij represents the 4x4 matrix with a 1 in row i & column j, and 0's elsewhere. 

We have 16 unit matrices: E11, E12, ..., E14, ... E41, E42, E43, E44. Then we have to look at all possible linear combinations that result in a rank-1 matrix. So for example: E11+E12,  E11+E21,  E11+E12+E13,  E33+E43, etc. The total number of rank-1 4x4 matrices over F_2 is 225. How can I find all of these quickly?

Hi number of officials,

As I enter the command defined 'int' I get an error message.

What can I do. Can you help me.

 

Commands Entered:

int(x^2+3*x, x);

 

ansvers:

Error, (in int) wrong number (or type) of arguments: invalid option value passed to indefinite integration: {}

I'm wondering why the fit routine gives weird results, when telling him to fit

with(Statistics):

X:=Vector[column](3, {1 = 616.3096, 2 = 616.3248, 3 = 616.3342}, datatype = sfloat, storage = rectangular, order = Fortran_order, shape = []);

Y:=Vector[column](3, {1 = 99590., 2 = 127564., 3 = 144223.}, datatype = sfloat, storage = rectangular, order = Fortran_order, shape = []);

g:=Fit(a+b*x+c*x^2,X,Y,x)

It's telling me: Warning, model is not of full rank

If I calculate the coefficients the direct way (well the fit should be exact,because we only have 3 values and 3 parameters) I get:

M:=Matrix(3, 3, {(1, 1) = 1, (1, 2) = 616.3096, (1, 3) = 379837.52305216, (2, 1) = 1, (2, 2) = 616.3248, (2, 3) = 379856.25909504, (3, 1) = 1, (3, 2) = 616.3342, (3, 3) = 379867.84608964}, datatype = anything, storage = rectangular, order = Fortran_order, shape = []);
V:=Vector[column](3, {1 = 99590, 2 = 127564, 3 = 144223}, datatype = anything, storage = rectangular, order = Fortran_order, shape = []);
with(LinearAlgebra):
C:=MatrixVectorMultiply(M^-1,V)

 

Hi there

I'm trying to isolate (y1-3)2+(x1-1)in the equation 25(y1-3)2+200+100(x1-1)2=0.

I have tried isolate and solve, but solve coplains about solving for expressions (but when inputting i:=(x1,y1)->(y1-3)2+(x1-1)2 it still doesn't work), and isolate can only isolate either (y1-3)2 or (x1-1). Not both.

How can I do this with as few lines as possible?

Thanks

- Alex

Use the Maple command DEplot to draw a direction field for the equation

make sure it at least covers the area −1 ≤ x ≤ 1,
−1 ≤ y ≤ 1

the equation is dy/dx+4y^3 −3y = 0.

then I have :

y(x):=dy/(dx)+4*y^(3)-3*y=0

with(DEtools); DEplot(y(x), x = -1 .. .1, y = -1 .. .1);

Error, (in DEtools/DEplot) vars must be declared as a list, e.g. [x(t),y(t),...]

 

so how can I change the command to make it work..

Hello friends,

I want to derive the function of a curve. It's easy in for example "xyExtract" but I don't know how to do it in maple. please help. thank you all.

Hi, I'm new to Maple and was trying to use it to solve 3 equations with 3 unknowns, in terms of another 2 parameters. This is what I put in, and the error that came up:

solve({(1-x)/(1+b) = b*y, (1-y)(1-a)/(-a*b+1) = b*z, (1-z)(1-b)/(-a*b+1) = a*x}, {x, y, z}); 

Error, (in SolveTools:-LinearSolvers:-Algebraic) unable to compute coeff

I want to solve the equations to get x,y,z in terms of a,b but I don't understand the error coming up - have I done something wrong or is it because not all the variables appear in all of the equations? As there are 3 equations and 3 unknowns there is a solution, but I want to check the answer I found on paper with something (the algebra got a bit messy!)

Any help greatly appreciated! :)

I am trying to simplify equation 18 using equations 8 and 9. It should look a little like equation 21, but instead I get the results in equations 19 and 20.  I tried using different substituions, but algsubs gets the closest answer. A few terms are going to zero after the substitution.

When I substitute Z(X) then Zbar(X) terms vanish, and visa versa.


Initialize the metric and tetrad

 

restart; with(Physics); with(Tetrads)

0, "%1 is not a command in the %2 package", Tetrads, Physics

(1.1)

X = [zetabar, zeta, v, u]

X = [zetabar, zeta, v, u]

(1.2)

ds2 := Physics:-`*`(Physics:-`*`(2, dzeta), dzetabar)+Physics:-`*`(Physics:-`*`(2, du), dv)+Physics:-`*`(Physics:-`*`(2, H(zetabar, zeta, v, u)), (du+Physics:-`*`(Ybar(zetabar, zeta, v, u), dzeta)+Physics:-`*`(Y(zetabar, zeta, v, u), dzetabar)-Physics:-`*`(Physics:-`*`(Y(zetabar, zeta, v, u), Ybar(zetabar, zeta, v, u)), dv))^2)

2*dzeta*dzetabar+2*du*dv+2*H(zetabar, zeta, v, u)*(du+Ybar(zetabar, zeta, v, u)*dzeta+Y(zetabar, zeta, v, u)*dzetabar-Y(zetabar, zeta, v, u)*Ybar(zetabar, zeta, v, u)*dv)^2

(1.3)

PDEtools:-declare(ds2)

Ybar(zetabar, zeta, v, u)*`will now be displayed as`*Ybar

(1.4)

NULL

vierbien = Matrix([[1, 0, -Ybar(zetabar, zeta, v, u), 0], [0, 1, -Y(zetabar, zeta, v, u), 0], [Physics:-`*`(H(zetabar, zeta, v, u), Y(zetabar, zeta, v, u)), Physics:-`*`(H(zetabar, zeta, v, u), Ybar(zetabar, zeta, v, u)), 1-Physics:-`*`(Physics:-`*`(H(zetabar, zeta, v, u), Y(zetabar, zeta, v, u)), Ybar(zetabar, zeta, v, u)), H(zetabar, zeta, v, u)], [Y(zetabar, zeta, v, u), Ybar(zetabar, zeta, v, u), -Physics:-`*`(Y(zetabar, zeta, v, u), Ybar(zetabar, zeta, v, u)), 1]])

vierbien = (Matrix(4, 4, {(1, 1) = 1, (1, 2) = 0, (1, 3) = -Ybar(zetabar, Zeta, v, u), (1, 4) = 0, (2, 1) = 0, (2, 2) = 1, (2, 3) = -Y(zetabar, Zeta, v, u), (2, 4) = 0, (3, 1) = H(zetabar, Zeta, v, u)*Y(zetabar, Zeta, v, u), (3, 2) = H(zetabar, Zeta, v, u)*Ybar(zetabar, Zeta, v, u), (3, 3) = 1-H(zetabar, Zeta, v, u)*Y(zetabar, Zeta, v, u)*Ybar(zetabar, Zeta, v, u), (3, 4) = H(zetabar, Zeta, v, u), (4, 1) = Y(zetabar, Zeta, v, u), (4, 2) = Ybar(zetabar, Zeta, v, u), (4, 3) = -Y(zetabar, Zeta, v, u)*Ybar(zetabar, Zeta, v, u), (4, 4) = 1}))

(1.5)

``

NULL

Setup(tetrad = rhs(vierbien = Matrix(%id = 18446744078408794830)), metric = ds2, mathematicalnotation = true, automaticsimplification = true, coordinatesystems = (X = [zetabar, zeta, v, u]), signature = "+++-")

[automaticsimplification = true, coordinatesystems = {X}, mathematicalnotation = true, metric = {(1, 1) = 2*H(X)*Y(X)^2, (1, 2) = 1+2*H(X)*Y(X)*Ybar(X), (1, 3) = -2*H(X)*Y(X)^2*Ybar(X), (1, 4) = 2*H(X)*Y(X), (2, 2) = 2*H(X)*Ybar(X)^2, (2, 3) = -2*H(X)*Ybar(X)^2*Y(X), (2, 4) = 2*H(X)*Ybar(X), (3, 3) = 2*H(X)*Y(X)^2*Ybar(X)^2, (3, 4) = 1-2*H(X)*Y(X)*Ybar(X), (4, 4) = 2*H(X)}, signature = `+ + + -`, tetrad = {(1, 1) = 1, (1, 3) = -Ybar(X), (2, 2) = 1, (2, 3) = -Y(X), (3, 1) = H(X)*Y(X), (3, 2) = H(X)*Ybar(X), (3, 3) = 1-H(X)*Y(X)*Ybar(X), (3, 4) = H(X), (4, 1) = Y(X), (4, 2) = Ybar(X), (4, 3) = -Y(X)*Ybar(X), (4, 4) = 1}]

(1.6)

``

Verification of Tetrad

 

I will try to verify the tetrad from (Kerr and Schild (1965)). However, the tetrad given in the paper seems to have the third tetrad with the wrong sign. I changed the sign and get the correct verification,

    e_[]

`𝔢`[a, mu] = (Matrix(4, 4, {(1, 1) = 1, (1, 2) = 0, (1, 3) = -Ybar(X), (1, 4) = 0, (2, 1) = 0, (2, 2) = 1, (2, 3) = -Y(X), (2, 4) = 0, (3, 1) = H(X)*Y(X), (3, 2) = H(X)*Ybar(X), (3, 3) = 1-H(X)*Y(X)*Ybar(X), (3, 4) = H(X), (4, 1) = Y(X), (4, 2) = Ybar(X), (4, 3) = -Y(X)*Ybar(X), (4, 4) = 1}))

(2.1)

g_[]

g[mu, nu] = (Matrix(4, 4, {(1, 1) = 2*H(X)*Y(X)^2, (1, 2) = 1+2*H(X)*Y(X)*Ybar(X), (1, 3) = -2*H(X)*Y(X)^2*Ybar(X), (1, 4) = 2*H(X)*Y(X), (2, 1) = 1+2*H(X)*Y(X)*Ybar(X), (2, 2) = 2*H(X)*Ybar(X)^2, (2, 3) = -2*H(X)*Ybar(X)^2*Y(X), (2, 4) = 2*H(X)*Ybar(X), (3, 1) = -2*H(X)*Y(X)^2*Ybar(X), (3, 2) = -2*H(X)*Ybar(X)^2*Y(X), (3, 3) = 2*H(X)*Y(X)^2*Ybar(X)^2, (3, 4) = 1-2*H(X)*Y(X)*Ybar(X), (4, 1) = 2*H(X)*Y(X), (4, 2) = 2*H(X)*Ybar(X), (4, 3) = 1-2*H(X)*Y(X)*Ybar(X), (4, 4) = 2*H(X)}))

(2.2)

Physics:-`*`(e_[a, mu], e_[a, nu]) = g_[mu, nu]

Physics:-Tetrads:-e_[a, mu]*Physics:-Tetrads:-e_[`~a`, nu] = Physics:-g_[mu, nu]

(2.3)

TensorArray(Physics:-Tetrads:-e_[a, mu]*Physics:-Tetrads:-e_[`~a`, nu] = Physics:-g_[mu, nu])

Matrix(4, 4, {(1, 1) = 2*H(X)*Y(X)^2 = 2*H(X)*Y(X)^2, (1, 2) = 1+2*H(X)*Y(X)*Ybar(X) = 1+2*H(X)*Y(X)*Ybar(X), (1, 3) = -2*H(X)*Y(X)^2*Ybar(X) = -2*H(X)*Y(X)^2*Ybar(X), (1, 4) = 2*H(X)*Y(X) = 2*H(X)*Y(X), (2, 1) = 1+2*H(X)*Y(X)*Ybar(X) = 1+2*H(X)*Y(X)*Ybar(X), (2, 2) = 2*H(X)*Ybar(X)^2 = 2*H(X)*Ybar(X)^2, (2, 3) = -2*H(X)*Ybar(X)^2*Y(X) = -2*H(X)*Ybar(X)^2*Y(X), (2, 4) = 2*H(X)*Ybar(X) = 2*H(X)*Ybar(X), (3, 1) = -2*H(X)*Y(X)^2*Ybar(X) = -2*H(X)*Y(X)^2*Ybar(X), (3, 2) = -2*H(X)*Ybar(X)^2*Y(X) = -2*H(X)*Ybar(X)^2*Y(X), (3, 3) = 2*H(X)*Y(X)^2*Ybar(X)^2 = 2*H(X)*Y(X)^2*Ybar(X)^2, (3, 4) = 1-2*H(X)*Y(X)*Ybar(X) = 1-2*H(X)*Y(X)*Ybar(X), (4, 1) = 2*H(X)*Y(X) = 2*H(X)*Y(X), (4, 2) = 2*H(X)*Ybar(X) = 2*H(X)*Ybar(X), (4, 3) = 1-2*H(X)*Y(X)*Ybar(X) = 1-2*H(X)*Y(X)*Ybar(X), (4, 4) = 2*H(X) = 2*H(X)})

(2.4)

Physics:-`*`(e_[a, mu], e_[b, mu]) = eta_[a, b]

Physics:-Tetrads:-e_[a, mu]*Physics:-Tetrads:-e_[b, `~mu`] = Physics:-Tetrads:-eta_[a, b]

(2.5)

NULL

TensorArray(Physics:-Tetrads:-e_[a, mu]*Physics:-Tetrads:-e_[b, `~mu`] = Physics:-Tetrads:-eta_[a, b])

Matrix(4, 4, {(1, 1) = 0 = 0, (1, 2) = 1 = 1, (1, 3) = 0 = 0, (1, 4) = 0 = 0, (2, 1) = 1 = 1, (2, 2) = 0 = 0, (2, 3) = 0 = 0, (2, 4) = 0 = 0, (3, 1) = 0 = 0, (3, 2) = 0 = 0, (3, 3) = 0 = 0, (3, 4) = 1 = 1, (4, 1) = 0 = 0, (4, 2) = 0 = 0, (4, 3) = 1 = 1, (4, 4) = 0 = 0})

(2.6)

``

gamma_[4, 2, 1]

diff(Y(X), zeta)-(diff(Y(X), u))*Ybar(X)

(2.7)

SumOverRepeatedIndices(Physics:-`*`(Physics:-`*`(D_[nu](e_[4, mu]), e_[2, mu]), e_[1, `~nu`]))

diff(Y(X), zeta)-(diff(Y(X), u))*Ybar(X)

(2.8)

NULL

``

For equation 2.8 we get the following:

SumOverRepeatedIndices(Physics:-`*`(Physics:-`*`(Riemann[`~sigma`, rho, mu, nu], e_[4, `~rho`]), e_[4, `~nu`]))

(-Physics:-Riemann[`~sigma`, 4, 4, mu]*Y(X)^2+(Physics:-Riemann[`~sigma`, 4, 1, mu]+Physics:-Riemann[`~sigma`, 1, 4, mu])*Y(X)-Physics:-Riemann[`~sigma`, 1, 1, mu])*Ybar(X)^2+((Physics:-Riemann[`~sigma`, 2, 4, mu]+Physics:-Riemann[`~sigma`, 4, 2, mu])*Y(X)^2+(-Physics:-Riemann[`~sigma`, 2, 1, mu]+Physics:-Riemann[`~sigma`, 3, 4, mu]+Physics:-Riemann[`~sigma`, 4, 3, mu]-Physics:-Riemann[`~sigma`, 1, 2, mu])*Y(X)-Physics:-Riemann[`~sigma`, 3, 1, mu]-Physics:-Riemann[`~sigma`, 1, 3, mu])*Ybar(X)-Physics:-Riemann[`~sigma`, 2, 2, mu]*Y(X)^2+(-Physics:-Riemann[`~sigma`, 2, 3, mu]-Physics:-Riemann[`~sigma`, 3, 2, mu])*Y(X)-Physics:-Riemann[`~sigma`, 3, 3, mu]

(1)

 

Now we replicate eqn 2.16. These are the conditions for e[4,mu] to be geodesic and shear-free. The outputs are eqn 3.5.

 

gamma_[4, 1, 1] = 0

diff(Ybar(X), zeta)-(diff(Ybar(X), u))*Ybar(X) = 0

(2)

gamma_[4, 2, 2] = 0

diff(Y(X), zetabar)-(diff(Y(X), u))*Y(X) = 0

(3)

gamma_[1, 4, 4] = 0

(diff(Ybar(X), u))*Y(X)*Ybar(X)-Y(X)*(diff(Ybar(X), zeta))-Ybar(X)*(diff(Ybar(X), zetabar))-(diff(Ybar(X), v)) = 0

(4)

gamma_[2, 4, 4] = 0

(diff(Y(X), u))*Y(X)*Ybar(X)-Y(X)*(diff(Y(X), zeta))-(diff(Y(X), zetabar))*Ybar(X)-(diff(Y(X), v)) = 0

(5)

gamma_[3, 4, 4] = 0

0 = 0

(6)

gamma_[4, 4, 4] = 0

0 = 0

(7)

shearconditions := {diff(Y(X), zetabar)-(diff(Y(X), u))*Y(X) = 0, diff(Ybar(X), zeta)-(diff(Ybar(X), u))*Ybar(X) = 0, (diff(Y(X), u))*Y(X)*Ybar(X)-Y(X)*(diff(Y(X), zeta))-(diff(Y(X), zetabar))*Ybar(X)-(diff(Y(X), v)) = 0, (diff(Ybar(X), u))*Y(X)*Ybar(X)-Y(X)*(diff(Ybar(X), zeta))-Ybar(X)*(diff(Ybar(X), zetabar))-(diff(Ybar(X), v)) = 0}:

 

Now we can define the rotation coefficients associated with rotation and expansion z = theta - i omega

 

gamma_[2, 4, 1] = Z(X)

-(diff(Y(X), zeta))+(diff(Y(X), u))*Ybar(X) = Z(X)

(8)

gamma_[1, 4, 2] = Zbar(X)

-(diff(Ybar(X), zetabar))+(diff(Ybar(X), u))*Y(X) = Zbar(X)

(9)

PDEtools:-declare(Z(X), Zbar(X))

Zbar(zetabar, zeta, v, u)*`will now be displayed as`*Zbar

(10)

Zdefinitions := {-(diff(Y(X), zeta))+(diff(Y(X), u))*Ybar(X) = Z(X), -(diff(Ybar(X), zetabar))+(diff(Ybar(X), u))*Y(X) = Zbar(X)}

{-(diff(Y(X), zeta))+(diff(Y(X), u))*Ybar(X) = Z(X), -(diff(Ybar(X), zetabar))+(diff(Ybar(X), u))*Y(X) = Zbar(X)}

(11)

We now show that the tetrad vectors are propogated parallel along each curve of the congruence of null geodesics which have e[4,~mu] as tangents.

 

   

We now use the tetrad form of the Ricci tensor. In order to use this in Maple we need to create a Ricci Tensor Tetrad function.

 

RicciT := proc (a, b) options operator, arrow; SumOverRepeatedIndices(Ricci[mu, nu]*e_[a, `~mu`]*e_[b, `~nu`]) end proc

proc (a, b) options operator, arrow; Physics:-SumOverRepeatedIndices(Physics:-`*`(Physics:-`*`(Physics:-Ricci[mu, nu], Physics:-Tetrads:-e_[a, `~mu`]), Physics:-Tetrads:-e_[b, `~nu`])) end proc

(12)

SlashD := proc (f, a) options operator, arrow; SumOverRepeatedIndices(D_[b](f)*e_[a, `~b`]) end proc

proc (f, a) options operator, arrow; Physics:-SumOverRepeatedIndices(Physics:-`*`(Physics:-D_[b](f), Physics:-Tetrads:-e_[a, `~b`])) end proc

(13)

SlashD(f(X), 1)

diff(f(X), zeta)-Ybar(X)*(diff(f(X), u))

(14)

SlashD(f(X), 2)

diff(f(X), zetabar)-Y(X)*(diff(f(X), u))

(15)

SlashD(f(X), 3)

(1+H(X)*Y(X)*Ybar(X))*(diff(f(X), u))-H(X)*((diff(f(X), zeta))*Y(X)+Ybar(X)*(diff(f(X), zetabar))+diff(f(X), v))

(16)

SlashD(f(X), 4)

-Y(X)*Ybar(X)*(diff(f(X), u))+(diff(f(X), zeta))*Y(X)+Ybar(X)*(diff(f(X), zetabar))+diff(f(X), v)

(17)

NULL

The geodesic and shear free condition given by Lemma 1 in (Goldberg and Sachs (1962)). Kerr uses the fourth tetrad instead of the third so we need to modify the Ricci tensor conditions. The equations (2) - (5) enforce the first Lemma.

 

   

 

Notice that none of the previous Ricci conditions can be used to solve for H.  We can use the remaining field equations to find the partial differential equations necessary to derive the metric.

 

  simplify(RicciT(1, 2), shearconditions) = 0

H(X)*(diff(diff(Y(X), zeta), zetabar))*Ybar(X)-H(X)*Ybar(X)*Y(X)*(diff(diff(Ybar(X), u), zetabar))-H(X)*Ybar(X)^2*(diff(diff(Y(X), u), zetabar))-H(X)*Y(X)^2*(diff(diff(Ybar(X), u), zeta))-2*H(X)*Y(X)*Ybar(X)*(diff(diff(Y(X), u), zeta))+H(X)*Y(X)^2*Ybar(X)*(diff(diff(Ybar(X), u), u))-H(X)*Y(X)*(diff(diff(Ybar(X), u), v))+H(X)*Y(X)*Ybar(X)^2*(diff(diff(Y(X), u), u))-H(X)*(diff(diff(Y(X), u), v))*Ybar(X)+H(X)*(diff(Ybar(X), zetabar))^2+(-3*H(X)*Y(X)*(diff(Ybar(X), u))-(diff(H(X), u))*Y(X)*Ybar(X)+(diff(H(X), zeta))*Y(X)+(diff(H(X), zetabar))*Ybar(X)+diff(H(X), v))*(diff(Ybar(X), zetabar))+H(X)*(diff(Y(X), zeta))^2+(-4*H(X)*(diff(Y(X), u))*Ybar(X)-(diff(H(X), u))*Y(X)*Ybar(X)+(diff(H(X), zeta))*Y(X)+(diff(H(X), zetabar))*Ybar(X)+diff(H(X), v))*(diff(Y(X), zeta))+2*H(X)*Y(X)^2*(diff(Ybar(X), u))^2-Y(X)*(-(diff(H(X), u))*Y(X)*Ybar(X)+(diff(H(X), zeta))*Y(X)+(diff(H(X), zetabar))*Ybar(X)+diff(H(X), v))*(diff(Ybar(X), u))+2*(H(X)*(diff(Y(X), u))*Ybar(X)+(1/2)*(diff(H(X), u))*Y(X)*Ybar(X)-(1/2)*(diff(H(X), zeta))*Y(X)-(1/2)*(diff(H(X), zetabar))*Ybar(X)-(1/2)*(diff(H(X), v)))*(diff(Y(X), u))*Ybar(X) = 0

(18)

-(diff(H(X), zetabar))*Ybar(X)*Z(X)-Y(X)*(diff(H(X), zeta))*Z(X)-H(X)*(diff(Y(X), zeta))^2+Z(X)*((diff(H(X), u))*Y(X)*Ybar(X)+2*H(X)*Z(X)-(diff(H(X), v))) = 0

-(diff(H(X), zetabar))*Ybar(X)*Z(X)-(diff(H(X), zeta))*Y(X)*Z(X)-H(X)*(diff(Y(X), zeta))^2-Z(X)*(-(diff(H(X), u))*Y(X)*Ybar(X)-2*H(X)*Z(X)+diff(H(X), v)) = 0

(19)

Zbar(X)*(-(diff(H(X), v))-(diff(H(X), zetabar))*Ybar(X)-(diff(H(X), zeta))*Y(X)+(diff(H(X), u))*Y(X)*Ybar(X)+H(X)*(diff(Ybar(X), zetabar)+2*Zbar(X))) = 0

-Zbar(X)*(-(diff(H(X), u))*Y(X)*Ybar(X)+(diff(H(X), zeta))*Y(X)+(diff(H(X), zetabar))*Ybar(X)-H(X)*(diff(Ybar(X), zetabar))-2*H(X)*Zbar(X)+diff(H(X), v)) = 0

(20)

Physics:-`*`(SlashD(H(X), 4), Z(X)+Zbar(X)) = Physics:-`*`(H(X), SlashD(Z(X), 4)+SlashD(Zbar(X), 4))

-(-(diff(H(X), v))-(diff(H(X), zeta))*Y(X)+Ybar(X)*((diff(H(X), u))*Y(X)-(diff(H(X), zetabar))))*(Z(X)+Zbar(X)) = H(X)*(-Y(X)*Ybar(X)*(diff(Z(X), u))+Ybar(X)*(diff(Z(X), zetabar))+Y(X)*(diff(Z(X), zeta))+diff(Z(X), v)-Y(X)*Ybar(X)*(diff(Zbar(X), u))+Ybar(X)*(diff(Zbar(X), zetabar))+Y(X)*(diff(Zbar(X), zeta))+diff(Zbar(X), v))

(21)

``

NULL

NULL


Download Deriving_the_Kerr_Metric.mw

Hi,

 

I'm trying to create interactive plots by using Explore to help demonstrate the effects parameters have on functions. I created one successfully to illustrate shifts and stretches of a polynomial:

 

transform(A,B,X,H,P,K):=Explore(plot(a*(b*x+h)^(p)+k,x=X),parameters=[a=A, b= B,h=H,p=P,k=K],placement=right)

 

However when I try to do the same with a solved ODE it returns an error message:

 

Explore(plot(1/(-p*x+x+1)^(1/(p-1)), x = -5 .. 5), parameters = [p = -20 .. 20], placement = right);

 

Executing this gives the error message: 

Warning, expecting only range variable x in expression 1/((-p*x+x+1)^(1/(p-1))) to be plotted but found name p
INTERFACE_PLOT(AXESLABELS(x, ""),

VIEW(-5. .. 5., DEFAULT, _ATTRIBUTE("source" = "mathdefault"))),

parameters = [p = -20 .. 20], placement = right

 

I'm not sure why it is having difficulty dealing with "p" when it had no difficulty with the first. Any help would be appreciated!

First 1267 1268 1269 1270 1271 1272 1273 Last Page 1269 of 2224