tMax = 10; f = Sin[t];
nds = NDSolve[{x'''''[t] == f, x''''[0] == 0, x'''[0] == 0, 
    x''[0] == 0, x'[0] == 0, x[0] == 0}, x, {t, 0, tMax}, 
   MaxSteps -> Infinity, 
   Method -> {"ImplicitRungeKutta", "DifferenceOrder" -> 20}, 
   InterpolationOrder -> All];
crackle -> Evaluate[x'''''[tMax] /. nds][[1]]
velocity -> Evaluate[x'[tMax] /. nds][[1]]
position -> Evaluate[x[tMax] /. nds][[1]]