Question: what's the reason of the different result of 1d and 2d modes?

Dear all,

I come up with a problem. When I execute the following commands in 2d mode, error occurs:

restart; with(MTM); with(plots); r := .125; omega := .3; l := 1; m := 4.5; IG := (1/2)*m*r^2; A := Matrix(4, 4, {(1, 1) = 0, (1, 2) = 0, (1, 3) = 1, (1, 4) = 0, (2, 1) = 0, (2, 2) = 0, (2, 3) = 0, (2, 4) = 1, (3, 1) = 0, (3, 2) = r*omega^2*(IG+m*r*(r+l))/(l*IG), (3, 3) = -2*ldt/l, (3, 4) = 0, (4, 1) = 0, (4, 2) = -r*omega^2*(IG+m*(r+l)^2)/(l*IG), (4, 3) = 2*ldt/l, (4, 4) = 0}); l := eig(A); complexplot([l[1], l[2], l[3], l[4]], ldt = -.5 .. .5, style = point, labels = ["Re", "Im"], color = 'red');

However, when I execute in 1d mode, it turns to be successful.

I wonder what's the reason!

Thank you!

ex.mw

belief111

Please Wait...