(bc)
/* GNU bc code, computes a(0) through a(7) */
define trunc(x) { auto sc, t; sc=scale; scale=0; t=x/1; scale=sc; return(t) }
scale = 200; pi = 4 * a(1); r = e(pi * sqrt(163)); s = 640320;
c0 = 1 + trunc(r - s^3);
c1 = -1 - trunc(((s^3 + c0) - r) * s^3);
c2 = 1 + trunc((r - (s^3 + c0 + c1/s^3)) * s^6);
c3 = -1 - trunc(((s^3 + c0 + c1/s^3 + c2/s^6) - r) * s^9);
c4 = 1 + trunc((r - (s^3 + c0 + c1/s^3 + c2/s^6 + c3/s^9)) * s^12);
c5 = -1 - trunc(((s^3 + c0 + c1/s^3 + c2/s^6 + c3/s^9 + c4/s^12) - r) * s^15);
c6 = 1 + trunc((r - (s^3 + c0 + c1/s^3 + c2/s^6 + c3/s^9 + c4/s^12 + c5/s^15)) * s^18);
c7 = -1 - trunc(((s^3 + c0 + c1/s^3 + c2/s^6 + c3/s^9 + c4/s^12 + c5/s^15 + c6/s^18) - r) * s^21);
|