function [g] = gauf72b(x)
% Aufgabe 7.2 b)
%
% s.f.
global ngc tau

%tau=10^2;
ngc = ngc + 1;
g = [1+4*tau*(x(1,:).^2+x(2,:).^2-1).*x(1,:);
     4*tau*(x(1,:).^2+x(2,:).^2-1).*x(2,:)];

return

