function dy=vzmet(t,y,par)
%izracuna desno stran sistema enacb za duseno nihanje
a=par(1);%koeficient vzmeti
b=par(2);%koeficient dusenja
x=y(1);%trenutni polozaj
v=y(2);%trenutna hitrost
dy=[v;-a*x-b*v];%Newtonova enacba