def sinox3(x): '''(x-sin(x))/(x*x*x)''' if (abs(x) > 1e-2): return (x - math.sin(x)) / (x*x*x) else: return 1./6. - x*x / 120. + x*x*x*x / 5040.