function [dydt] = odes2(t,y)

m=1;k=2;vb=0.1;mis=0.15;mik=0.1;g=9.81;Fn=m*g;Fs=mis*Fn;Fc=mik*Fn;vs=0.001;
sigma1=sqrt(10e5);zmax=1.839e-6;zba=1.287e-6;
sigma2=0.1;
sigma0=10e5;

if(vb-y(2))<0
    alfa=0
elseif y(3)<zba
    alfa=0
elseif zba<y(3) && y(3)<zmax
    alfa=0.5*(sin(pi*(y(3)-(zmax+zba)/2)/(zmax-zba))+1)
elseif zmax<y(3)
    alfa=1
end

dydt(1)=y(2);
dydt(3)=(vb-y(2))-(alfa*(sigma0*abs(vb-y(2)))/(Fc+(Fs-Fc)*exp(-((vb-y(2))/vs)^2)))*y(3);
dydt(2)=1/m*(sigma1*dydt(3)+sigma0*y(3)+sigma2*(vb-y(2))-k*y(1));

dydt=dydt(:);

end