Maple Questions and Posts

These are Posts and Questions associated with the product, Maple


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.

Write a recursive Maple procedure, called “decToBin”, that converts an integer from decimal (base 10) to binary (base 2). Ensure that only integers are allowed as arguments to the procedure. Passing a negative integer should result in no output. Test your procedure by outputting the result of the following:
decToBin(10)
decToBin(163)
decToBin(213)
decToBin(-1)
decToBin(“A”) # should result in an “invalid input” error

the calculation is like the following command, the result in the picture

restart;
with(VectorCalculus);
SetCoordinates(spherical[r, theta, phi]);
Fv := rho*VectorField(`<,>`(v[r](r, theta, phi), v[theta](r, theta, phi), v[phi](r, theta, phi)));
Divergence();
Divergence(Fv);

divergence

 

1) when the Divergence act on the Fv, then it will be expanded, which is lengthy and not like most book's formulation , especially when I want to continue for a Conversation law like in fluid mechanics, this will be too long and a messy for later check.

could there be a way to not expand this result, just as the eq(3) like.

2) when I want to calculate the Divergence of Fv, I must construct a VectorField at first, but this is in components way, is there a quick way for Vector Field Function

 

 A'(t)=20+0.03*A(t)-0.00005*A(t)*Q(t)-0.0002*A(t)*Q(t)-0.0743*A(t)-12,

B'(t)=0.03*B(t)+0.00005*A(t)*Q(t)+0.0002*A(t)*Q(t)-0.0743*B(t)-0.01072*B(t),

C'(t)=0.01072*B(t)-(0.0743+0.0123)*C(t),

P'(t)=15+0.03*P(t)-0.0001*P(t)*B(t)-0.0743*A(t)-10,

Q'(t)=0.03*Q(t)+0.0001*P(t)*B(t)-0.0743*Q(t)-0.25*Q(t),

R'(t)=0.25*Q(t)-(0.0743+0.0123)*R(t),

X'(t)=25-0.0003*X(t)*B(t)-0.0743*X(t)-10,

Y'(t)=0.0003*X(t)*B(t)-0.0743*Y(t)-0.3*Y(t),

Z'(t)=0.3*Y(t)-(0.0743+0.0123)*Z(t).

Initial conditions

A(0)=150,B(0)=30, C(0)=10,P(0)=130,Q(0)=45,R(0)=10,X(0)=85,Y(0)=30, Z(0)=20.

Can I plot the graph between {A(t),B(t),C(t),P(t),Q(t),R(t),X(t),Y(t),Z(t)}and time t{0,100}.

Hi all,

I've been working on this problem for a while and only need to figure out why the function vector is returning zeros every time I call the procedure. The problem statement is in the file I attached to this.

If any of you could help me, I would be a happy guy.

 

File: 1.2.mw

Just as the title says: How do I avoid Maple Prompt to appear automatically?

When writing a code in Maple Prompt and hitting enter, a new [> is created and I can't delete it. Which is unfortunate, for most of the time I want to write some TEXT afterwards, not more code.

Sometimes it doesn't happen, but I can't find the logic in it. Any ideas? Much obliged!

For an example, see this post:

http://www.mapleprimes.com/questions/203735-Print-To-Pdf-Results-In-Blank-Space

2nd picture. I asked it as a BTW-question.

 

Thanks in advance

hi friends

i have a problem in maple with an error

  

Read'dsnumsort.map';
dsnumsort(Ns(0),[x,y,z]);
for i from 0 to 1000 do;
T:=i/25;
NsT:=Ns(T):
X[i]:=rhs(NsT[C1]); Vx[i]:=rhs(NsT[V1]);
Y[i]:=rhs(NsT[C2]); Vy[i]:=rhs(NsT[V2]);
Z[i]:=rhs(NsT[C3]); Vz[i]:=rhs(NsT[V3]);
KepVec[i]:=convert(crossprod([X[i],Y[i], Z[i]],[Vx[i],Vy[i], Vz[i]]),list);
KepAbs[i]:=norm(KepVec[i],2);
od:

pacecurve({[seq(X[i],Y[i],Z[i],i=0..1000)],[[-1/(2),0,0],[1/(2),0,0]]},labels=['x','y','z']);

but i see:

Error, invalid input: seq expects 2 arguments, but received 4

can you helpe me?

Thanks

Hello!
Please help solve this equation!

restart;

Dd:=4.9*10^(-10):
n:=4.5:
m:=3.5:
ny:=20:
B:=1.47*10^(-11):
H1:=70:
H2:=20:
H:=100:
b1:=40:
b2:=80:
delt_y:=H/ny:
A_max:=b1*H1+b2*H2+10*b1:
M:=1000:




b:=y->piecewise(`and`(y>=0,y<=H1),b1,`and`(y>H1,y<H1+H2),b2,`and`(y>=H1+H2,y<=H),b1):
y:=0: eq:=0:

for j from 1 to ny do
y:=(j-1)*delt_y+delt_y/2;

eq:=evalf(eq+b(y)*sign(y-y0)*abs((y-y0))^(1/n)*delt_y);
end do;


#eq=0, y0-?
sol1:=fsolve({eq=0},{y0});

Hi

I'm trying to use the command: DocumentTools[Retrieve](path,label), but i recieve the message "unable to retrieve label reference" no mather what i do. When i use the GUI: Insert - reference and get the value it works fine, but in that way I can't assign the value to a variable in the document.

here is an example:

The document i want to retrieve a value from:

restart;

a := 5;

#a gets the label (1)

 

The document i want to retrieve the value to:

restart;
with(DocumentTools);

path := FileTools:-JoinPath(["test.mw"], base = worksheetdir); = "C:\Users\Nicolai\CloudStation\Nicolai\Skole stof\Sem 7\test.mw"

 

Retrieve(path, "());
Retrieve(path, "(1)");
Retrieve(path, "a");

#None of the above works

 

Is it even possible to do such thing?

with best regards

Nicolai

with(Physics) :
t:=Intc(Dirac(k1+k2+k3)*phi(k1)*phi(k2)*phi(k3),k1,k2,k3) ;


# how to force
Fundiff(t,phi(-k));

# to find
3*Intc(Dirac(k2+k3-k)*phi(k2)*phi(k3),k2,k3) ;

 

First 1255 1256 1257 1258 1259 1260 1261 Last Page 1257 of 2224