new force function. velver algorithm moved to velver.m