Carl Love

Carl Love

28140 Reputation

25 Badges

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

MaplePrimes Activity


These are replies submitted by Carl Love

@Kitonum This model seems rather unlikely to me. The data are nondecreasing, but the model has a local maximum at x=1.83---just to the right of the end of the dataset. The marginal improvement in the fit that you get from the extra term (as measured by the first value returned by DataFit) seems too small to justify using the extra term.

@Adri van der Meer That should be

plot([X(t), Z(X(t)), t=0..10], numpoints= 200, color= blue, legend= "z(x)");

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?

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

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