%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 out=sistem(t,u,flag,beta)
switch flag
  case ''                                 % Return dy/dt = f(t,y).
    out = f(t,u,beta);
  
  otherwise
    error(['Unknown flag ''' flag '''.']);
  end
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))];
