Blob
1 scale = 502 define e(x) {3 auto a, b, c, d, e, g, w, y, t, r5 r = ibase6 ibase = A8 t = scale9 scale = t + .434*x + 111 w = 012 if(x<0) {13 x = -x14 w = 115 }16 y = 017 while(x>2) {18 x /= 219 y++20 }22 a = 123 b = 124 c = b25 d = 126 e = 127 for(a=1; 1; a++) {28 b *= x29 c = c*a+b30 d *= a31 g = c/d32 if(g == e) {33 g = g/134 while(y--) {35 g *= g36 }37 scale = t38 if(w==1) {39 ibase = r40 return 1/g41 }42 ibase = r43 return g/144 }45 e = g46 }47 }49 define l(x) {50 auto a, b, c, d, e, f, g, u, s, t, r, z52 r = ibase53 ibase = A54 if(x <= 0) {55 z = 1-10^scale56 ibase = r57 return z58 }59 t = scale61 f = 162 scale += scale(x) - length(x) + 163 s = scale64 while(x > 2) {65 s += (length(x)-scale(x))/2 + 166 if(s>0) {67 scale = s68 }69 x = sqrt(x)70 f *= 271 }72 while(x < .5) {73 s += (length(x)-scale(x))/2 + 174 if(s>0) {75 scale = s76 }77 x = sqrt(x)78 f *= 279 }81 scale = t + length(f) - scale(f) + 182 u = (x-1)/(x+1)84 scale += 1.1*length(t) - 1.1*scale(t)85 s = u*u86 b = 2*f87 c = b88 d = 189 e = 190 for(a=3; 1; a=a+2){91 b *= s92 c = c*a + d*b93 d *= a94 g = c/d95 if(g==e) {96 scale = t97 ibase = r98 return u*c/d99 }100 e = g101 }102 }104 define s(x) {105 auto a, b, c, s, t, y, p, n, i, r107 r = ibase108 ibase = A109 t = scale110 y = x/.7853111 s = t + length(y) - scale(y)112 if(s<t) {113 s = t114 }115 scale = s116 p = a(1)118 scale = 0119 if(x>=0) {120 n = (x/(2*p)+1)/2121 }122 if(x<0) {123 n = (x/(2*p)-1)/2124 }125 x -= 4*n*p126 if(n%2 != 0) {127 x = -x128 }130 scale = t + length(1.2*t) - scale(1.2*t)131 y = -x*x132 a = x133 b = 1134 s = x135 for(i=3; 1; i+=2) {136 a *= y137 b *= i*(i-1)138 c = a/b139 if(c==0){140 scale = t141 ibase = r142 return s/1143 }144 s += c145 }146 }148 define c(x) {149 auto t, r151 r = ibase152 ibase = A153 t = scale154 scale = scale+1155 x = s(x + 2*a(1))156 scale = t157 ibase = r158 return x/1159 }161 define a(x) {162 auto a, b, c, d, e, f, g, s, t, r, z164 r = ibase165 ibase = A166 if(x==0) {167 return 0168 }169 if(x==1) {170 z = .7853981633974483096156608458198757210492923498437764/1171 ibase = r172 if(scale<52) {173 return z174 }175 }176 t = scale177 f = 1178 while(x > .5) {179 scale++180 x = -(1 - sqrt(1.+x*x))/x181 f *= 2182 }183 while(x < -.5) {184 scale++185 x = -(1 - sqrt(1.+x*x))/x186 f *= 2187 }188 s = -x*x189 b = f190 c = f191 d = 1192 e = 1193 for(a=3; 1; a+=2) {194 b *= s195 c = c*a + d*b196 d *= a197 g = c/d198 if(g==e) {199 scale = t200 ibase = r201 return x*c/d202 }203 e = g204 }205 }207 define j(n,x) {208 auto a,b,c,d,e,g,i,s,k,t,r210 r = ibase211 ibase = A213 t = scale214 k = 1.36*x + 1.16*t - n215 k = length(k) - scale(k)216 if(k>0) {217 scale += k218 }220 s = -x*x/4221 if(n<0) {222 n = -n223 x = -x224 }225 a = 1226 c = 1227 for(i=1; i<=n; i++) {228 a *= x229 c *= 2*i230 }231 b = a232 d = 1233 e = 1234 for(i=1; 1; i++) {235 a *= s236 b = b*i*(n+i) + a237 c *= i*(n+i)238 g = b/c239 if(g==e) {240 scale = t241 ibase = r242 return g/1243 }244 e = g245 }246 }