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