function [x] = prod01(n)

if(n==1)
    x=1;
else
    x=1;
    for i=2:n
        x=x*(1+ (atan(n-1)/n));
    end        

end
end