Carl Love

Carl Love

28130 Reputation

25 Badges

13 years, 309 days
Himself
Wayland, Massachusetts, United States
My name was formerly Carl Devore.

MaplePrimes Activity


These are replies submitted by Carl Love

It is not really working in IE either, but it's working better than in Firefox.

It's not really working in IE either, but it's working better than in Firefox.

@NoThik Only the first vector returned by IntegerLinearSolve satisfies the equation A.X=B. The remaining vectors are a basis for the nullspace of A.

@NoThik Only the first vector returned by IntegerLinearSolve satisfies the equation A.X=B. The remaining vectors are a basis for the nullspace of A.

@Alexander zhang I'm still thinking about it. I need to refresh my knowledge of surface integrals.

I can't view your image. Please upload a worksheet, not an image.

The diff_table is symmetric, so it won't let you control the order of the derivatives: F[x,y] is identical to F[y,x]. Look closely at the order that it places the derivatives:

lprint(F[x,y]);
diff(diff(f(x, y, w, z, a, b), x), y)
lprint(F[x,b]);
diff(diff(f(x, y, w, z, a, b), b), x)

In the first case the F[x] is an integral part of the expression, and hence it can be substituted. In the second case it is not an integral part of the expression, and hence it cannnot be substituted.

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


 

restart

kernelopts(version);

`Maple 16.00, X86 64 WINDOWS, Sep 16 2012, Build ID 773592`

(1)

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

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

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})

(2)

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

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

(3)

AG1_R0 := simplify(P_1_0.AG1_R1)

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

(4)

V_G1_R0 := map(diff, AG1_R0, t)

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})

(5)

``

``


Download gamma.mw

@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.

First 599 600 601 602 603 604 605 Last Page 601 of 710