%parametry

%polomìr øemenice motoru
R1=0.1;
%polomìr øemenice pilového kotouèe
R2=0.05;
%tuhost øemenu
k=1e4;
%tlumení øemenu
b=100*0;
%moment setrvaènosti øemenice motoru
J1=12;
%moment setrvaènosti øemenice pilového kotouèe
J2=6;

%vlastni frekvence netlumene soustavy
omega2 = sqrt ((J1*R2*k + J2*R1*k)/(J1*J2));
f = 1/(2*pi)*omega2

% symbolicky vypocet vlastnich cisel
syms k R1 R2 J1 J2

M = [J1 0; 0 J2];
K = k*[R1 -R1;-R2 R2];

lamda = eig(inv(M)*K)