此为一种拟牛顿法的一种方式,为matlab程序