%SISTEM vrne desno stran sistema enacb za rot. vrv
% dudt=SISTEM(t,u,beta)
% t je vrednost parametra t=s/H
% u je vektor vrednosti funkcij u=[u1,u2,u3,u4]
% beta=g/(w^2*H)
function dudt=f(t,u,beta)
dudt=[-u(3)*cos(u(2))+beta*sin(u(2));
(u(3)*sin(u(2))+beta*cos(u(2)))/u(1);
cos(u(2));
sin(u(2));
-(u(7)*cos(u(2))-u(3)*u(6)*sin(u(2)))+beta*u(6)*cos(u(2));
1/u(1)*(u(7)*sin(u(2))+u(3)*u(6)*cos(u(2))-beta*u(6)*sin(u(2)))-u(5)/u(1)^2*(u(3)*sin(u(2))+beta*cos(u(2)));
-u(6)*sin(u(2));
u(6)*cos(u(2));
-(u(11)*cos(u(2))-u(3)*u(10)*sin(u(2)))+beta*u(10)*cos(u(2));
1/u(1)*(u(11)*sin(u(2))+u(3)*u(10)*cos(u(2))-beta*(u(10))*sin(u(2)))-u(9)/u(1)^2*(u(3)*sin(u(2))+beta*cos(u(2)));
-u(10)*sin(u(2));
u(10)*cos(u(2))];
