vector algebra problem
- Subject: vector algebra problem
- From: Charles James Leonardo Quarra Cappiello
- Date: Wed, 06 Nov 2002 01:18:13 -0400
Hi,
Im having problems trying to derive a scalar function defined with dot
products:
(C1) LOAD(VECT);
(D1) /usr/local/share/maxima/5.9.0rc1/share/vector/vect.mac
(C2) LOAD(EIGEN);
Warning - you are redefining the MACSYMA function EIGENVALUES
Warning - you are redefining the MACSYMA function EIGENVECTORS
(D2) /usr/local/share/maxima/5.9.0rc1/share/matrix/eigen.mac
(C3) V1(t):= O1 + t*Ov1;
(D3) V1(t) := O1 + t Ov1
(C4) V0(t):= O0 + t*Ov0;
(D4) V0(t) := O0 + t Ov0
(C5) U1(t):= P1 + t*Pv1;
(D5) U1(t) := P1 + t Pv1
(C6) U0(t) := P0 + t*Pv0;
(D6) U0(t) := P0 + t Pv0
(C7) s0(t) := ((U1(t).U1(t))*(U0(t).V0(t) - U0(t).V1(t)) -
(U1(t).U0(t))*(U1(t).V0(t) - U1(t).V1(t)))/((U1(t) . U0(t)) - (U1(t) .
U1(t)) (U0(t) . U0(t)));
(D7) s0(t) := ((U1(t) . U1(t)) (U0(t) . V0(t) - U0(t) . V1(t))
- (U1(t) . U0(t)) (U1(t) . V0(t) - U1(t) . V1(t)))
/(U1(t) . U0(t) - U1(t) . U1(t)(U0(t) . U0(t)))
(C8) s1(t) := ((U0(t).U1(t))*(U0(t).V0(t) - U0(t).V1(t)) -
(U0(t).U0(t))*(U1(t).V0(t) - U1(t).V1(t)))/((U1(t) . U0(t)) - (U1(t) .
U1(t)) (U0(t) . U0(t)));
(D8) s1(t) := ((U0(t) . U1(t)) (U0(t) . V0(t) - U0(t) . V1(t))
- (U0(t) . U0(t)) (U1(t) . V0(t) - U1(t) . V1(t)))
/(U1(t) . U0(t) - U1(t) . U1(t)(U0(t) . U0(t)))
(C9) D2(t):= ((V1(t) - V0(t)) + s1(t)*U1(t) - s0(t)*U0(t) ).((V1(t) - V0(t))
+ s1(t)*U1(t) - s0(t)*U0(t) );
(D9) D2(t) := (V1(t) - V0(t) + s1(t) U1(t) + (- s0(t)) U0(t))
. (V1(t) - V0(t) + s1(t) U1(t) + (- s0(t)) U0(t))
(C10) diff(D2(t),t);
Improper name or value in functional position:
ERREXP1
#0: s1(t=t)
#1: D2(t=t)
-- an error. Quitting. To debug this try DEBUGMODE(TRUE);)
After that, Maxima crash in the next command (whatever it is) and i have to
restart. Im doing something wrong here?
Greetings,
Charlls
_________________________________________________________________
Help STOP SPAM with the new MSN 8 and get 2 months FREE*
http://join.msn.com/?page=features/junkmail