van=facebook.com/LearningCalculus$

product((1+1/k^2),k, 1,n)=Gamma((1 - %i) + n)*Gamma((1 + %i) + n)*Sinh(%pi)/Gamma(1 + n)^2;

y(x):= (Sinh(%pi)*Gamma(x-%i+1)*Gamma(x+%i+1))/Gamma(x+1)^2;

plot2d(y(x),[x,-%pi,%pi]);