function y =hhorner(n)
% argument d'entrée:
% n : le degré du polynôme
% argument de sortie:
% y : la valeur de P_n(?^2)

alphacarre =  0.006694380022903415749574948586;

% calcul des coefficients du polynôme
a=zeros(n+1,1);
for i=1: n+1
    j= i-1;
    a(i) =factorial(2*j+1)/((2*j+3)*(2*j+5)*factorial(j)^2*2^(2*j-2));
end
y = a(n+1);
for i=n+1:-1:2
    
    y = a(i-1) + alphacarre*y;
end