Comment #4 on issue 51 by notmyadd...@teewars.org: RootOf for polynomial
equations
http://code.google.com/p/sympy/issues/detail?id=51
asmeurer ask me to add a comment on this.
I tried to get eigenvals
E = Symbol('E')
J = Symbol('J')
F = Symbol('F')
L = Symbol('L')
e = Matrix([[162900*E*J/(211*F*L**2), -546000*E*J/(211*F*L**3), 57600*E*J/
(211*F*L**2), -90000*E*J/(211*F*L**3), -15000*E*J/(211*F*L**2), -15000*E*J/
(211*F*L**3), 0, -29400*E*J/(211*F*L**2)], [-1552*E*J/(211*F*L), 39980*E*J/
(211*F*L**2), -1528*E*J/(211*F*L), -10800*E*J/(211*F*L**2),
-10240*E*J/(211*F*L),
-1800*E*J/(211*F*L**2), 0, -11968*E*J/(211*F*L)], [60480*E*J/(211*F*L**2),
-220200*E*J/(211*F*L**3), 167220*E*J/(211*F*L**2), -558000*E*J/(211*F*L**3),
33600*E*J/(211*F*L**2), -93000*E*J/(211*F*L**3), 0,
-55680*E*J/(211*F*L**2)],
[-376*E*J/(211*F*L), -29360*E*J/(211*F*L**2), -3024*E*J/(211*F*L),
31100*E*J/
(211*F*L**2), -12400*E*J/(211*F*L), -12400*E*J/(211*F*L**2), 0, -24304*E*J/
(211*F*L)], [10080*E*J/(211*F*L**2), -775200*E*J/(211*F*L**3), 59520*E*J/
(211*F*L**2), -93000*E*J/(211*F*L**3), 153300*E*J/(211*F*L**2), -543000*E*J/
(211*F*L**3), 0, -51480*E*J/(211*F*L**2)], [-208*E*J/(211*F*L), -42280*E*J/
(211*F*L**2), -2032*E*J/(211*F*L), -23200*E*J/(211*F*L**2),
-15120*E*J/(211*F*L),
31300*E*J/(211*F*L**2), 0, -35712*E*J/(211*F*L)], [-208*E*J/(211*F*L),
-42280*E*J/
(211*F*L**2), -2032*E*J/(211*F*L), -23200*E*J/(211*F*L**2),
-15120*E*J/(211*F*L),
-138200*E*J/(633*F*L**2), 500*E*J/(3*F*L**2), -115576*E*J/(633*F*L)], [0,
0, 0, 0,
0, -7000*E*J/(F*L**3), -5000*E*J/(F*L**3), 700*E*J/(F*L**2)]])
e.eigenvals()
but fails with sympy.polys.polynomial.PolynomialError: Can't decompose 1/F
What I get in maple is:
[ /13 2 (1/2)\ /13 2 (1/2)\
[100 |-- + - 31 | J E 100 |-- - - 31 | J E
[ \3 3 / \3 3 / 1 / / 6
[------------------------, ------------------------, ---- \20 RootOf\211 _Z
[ 2 2 2
[ L F L F L F
2 3 4
+ 685546875 - 856406250 _Z + 263728125 _Z - 30137500 _Z + 1452565 _Z
5 \ \ 1 / / 6
- 29722 _Z , index = 1/ J E/, ---- \20 RootOf\211 _Z + 685546875
2
L F
2 3 4 5
- 856406250 _Z + 263728125 _Z - 30137500 _Z + 1452565 _Z - 29722 _Z ,
\ \ 1 / / 6
index = 2/ J E/, ---- \20 RootOf\211 _Z + 685546875 - 856406250 _Z
2
L F
2 3 4 5 \ \
+ 263728125 _Z - 30137500 _Z + 1452565 _Z - 29722 _Z , index = 3/ J
E/,
1 / / 6 2
---- \20 RootOf\211 _Z + 685546875 - 856406250 _Z + 263728125 _Z
2
L F
3 4 5 \ \ 1
/ /
- 30137500 _Z + 1452565 _Z - 29722 _Z , index = 4/ J E/, ---- \20
RootOf\
2
L F
6 2 3
211 _Z + 685546875 - 856406250 _Z + 263728125 _Z - 30137500 _Z
4 5 \ \ 1 / / 6
+ 1452565 _Z - 29722 _Z , index = 5/ J E/, ---- \20 RootOf\211 _Z
2
L F
2 3 4
+ 685546875 - 856406250 _Z + 263728125 _Z - 30137500 _Z + 1452565 _Z
]
]
5 \ \]
- 29722 _Z , index = 6/ J E/]
]
]
--
You received this message because you are listed in the owner
or CC fields of this issue, or because you starred this issue.
You may adjust your issue notification preferences at:
http://code.google.com/hosting/settings
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to sympy-iss...@googlegroups.com.
To unsubscribe from this group, send email to
sympy-issues+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.