% funkcia ‘kyvadlo.m’ – kmitanie vychyleneho kyvadla s tlmením
function xdot = kyvadlo(t,x)
l=0.6; B=1e-2; g=9.81; m=0.2038;
xdot = [x(2); -B/m*x(2)-g/l*sin(x(1))];