用一阶梯度法

				function y=dfm(x)
				    %y(1)=8*x(1)-4*x(2)+1;
				    %y(2)=6*x(2)-4*x(1);
				    %y=[y(1);y(2)];
				    
				    
				    y(1)=-400*x(1)*(x(2)-x(1)^2)+2*(x(1)-1);
				    y(2)=200*(x(2)-x(1)^2)+20.2*(x(2)-1)+19.8*(x(4)-1);
				    y(3)=-360*x(3)*(x(4)-x(3)^2)+2*(x(3)-1);
				    y(4)=180*(x(4)-x(3)^2)+20.2*(x(4)-1)+19.8*(x(2)-1);
				    y=[y(1);y(2);y(3);y(4)];