matlab牛顿法

clear all clc format short n 0; x [ 1 5 0]"; e 1; for i 1:7 n n+1 g [400 x 1 ^3 x 1 x 2 +2 x 1 1 200 x 1 ^2+200 x 2 ]"; H [1200 x 1 ^2 400 x 2 +2 400 x 1 ; 400 x 1 200]; d inv H g x x+d e norm g 2 ; end