% vlastní fekvence motoru
clc
I1=0.067017;                        %momenty setrvaènosti
I2=0.000132;
I3=0.016;
k1=140;                             %tuhosti spojek
k2=k1;


M=[ I1   0  0 
    0	I2  0                       %matice setrvaèností
    0    0  I3];

K=[     k1    -k1           0
       -k1    k1+k2       -k2       %matice tuhostí
        0      -k2         k2];

[V,L]= eig(K,M);

V ;                      %matice vektorù
L  ;                     %matice frekvence

lam1=L(1,1);             
lam2=L(2,2);             %vypsani jednotlivých lambda z matice frekvencí
lam3=L(3,3);


Om1=((sqrt(lam1))/(2*pi)),
Om2=((sqrt(lam2))/(2*pi)),   %uprava lambd na herze
Om3=((sqrt(lam3))/(2*pi)),






