phi1                = -pi;
dtphi1              = 0;    % pocetna brzina prvog tega
phi2                = -pi/2;
dtphi2              = 0;    % pocetna brzina drugog tega
g                   = 9.81; 
m1                  = 1; 
m2                  = 1; 
l1                  = 2; 
l2                  = 1;
vreme               = 10;
fps                 = 100;
movie               = true;

clear All; clf;
ivp = [phi1; dtphi1; phi2; dtphi2; g; m1; m2; l1; l2];
Dvostruko_klatno(ivp, vreme, fps, movie);

%[t y] = ode45(@Dvostruko_klatno_ODE,[0 vreme], ivp);
%figure(1)
%hold on
%plot(t,y(:,1), 'b-','LineWidth',2)
%plot(t,y(:,2), 'g-','LineWidth',2)
%plot(t,y(:,3), 'r-','LineWidth',2)
%plot(t,y(:,4), 'k-','LineWidth',2)
%xlabel('vreme');
%ylabel('ugao');
%hold off