function DveMasi

%vnos parametrov naloge
m1=10; %masa1 v kilogramih
m2=1;
m3=5; %masa2 v kilogramih
L=2;%dol�ina palice (m)
F=50; %sila (N)

m=m1+m2+m3;
a=(m2*(L/2)+(m3*L))/m;
b=(L/2)-a;
koef=a*a*m1+b*b*m2+((b+(L/2))^2)*m3;

%PARAMETRI �ASOVNE INTEGRACIJE
tstart = 0; %za�etni �as
tfinal = 50; %kon�ni �as
y0 = [0,0,0,0,0,0]; %za�etne vrednosti vektorja neznank (stolpec vrednosti)


 % NUMERI�NO RE�EVANJE PROBLEMA
[t,y] = ode23t(@desne_strani,[tstart tfinal],y0,[],F,m,b,koef)
osi=[min(y(:,1))-L/2 max(y(:,1))+L/2 min(y(:,2))-L/2 max(y(:,2))+L/2]

% %GRAFI�NI PRIKAZ REZULTATOV
 figure(13);
 clf
 plot(t,y(:,1:3))
 legend("xs","ys","zasuk")
 axis equal

 figure(12);
 clf
 plot(t,y(:,1:3))
 legend("xs","ys","zasuk")
 axis equal

figure(14);
for i=1:size(t,1)
    
    palica(y(i,1),y(i,2),a,b,0.3,0.1,y(i,3))
    
    axis(osi);   
    film(i)=getframe;
end

% --------------------------------------------------------------------------
%OPIS DIFERENCIALNIH ENA�B
function dydt =desne_strani(t,y,F,m,b,koef)
dydt = [
    y(4)
    y(5)
    y(6)
    F*cos(t/10)/m
    F*sin(t/10)/m
    F*b*sin((t/10)-y(3))/koef
];