c = 1.5;

phi = @(x)( [ x(1) + x(2) + c*x(1)*(x(1)-1) ; x(2)+c*x(1)*(x(1)-1)]);


figure
Wu_g (  [1;0],11,  phi,'r-');
hold on;
Wu_g2( [1;0] ,5,  phi,'r-');
hold on;

Wu_g(  [1;0],11, @(x)( fsolve(@(s)(x - phi(s)) , [0;0])));
hold on;
Wu_g2( [1;0],4, @(x)( fsolve(@(s)(x - phi(s)) , [0;0])));