MaplePrimes Questions

I am having problems with VectorCalculus[Norm] with argument passing.

In certain cases I get the error message:

     "Error, (in VectorCalculus:-GetCoordinates) the first argument, when present, must be a Vector"

In all cases, the first argument IS a vector.

Here is the case I need to succeed:

VectorCalculus[Norm](DeleteRow(Vector[column](5, [x, y, z, w, 1]), RowDimension(Vector[column](5, [x, y, z, w, 1]))))

or more completely:

here is the test Document with several cases of the problem explored:

Failure_on_Calculus_Verctor_Norm.mw

Background.  Calculus Vector Norm doesn't know about Homogenous systems so it inculdes an unwanted term in its norm.  Trimming off the last row of the vector fixes this problem.

(Note: as I am new, so if there is an existing cononical notation that will take care of that without further ado, then I could leave CalculusVector[Norm] to its own devices and move on.

Also, if I knew how to peek at the underside of such gotchas, I might have been able to sort it out myself.  Any leads in this direction would be taken to heart.  Thanks)

 I want to find AutomorphismGroup  of star S4 (ie, K1,3).

with(GraphTheory):
with(GroupTheory):
with(SpecialGraphs):
S4:=CompleteGraph(1,3);
G1 := AutomorphismGroup(S4); # out is GroupTheory:-PermutationGroup({module() ... end module,module() ... end module},degree=4)

 

 

the S_4 has six automorphisms: (1, 2, 3), (1, 3, 2), (2, 1, 3), (2, 3, 1), (3, 1, 2,), (3, 2, 1):

GraphAutomorphismStar

For the output of maple ,I cannot understand.

Maybe I don’t understand the group enough. May somebody explain in detail,  Thanks!

There is a whole chapter in programming guide on OO in Maple.

How much is OO actually used by Maple programmers? Is OO used much internally at Maplesoft itself in implementing internal code? I do not see questions about OO in Maple in the forum. May be it has not taken off? 

Any one knows of packages in Maple written in OO style, using the Object module? I think OO was added in Maple 16, which is 2012?

p.s. I myself like OO programming. Done alot of it in Java. I think it can be useful for large applications.

 

 

Hi I have product of function Q along with its derivative.  

Q'j(x)*Qi(x) for j=0,1 and i=0,...n.

how can i make loop  to collect all coefficeints for this function for j=0,1 and i=0,1,..n.

i hope some one will help me.

Thanks

Hello, i'm a newbie in maple. I need to simulate the movement of the conical pendulum.

The equations of motion  https://www.cyberforum.ru/cgi-bin/latex.cgi?\ddot{\theta }=cos\theta sin\theta {\dot{\varphi }}^2-\frac{gsin\theta}{l} and https://www.cyberforum.ru/cgi-bin/latex.cgi?\ddot{\varphi }=-\frac{2cos\theta \dot{\varphi }\dot{\theta }}{sin\theta }

I include a link to a description of the conical pendulum.http://olewitthansen.dk/Physics/Conical_pendulum.pdf

Thanks in advance! 

Hi,

How can I find the RMS of a function like x(t) in maple? the help is not clear

x(t) := -3.703703704*10^(-7)*(0.000111668023*cos(1000/33*sqrt(1122)*t) - 0.0001214712007*sin(1000/33*sqrt(1122)*t) - 0.0002325581396*sqrt(561)*sqrt(2)*(-0.0004467462845*sqrt(1122)*sin(1000/33*sqrt(1122)*t) + 0.0004467462845*sqrt(1122)*cos(1000/33*sqrt(1122)*t)))/((2.074226433*10^14*cos(1000/33*sqrt(1122)*t) + 2.074226433*10^14*sin(1000/33*sqrt(1122)*t))*(4.895037587*10^(-11) + 0.01685634229*(0.00001474262739*cos(1000/33*sqrt(1122)*t) + 0.00001474262739*sin(1000/33*sqrt(1122)*t))^2)^2);

Thanks,

Baharm31

 

Sorry, I have not used Maple for a long time...

I would like to know the derivative of Jy w.r.t. py1 in the following. I expect the answer to be 2 Wy6 py1-2 Wy6 py10 but I get 0. Am I not using the subscripts properly that Maple is not understanding me?


 

J__y := (1/2)[W__y1(a[c]^y-(diff(S[1](t), t, t))^y)^2+W__y2((diff(S[1](t), t))^y)^2+W__y3(S[1]^y)^2+W__y4(`#mscripts(mi("φ",fontstyle = "normal"),mi("c"),none(),none(),mo("."),none(),none())`-`#mscripts(mi("φ",fontstyle = "normal"),mi("S"),none(),none(),mo("."),none(),none())`)^2+W__y5(`φ__S`)^2+W__y6(p__y1-p__y10)^2+W__y7(p__y2-p__y20)^2+W__y8(p__y3-p__y30)^2+W__y9(p__y4-p__y40)^2]

(1/2)[W__y1(a[c]^y-(diff(diff(S[1](t), t), t))^y)^2+W__y2((diff(S[1](t), t))^y)^2+W__y3(S[1]^y)^2+W__y4(`#mscripts(mi("φ",fontstyle = "normal"),mi("c"),none(),none(),mo("."),none(),none())`-`#mscripts(mi("φ",fontstyle = "normal"),mi("S"),none(),none(),mo("."),none(),none())`)^2+W__y5(phi__S)^2+W__y6(p__y1-p__y10)^2+W__y7(p__y2-p__y20)^2+W__y8(p__y3-p__y30)^2+W__y9(p__y4-p__y40)^2]

(1)

