matlab程序阻尼牛顿法