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]]