e=initell([0,0,1,-7,6]) level=5077 f = anell(e,5000) omega1=e[15] omega2=e[16] delta = e[12] smb(a,b,q,s1,s2,q1,q2,p1,p2,n,lambda,c1,c2,res) =\ print1("Computing symbol for ");print(a);\ b= (-1/a)% level;\ q=exp(-2*pi/level);\ s1=0;s2=0;\ q1 = exp(2*pi*i*a/level)*q;\ q2 = exp(2*pi*i*b/level)*q;\ p1=q1;p2=q2;\ for(n=1,5000,\ s1 = s1 + p1*f[n]/n;\ s2 = s2 + p2*f[n]/n;\ p1 = p1 * q1;\ p2 = p2 * q2);\ lambda=s1-s2 ;\ c2 = imag(lambda)/imag(omega2);\ print(c2);\ c2=round(c2);\ lambda=lambda-c2*omega2;\ c1=real(lambda)/real(omega1);\ print(c1);\ c1=round(c1);\ if(delta>0,res=c1+c2*i,\ res=c1 + c2*(1/2+1/2*i));\ res setup(n,u) =\ table = vector(level,dummy,x);\ put(1,0);\ for(n=2,level-1,\ u = (1-1/n)% level;\ if(table[u]<>x,put(n,-table[u]+table[n-1]),\ if(table[n]<>x,,put(n,smb(n)));\ put(u,-table[n]+table[n-1]))) put(n,s,u) =\ print1(" Entering the symbol ");print(n);\ print1(n); print1(" ");print1(table[n]); print1(" --> ");\ table[n]=s;\ print(table[n]);\ print1(level-n);print1(" ");print1(table[level-n]);print1(" --> ");\ table[level-n]= conj(s);\ print(table[level-n]);\ u = (1/n) % level;\ print1(u);print1(" ");\ print1(table[u]);print1(" --> ");\ table[u]= -conj(s);\ print(table[u]);\ print1(level-u);print1(" ");\ print1(table[level-u]);print1(" --> ");\ table[level-u]= -s;\ print(table[level-u]);