3. MULTIPLE ZETA VALUES and LIN_DEP
>
> Zeta(alpha,beta)=Sum(1/n^alpha *Sum(1/k^beta,k=1..n-1),n=1..infinity);
> Zeta(-alpha,beta)=Sum((-1)^n/n^alpha* Sum(1/k^beta,k=1..n-1),n=1..infinity);
> Digits:=10;lin_dep([ez(3,1),Zeta(4)]);
> lin_dep([ez(2,1),ez(-2,1)]);
> lin_dep([ez(5,1),Zeta(6),Zeta(3)^2]);
> lin_dep([ez(3,1,3,1),Zeta(8)]);
> for k from 4 to 12 by 3 do ts:=time(): evalf(ez(3,1,3,1),k),time()-ts;od;
>