Ask Your Question
0

Unable to make sense of Maxima expression

asked 2011-06-01 06:02:21 +0100

omerlevy gravatar image

updated 2011-06-01 19:49:34 +0100

Kelvin Li gravatar image

Hey everybody,

I'm trying to do some symbolic calculations. When trying to solve the following equation set for v_0, v_1, w_0, and w_1, I get "TypeError: Unable to make sense of Maxima expression".

The equations are added below:

v_0 == 1/6*(2*(s0_0 + s1_0 + s2_0)*t0_1 + 2*(s0_0 + s1_0 + s2_0)*t1_1 -
3*((t0_0 + t1_0)*t1_1 + t0_0*t0_1 + t0_1*t1_0)*w_0 - sqrt(-18*(t0_1^4 -
2*t0_1^3*t1_1 + 2*t0_1^2*t1_1^2 - 2*t0_1*t1_1^3 + t1_1^4)*w_1^2 -
9*(t0_0^2*t0_1^2 - 2*t0_0*t0_1^2*t1_0 + t0_1^2*t1_0^2 + (t0_0^2 -
2*t0_0*t1_0 + t1_0^2)*t1_1^2 + 2*(t0_0^2*t0_1 - 2*t0_0*t0_1*t1_0 +
t0_1*t1_0^2)*t1_1)*w_0^2 - 8*((s0_1 + s1_1)*s2_1 + (s0_0 + s1_0)*s2_0 -
s0_0^2 + s0_0*s1_0 - s0_1^2 + s0_1*s1_1 - s1_0^2 - s1_1^2 - s2_0^2 -
s2_1^2)*t0_1^2 - 16*((s0_1 + s1_1)*s2_1 + (s0_0 + s1_0)*s2_0 - s0_0^2 +
s0_0*s1_0 - s0_1^2 + s0_1*s1_1 - s1_0^2 - s1_1^2 - s2_0^2 -
s2_1^2)*t0_1*t1_1 - 8*((s0_1 + s1_1)*s2_1 + (s0_0 + s1_0)*s2_0 - s0_0^2
+ s0_0*s1_0 - s0_1^2 + s0_1*s1_1 - s1_0^2 - s1_1^2 - s2_0^2 -
s2_1^2)*t1_1^2))/(t0_1 + t1_1)

v_1 == -1/3*(3*(t0_1^2 + t1_1^2)*w_1 -
(s0_1 + s1_1 + s2_1)*t0_1 - (s0_1 + s1_1 + s2_1)*t1_1)/(t0_1 + t1_1),

w_0 == -1/6*(3*(t0_0*t0_1 + t0_1*t1_0)*t1_1*v_0 - ((s0_0 + s1_0 +
s2_0)*t0_0*t0_1 + (s0_0 + s1_0 + s2_0)*t0_1*t1_0)*t1_1 + sqrt((2*(s0_1 +
s1_1)*s2_1 + s0_1^2 + 2*s0_1*s1_1 + s1_1^2 +
s2_1^2)*t0_0*t0_1^3*t1_0*t1_1 + (2*(s0_1 + s1_1)*s2_1 + s0_1^2 +
2*s0_1*s1_1 + s1_1^2 + s2_1^2)*t0_0*t0_1*t1_0*t1_1^3 + 9*(t0_0^2*t0_1^2
- 2*t0_0*t0_1^2*t1_0 + t0_1^2*t1_0^2)*t1_1^2*v_0^2 - 6*((s0_0 + s1_0 +
s2_0)*t0_0^2*t0_1^2 - 2*(s0_0 + s1_0 + s2_0)*t0_0*t0_1^2*t1_0 + (s0_0 +
s1_0 + s2_0)*t0_1^2*t1_0^2)*t1_1^2*v_0 + ((2*(s0_0 + s1_0)*s2_0 + s0_0^2
+ 2*s0_0*s1_0 + s1_0^2 + s2_0^2)*t0_0^2*t0_1^2 + (2*(s0_0 + s1_0)*s2_0 +
s0_0^2 + 2*s0_0*s1_0 + s1_0^2 + s2_0^2)*t0_1^2*t1_0^2 + 2*(2*(s0_1 +
s1_1)*s2_1 + 2*(s0_0 + s1_0)*s2_0 - 5*s0_0^2 + 2*s0_0*s1_0 - 5*s0_1^2 +
2*s0_1*s1_1 - 5*s1_0^2 - 5*s1_1^2 - 5*s2_0^2 -
5*s2_1^2)*t0_0*t0_1^2*t1_0)*t1_1^2 + 9*(t0_0*t0_1^3*t1_0*t1_1 -
2*t0_0*t0_1^2*t1_0*t1_1^2 + t0_0*t0_1*t1_0*t1_1^3)*v_1^2 - 6*((s0_1 +
s1_1 + s2_1)*t0_0*t0_1^3*t1_0*t1_1 - 2*(s0_1 + s1_1 +
s2_1)*t0_0*t0_1^2*t1_0*t1_1^2 + (s0_1 + s1_1 +
s2_1)*t0_0*t0_1*t1_0*t1_1^3)*v_1))/(t0_0*t0_1*t1_0*t1_1)

w_1 == -1/6*(3*(t0_1 + t1_1)*v_1 - (s0_1 + s1_1 + s2_1)*t0_1 - (s0_1 + s1_1 +
s2_1)*t1_1)/(t0_1*t1_1)]

And the exact error is:

TypeError: unable to make sense of Maxima expression
'[If(and(-pi/2<parg(sqrt(s2_1^2+(-s1_1-s0_1)*s2_1+s2_0^2+(-s1_0-s0_0)\
*s2_0+s1_1^2-s0_1*s1_1+s1_0^2-s0_0*s1_0+s0_1^2+s0_0^2)*(sqrt(2)*t0_0*t0_ ...
(more)
edit retag flag offensive close merge delete

1 Answer

Sort by ยป oldest newest most voted
0

answered 2011-06-01 09:07:15 +0100

niles gravatar image

updated 2011-06-01 12:42:15 +0100

First, I noticed that there's a comma at the end of

v_1 == ...

and a square bracket at the end of

w_1 == ...

These are copy/paste typos, right?

Second, do the equations get solved successfully if you drop 1/2 of the terms or so? If so, that might point to memory overflow problems.

Third, could you also post exactly what solve command you're using?

edit flag offensive delete link more

Comments

Indeed copy/paste typos. Those terms were created dynamically, of course, and they're already at a minimum (yes, I know it doesn't look very much like it). I'm using the regular solve() command.

omerlevy gravatar imageomerlevy ( 2011-06-01 14:05:26 +0100 )edit

what I mean is enter the equations by hand (as I did, pasting from your post), but just delete most of the terms and see if the equations can stiil be solved. I deleted about 2/3 of the terms and the solve command worked without complaint; I didn't have time to test more.

niles gravatar imageniles ( 2011-06-01 17:23:59 +0100 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 2011-06-01 06:02:21 +0100

Seen: 735 times

Last updated: Jun 01 '11