y=2;
for n=1:N
    y=y*((2*n)*(2*n))/((2*n-1)*(2*n+1));
end