format(25)
p=1000000000100011
for i=1:p
calculpi=calculpi+1/(i*p)^2;
end
for j=(1:(p-1))
calculpi=calculpi+1/(p^2-j)^2;
end
calculpi=calculpi*(p^2-p+1)
calculpi/%pi^2*6
Petit programme sous scilab abouti
format(25)
b=997
A=zeros(b,b);
for i=1:b
for j=1:b
A(i,j)=modulo(i*j,b);
end
end
calculpi=0;
B=matrix(A,1,b^2);
indice=find(B==0);
for i=(1:b^(4))
C=(1:b^2)+(i-1)*b^2;
D=C(indice);
calculpi=calculpi+sum(D.^-2);
end
calculpi=calculpi*(b^2-b+1)
(%pi^2/6)/calculpi
J’avais découvert comment calculer une somme polynomiale en ne calculant pas tous les termes
Je voudrais aller plus loin pour faire des sommes avec des polynômes comportant des monomes à puissance négative
Bienvenue sur mon tout nouveau blog. Si vous voulez retrouvez mes anciens écrits allez sur amazon