def METHOD_NAME(r, h):
    d = np.sqrt(np.sum(r**2,axis=-1))
    q = d/h
    val = np.zeros_like(q)
    tag = np.where((q > 0) & (q <= 2))
    val[tag] = -35/(4*np.pi*h**3) * q[tag] * (1-q[tag]/2)**3 
    val[tag] /= d[tag]
    return val