Carl Love

Carl Love

28055 Reputation

25 Badges

12 years, 361 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

@Bendesarts Here are my results. Download the attached worksheet and run it on your system.


restart:

kernelopts(version);

`Maple 16.02, X86 64 WINDOWS, Nov 18 2012, Build ID 788210`

constants:= ({constants} minus {gamma})[]:

`evalf/gamma`:= proc() end proc:

`evalf/constant/gamma`:= proc() end proc:

unprotect(gamma);

P_1_0:=Matrix([<cos(gamma(t)),-sin(gamma(t)),0>,<sin(gamma(t)), cos(gamma(t)),0 >,<0, 0, 1>]);

AG1_R1:=<(1/2)*r, 0, 0>;

AG1_R0:=simplify(P_1_0.AG1_R1);

V_G1_R0:=map(diff, AG1_R0, t);

 

P_1_0 := Matrix(3, 3, {(1, 1) = cos(gamma(t)), (1, 2) = sin(gamma(t)), (1, 3) = 0, (2, 1) = -sin(gamma(t)), (2, 2) = cos(gamma(t)), (2, 3) = 0, (3, 1) = 0, (3, 2) = 0, (3, 3) = 1})

AG1_R1 := Vector(3, {(1) = (1/2)*r, (2) = 0, (3) = 0})

AG1_R0 := Vector(3, {(1) = (1/2)*cos(gamma(t))*r, (2) = -(1/2)*sin(gamma(t))*r, (3) = 0})

V_G1_R0 := Vector(3, {(1) = -(1/2)*sin(gamma(t))*(diff(gamma(t), t))*r, (2) = -(1/2)*cos(gamma(t))*(diff(gamma(t), t))*r, (3) = 0})

 


Download gamma.mw

Physics is not listed at ?index,threadsafe , so it must be considered unsafe.

@Bendesarts It works for me in Maple 16. What error are you getting?

Note that MatrixVectorMultiply(A, B) can be shortened to simply A.B

@Bendesarts It works for me in Maple 16. What error are you getting?

Note that MatrixVectorMultiply(A, B) can be shortened to simply A.B

@Bendesarts Would you please show me an example of it not working? With the above commands, I am able to do diff(gamma(t), t) and simplify(gamma(t)) with no problem.

@Bendesarts Would you please show me an example of it not working? With the above commands, I am able to do diff(gamma(t), t) and simplify(gamma(t)) with no problem.

There is no attached file on your question.

@J4James Sorry, I didn't notice the attached file. I got it now.

What is it that the output doesn't match? Can you show a graph of the alternative?

It's an interesting problem. The simplified form that you got with solve andsubs contains sqrt(s[4]), even though the expressions to be simplified and the side relations are purely polynomial. The simplify-with-side-relations will never introduce sqrt like that. There may be some options though. I haven't been able to get it to work yet, but there may be a way where you change s[4] to s[4]^2.

Thanks for the addition. I had stretched out the computation for pedagogical reasons---to emphasize the linear algebra---thinking that the OP was taking a course in it.

Thanks for the addition. I had stretched out the computation for pedagogical reasons---to emphasize the linear algebra---thinking that the OP was taking a course in it.

No MaplePrimes 2D math mode, please, because it looks like this: int(1/(1+cot(x)), x);

No MaplePrimes 2D math mode, please, because it looks like this: int(1/(1+cot(x)), x);

There was no attached file on your post.

Generally only real symmetric or Hermitian matrices are classified by the "nature" of their definiteness. Those matrices only have real eigenvalues.

First 598 599 600 601 602 603 604 Last Page 600 of 709