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