function y = fn(n,x)
y = (1-(1-x).^n).^(1/n);
return
if n==1
    y = x;
else
    y = x + fn(n-1,x) - x.*fn(n-1,x);
end
return