def METHOD_NAME(x):
    small=1e-4
    if np.abs(x) < small:
        c = [-2.0/6.0,4.0/120.0,-6.0/5040.0,8.0/362880.0]
        x2=x*x
        return x*(c[0]+x2*(c[1]+x2*(c[2]+x2*c[3])))
    else:
        return np.cos(x)/x-np.sin(x)/(x*x)