Advanced Search

Carl Love

Carl Love

28110 Reputation

25 Badges

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

MaplePrimes Activity


These are replies submitted by Carl Love

repeated x-coordinates...

Carl Love 28110
October 23 2013
0 0

A problem is that some of your first coordinates are repeated. Can you get another digit of accuracy so that that is not the case?

IE not really working either...

Carl Love 28110
October 23 2013
0 2

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

Not really working in IE either...

Carl Love 28110
October 23 2013
0 0

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

Basis for the nullspace...

Carl Love 28110
October 23 2013
0 0

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

Basis for the nullspace...

Carl Love 28110
October 23 2013

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

Still thinking...

Carl Love 28110
October 23 2013

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

Can't view image...

Carl Love 28110
October 23 2013

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

It's the order of the derivatives...

Carl Love 28110
October 23 2013

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.

Here's my results...

Carl Love 28110
October 22 2013
0 0

@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

Here's my results...

Carl Love 28110
October 22 2013

@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

Not listed as thread safe...

Carl Love 28110
October 22 2013

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

It works for me...

Carl Love 28110
October 22 2013
0 0

@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

It works for me...

Carl Love 28110
October 22 2013

@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

Example?...

Carl Love 28110
October 22 2013
0 0

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

Example?...

Carl Love 28110
October 22 2013

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

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