please explain the solution

Hi everybody,

Can't figure the solution of this code (should get 3 answres in form: w=number*k)

syms k;
syms w1;
m=240;I=10000;
Tuq = [0 1 -10;
       1 0 -5;
       1 0 5; 
      -1/sqrt(5) 2/sqrt(5) 21/sqrt(5);
      -1/sqrt(5) -2/sqrt(5) -21/sqrt(5)];
K = Tuq'*eye(5)*Tuq;
M = [ m 0 0 ; 0 m 0 ; 0 0 I ];
eq = det(K*k-w1*M);
s = solve(eq,w1);
w = real(sqrt(s));

Thank you very much.

Answer 1

[Ridwan Alam]

K = round(Tuq'*eye(5)*Tuq,2);
M = [ m 0 0 ; 0 m 0 ; 0 0 I ];
eq = det(K*k-w1*M);
s = solve(eq,w1);

will prevent sym() to consider the higher precision values.

Output:

vpa(s)
 
ans =
 
 0.009983*k
  0.03349*k
     0.01*k