    function z = fb(n,x,y)

    z=zeros(n,1);
    z(1)=10*y(1)*(1-y(2));
    z(2)=y(2)*(y(1)-1);