Li-Fan Chen says to YSITD
通常我們用的方法是 new_x <- x-n*v(x)