function dxdt = sedacka(t,x) F=800; dxdt = [x(2); -24*x(2)-80*x(1)+0.01*F]; return