^
q*2^
2q;
//n=nombres de termes et a fonction definissant a(n) //S_n(P_n) =seriealt(n,a) //S_n(P_n) approche sum((-1)^k*a(k),k,0,+infinity) //avec P_n=poly de chebyshev seriealt3(n,a):={ local k,d,c,p,q,S; q:=ceil(n/3); n:=3*q; d:=(-1)^q*2^(2*q); p:=0; c:=d-p; S:=a(0)*c; for (k:=1;k<q;k++) { c:=-p-c; S:=S+a(k)*c; } p:=1; c:=-c-p; S:=S+a(q)*c; for (k:=q+1;k<n;k++) { p:=p*(k-n-1)/(k-q); c:=-p-c; S:=S+a(k)*c; } return evalf(S/d); };