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)];