Matlab - operaciones con variables simbolicas (comando sym)

clc

clear all

syms r m k x1 x2

M=[20 0;0 10] %matriz de masa

K=[2 -1;-1 1]*1e5 %matriz de rigidez

[xx,omega]=eig(inv(M)*K); %valores propios

det=det(M*r^2 + K)

fr_natural=sqrt(omega)

r2=-omega

A=[M*r^2 + K];

fila1=A(1,:)*[x1 x2]'

syms r m k x1 x2  real

>> syms_real1

M =

    20     0

     0    10

K =

      200000     -100000

     -100000      100000

det =

200*r^4 + 4000000*r^2 + 10000000000

fr_natural =

  130.6563         0

         0   54.1196

r2 =

   1.0e+04 *

   -1.7071         0

         0   -0.2929

fila1 =

x1*(20*r^2 + 200000) - 100000*x2