Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

So this is my minimal working code. Everything works, but I cannot get the arrow size fixed you can see the animation propperly. Adding wid=1/2 gives an error message.

dsolve.mw

Hi all, I want the bewst for you.

 

Could anyone help me with this bad equation, please?

 

Regards.

This message is for those who prefer use Maple in 1-D Math Input.

In 1-D Math Input, in previous versions of Maple, it was very comfortable to have the freedom to write ONE large statement in Continuing ON SEVERAL LINES for clarity and better reading from a *.wm file.

Here are simple examples, but in reality, I work on very complex cases.

 

Example #1:

A := u*( x + ln(x +1) + cos(x))

     + v*(1 + sqrt(x))

     + w*(sin(x) + tan(x) + x);

 

Example #2:

Vector([cos(s)*cos(t)

          , cos(s)*sin(t)

          , sin(s)]);

 

Example #3:

display(

  plot(f(x), x = x1..x2)

, plot(g(x), x = x1..x2)

, plot(h(x), x = x1..x2)

);

 

In Maple 2015 with a *.wm file, when you try to execute these example in 1-D Math Input, an error is returned and unfortunately forces you to write everything on the same line, what makes reading tiring.

 

Is a bug or a voluntary deactivation?

Can you help me, please?

 

Guy.

I have drawn a 3d plot of a sphere in maple, but when I try to export this one in an .eps file, only the coordinate system of the plot is shown, but not the colourful plot. The export works if I use jpg and similar file types, but not with .eps, does anybody know if there is a way out of this problem?

Hi,

Is there a way to export worksheet to pdf format but not in A4 size since some lines are going beyond the page?

I have a problem solving a system of PDEs.

The system of PDEs are

PDE01 := -(l^2+1^2)*(diff(v(l, t), t))+(l^2+1^2)*(diff(R(l, t), l, l))+4*l*(diff(R(l, t), l))+4*l*v(l, t)/(l^2+1^2)^(1/4)-6*R(l, t)/(l^2+1^2)+(l^2+1^2)^(1/2)*(-1.1+sqrt(.1))^2*sqrt(24)*u(l, t) = 0

PDE02 := diff(R(l, t), t) = v(l, t)

PDE03 := diff(u(l, t), t)-sqrt((1.1^2-1)/1.1^2)*(diff(u(l, t), l))-2*l*sqrt(1.1^2-1)*u(l, t)/(l^2+1^2) = 0

the initial condisions are

v(l, 0) = 0, R(l, 0) = 0, u(l, 0) = sqrt((l^2+1^2)^(1/2))*10^(-5)*exp(-(l-10)^2/.5^2)

and the BCs are

bdry00 := {((30^2+1^2)/30^2)^(1/4)*v(-30, t) = -((30^2+1^2)/30^2)^(1/2)*(D[1](R))(-30, t), ((30^2+1^2)/30^2)^(1/4)*v(30, t) = -((30^2+1^2)/30^2)^(1/2)*(D[1](R))(30, t), u(-30, t) = sqrt(30^2+1^2)*10^(-5)*exp(-40000), u(30, t) = sqrt(30^2+1^2)*10^(-5)*exp(-10000)}

to solve the system,

I enter

pde := pdsolve({PDE01, PDE02, PDE03}, {bdry00, init00}, time = t, numeric, range = -30 .. 30, timesstep = 1/60, spaceste = 1/254)

then, I failed to get the result constantly.

I tried several cases changing the initial conditions...

Can you let me know what I am doing wrong?

 

I have a problem with IsMatrixShape. I have in my part of formulation this matrix expression: QTIbQ

While Ib is a symmetric matrix, this matrix expression is clearly symmetric. However, when I try to check this issue with IsMatrixShape command, it returns false. I am extremely confused. Can anyone help me? Thanks in advance.

prob1.mw    prob2.mw

Hi, I don't know why this happened...

If I write the command in new sheet, it runs correctly( prob2.mw).

 

Second, help me to evaluate N, please...  (prob1.mw)

 

Regards :-)

 

 

 

 

 

