A weighted mean-square method of “cubication” for non-linear oscillators