function y = fun(z)
n = 2;
x = 0.2;
for i = 1:n
    x = x.*(1+z.*x);
end
y = x-1;