Another application for the study of rational numbers in operations, generating fraction, etc.

 

Numeros_Racionales.mw

(in spanish)

 

Atte.

L.AraujoC.


restart; with(Physics); with(Tetrads); with(PDETools)

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

(1)

coords := zetabar, zeta, v, u

zetabar, zeta, v, u

(2)

X = [coords]

X = [zetabar, zeta, v, u]

(3)

ds2 := Physics:-`*`(Physics:-`*`(2, dzeta), dzetabar)+Physics:-`*`(Physics:-`*`(2, du), dv)+Physics:-`*`(Physics:-`*`(2, H(coords)), Physics:-`^`(du+Physics:-`*`(Ybar(coords), dzeta)+Physics:-`*`(Y(coords), dzetabar)-Physics:-`*`(Physics:-`*`(Y(coords), Ybar(coords)), 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

(4)

PDEtools:-declare(ds2)

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

(5)

vierbien := Matrix([[1, 0, -Ybar(coords), 0], [0, 1, -Y(coords), 0], [Physics:-`*`(H(coords), Y(coords)), Physics:-`*`(H(coords), Ybar(coords)), 1-Physics:-`*`(Physics:-`*`(H(coords), Y(coords)), Ybar(coords)), H(coords)], [Y(coords), Ybar(coords), -Physics:-`*`(Y(coords), Ybar(coords)), 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})

(6)

``

Physics:-Setup(coordinatesystem = (X = [zetabar, zeta, v, u]), metric = ds2, tetrad = vierbien, mathematicalnotation = true, automaticsimplification = true, signature = "+++-")

RicciT := proc (a, b) options operator, arrow; Physics:-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

(7)

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

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

(8)

SlashD(H(X), 4) = H(X)[4]

(diff(H(X), zetabar))*Ybar(X)+(diff(H(X), zeta))*Y(X)+diff(H(X), v)-(diff(H(X), u))*Y(X)*Ybar(X) = H(X)[4]

(9)

Gamma := proc (a, b, c) options operator, arrow; -gamma_[a, b, c] end proc

proc (a, b, c) options operator, arrow; Physics:-`*`(Physics:-Tetrads:-gamma_[a, b, c], -1) end proc

(10)

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

(11)

``

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

(12)

``

Gamma(3, 4, 4) = 0

0 = 0

(13)

``

Gamma(4, 4, 4) = 0

0 = 0

(14)

Gamma(4, 1, 1) = 0

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

(15)

``

Gamma(4, 2, 2) = 0

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

(16)

NULL

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}

{-(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}

(17)

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), zetabar))*Ybar(X)+(diff(H(X), zeta))*Y(X)+diff(H(X), v)-(diff(H(X), u))*Y(X)*Ybar(X))*(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)

0 = 0

0 = 0

(19)

The values in the paraenthesis  should substitute H[4]. This sequence works in Maple 18 but not in Maple 2015

 

NULL


Download Question_algsubs_3.27.15.mw

I'm attempting to plot the Nyquist plot for a complex system such that Maple cannot determine the frequency limits automatically, and suggested I use 'range' to specify the frequency limits. How do you do this? For example if you have a system G=1/(s+1), how do you plot the Nyquist plot in the frequency range 0 to 1rad/s?  Thanks

Helo friends. Hope you will be fine. I need the command for take the conjugate of exp(I*x). I need the result exp(-I*x), x treated as real number.

 

PhD (Scholar)
Department of Mathematics

In a plot command, how do I indicate not to show the legend ?

(I know about unchecking "show legend" with a right click, and I have searched mapleprimes.)

 

Thanks, all !

Write a Maple code that performs the Gaussian elimination for an nxn matrix, converting it to an upper triangular matrix. 

(Hint: you will need to use three for .. do loops.)

How to add values in vector? This is what I am trying in my code below. I know there are more than one values satifying my conditon.

What is the best way constructing such vectors. If there are other options, please let me know. Thanks.

First 1286 1287 1288 1289 1290 1291 1292 Last Page 1288 of 2256