Binomial(2k,k).mws
>
with(plots):
>
plot([binomial(2*x, x)*2^(-2*x)*sqrt(Pi*x),1], x = 1..200, thickness = 2);
>
>