function xp = kolok1zad2( t,x )
b1 = 10;
b2 = 20;
m1 = 20;
m2 = 10;
k1 = 10;
k2 = 20;
k3 = 40;

y = (2*t).*(t>=0 & t<2) + (-2*t+8).*(t>=2 & t<4) + (2*t-8).*(t>=4 & t<6) + (-2*t+16).*(t>=6 & t<8) + (2*t-16).*(t>=8 & t<10) + (-2*t+24).*(t>=10 & t<12);
f = min(y,3);

xp = [  x(3)
        (1/m1)*(-2*b1*x(3)-k1*x(1)-k2*(x(1)-x(2)))
        x(4)
        (1/m2)*(k2*(x(1)-x(2))-f-k3*x(2)-b1*x(2))
     ];

end

% [t,x] = ode23(@kolok1zad2, [0 12], [0; 0; 0; 0]);
% plot(t,[x(:,1) x(:,4)])
% a = size(t);
% %distance shit
% put1 = sum(abs(x(2:a(1),1) - x(1:a(1)-1:1)));
% put2 = sum(abs(x(2:a(1),2) - x(1:a(1)-1:2)));
% visepresao = max(put1,put2);
% %w00t!