diff(p[y*i](t), t) := K[yp*i](p[y*i]-p[y0*i])-K[py*i]*(diff(Jy, p[y*i]))

K[yp*i](p[y*i]-p[y0*i])-K[py*i]*(diff(Jy, p[y*i]))

(2)

Diff(J__y, p__y1)

Diff((1/2)[W__y1(a[c]^y-(diff(diff(S[1](t), t), t))^y)^2+W__y2((diff(S[1](t), t))^y)^2+W__y3(S[1]^y)^2+W__y4(`#mscripts(mi("φ",fontstyle = "normal"),mi("c"),none(),none(),mo("."),none(),none())`-`#mscripts(mi("φ",fontstyle = "normal"),mi("S"),none(),none(),mo("."),none(),none())`)^2+W__y5(phi__S)^2+W__y6(p__y1-p__y10)^2+W__y7(p__y2-p__y20)^2+W__y8(p__y3-p__y30)^2+W__y9(p__y4-p__y40)^2], p__y1)

(3)

NULL

``


 

Download cost_function_partial_derivative.mw

f1 := 0.001;
f2 := 0.002;
c1 := 0.002;
c2 := 0.005;
w1 := 0.1;
W := 0.14;
p1 := 0.65;
p2 := 0.28;
p := 1.167;
r := 0.004;
alpha := 0.2;
ze := 0.14;
mu := 0.05;
ga := 0.01;
cp := 9.3;
Rp := 100;
sigma2 := 0;
l[1] := 3*W^2*(-1+W^2)/(4*(-1+4*W^2));
l[2] := Rp^2*w1^2*mu*ga*cp/(2*(Rp^2*w1^2*cp^2+1));
l[3] := (Rp.w1.mu.ga)/(2*(Rp^2*w1^2*cp^2+1));
l[4] := W*(-1-7*W^2+8*W^4)/(8*(-2+8*W^2));
eq1 := (-x*t+l[2]*x+l[1]*x*y^2+(3/2)*alpha*ze^2*x)^2+((1/2)*c1*x+l[3]*x)^2-(1/4)*f1^2;
eq2 := (-l[4]*y^3+l[1]*x^2*y)^2+(1/4)*y^2*c2^2-f2^2/(4*W^2);

I need to solve these two equation for the variable x and y in terms of  t because I want to plot x with a range for t 

and the same plot y with a range for t and if it is possible to plot x vs y for the previuos same range for t

ode.mw

I found an assignment question concerning an ODE from 20 years ago. This was before I had discovered Maple.

So I put it through Maple and got a complex solution, But the lecturer's solution was non-imaginary.

For a) gamma <> Beta and b) gamma = beta.

Actually for b) I couldn't follow his solution

 

Consider

restart;
L := [1,2,3,4];
x[j];
seq(%, x[j] in L);
x[j];

What would you expect the output of the last line to be?

Answer:   Maple says 4.   Why?

 

Hi, 

In Maple it's possible to say that some quantity x has a value of 1 and that its unit is the meter (for instance).
But is it possible to say that x is a quantity whose fundamental unit is the 'length' (or any other fundamental unit, I insist on this point: I don't want to use derived units such as 'force', 'speed', ...)?

I think that the Units package doesn't allow to do this (or maybe I didn't figure out how).

Do you have any suggestions about that?

TIA


This is my notebook

 

restart

with(DifferentialGeometry)

DGsetup([t, rho, theta, phi, `#mover(mi("&theta;",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("&phi;",fontstyle = "normal"),mo("&uminus0;"))`, `#mover(mi("&theta;",fontstyle = "normal"),mo("&tilde;"))`, `#mover(mi("&phi;",fontstyle = "normal"),mo("&tilde;"))`], M)

M > 

``

M > 

evalDG(`&w`(dt, `d&rho;`))

_DG([["form", M, 1], [[[1], `d&rho;`]]])

(1)
M > 

evalDG(`&s`(dt, `d&rho;`))

_DG([["tensor", M, [["cov_bas"], []]], [[[1], `d&rho;`]]])

(2)
M > 

 

 

 

``

I don't understand why computing the wedge product of two differential forms gives me an output which prints the product of both. Why  don't it just return d rho ^ dt?

Hi Everyone, So i have a pretty ugly function which i am interesting in plotting. When i plot it seems that all values are real however when i do evalf with specific values of variable i am plotting over it shows that it is complex. Any thoughts on how to get it to stop plotting once it becomes complex? See attached maple file.

ComplexPlotQuestion.mw

Thanks. 

1. Plot the solid that lies between the cone z=(x^2+y^2)^(1/2) and the plane z = 4.

2. Plot the solid that lies above the cone z=(x^2+y^2)^(1/2)  and below the sphere z=x^2+y^2+z^2.

3. Plot the solid that is inside the hemisphere z=sqrt(25-x^2-y^2) , outside the cylinder x^2+y^2 = 4 and above the plane xy − plane.

4. Plot the solid bounded by the hemisphere z=sqrt(16-x^2-y^2)and the paraboloid z=x^2+y^2.

5. Plot the solid that lies within the sphere x^2+y^2_z^2=4 , above the xy − plane and below the cone z=sqrt(x^2+y^2).

Dear Community,

I have a simulation in MapleSim, and I have a probe in the system, for which measured data are also available. How can I load these measured data into MapleSim and compare them to the calculated results (probe data) visually? (Preferably in the "Analysis Window: Simulation Results" which pops up by pressing F6)

Tx for the kind help in advance,

best regards,

Andras

 

First 545 546 547 548 549 550 551 Last Page 547 of 2431