function RHS=RHSfunction(t,w,k,c) 

RHS=zeros(6,1);


RHS(1)=w(2);

RHS(2)=-k*w(1)-c*w(2);


RHS(3)=w(4);

RHS(4)=-k*w(3)-c*w(4)-w(1);


RHS(5)=w(6);

RHS(6)=-k*w(5)-c*w(6)-w(2);