result of j1 is never more precise than double 1 results were exact to double the precision 20 largest ulp errors (stored in a double) 500062037853668740.00000 ulp for j1(-7.0155866698156188) = -3.1101236429065784e-15) 500062037853668740.00000 ulp for j1(7.0155866698156188) = 3.1101236429065784e-15) 72766528853874528.00000 ulp for j1(7.015586669815618) = 3.3493639231301612e-15) 72766528853874528.00000 ulp for j1(-7.015586669815618) = -3.3493639231301612e-15) 18291900861198940.00000 ulp for j1(7.0155866698156215) = 2.6316430824594121e-15) 18291900861198940.00000 ulp for j1(-7.0155866698156215) = -2.6316430824594121e-15) 10098223402242856.00000 ulp for j1(3.8317059702075125) = -1.8596828691709446e-16) 10098223402242856.00000 ulp for j1(-3.8317059702075125) = 1.8596828691709446e-16) 2782959915422339.50000 ulp for j1(3.831705970207512) = 1.8596828691709444e-16) 2782959915422339.50000 ulp for j1(-3.831705970207512) = -1.8596828691709444e-16) 1259683190873853.20000 ulp for j1(-7.0155866698155922) = 5.9810070055895736e-15) 1259683190873853.20000 ulp for j1(7.0155866698155922) = -5.9810070055895736e-15) 959027438521293.50000 ulp for j1(3.8317059702075134) = -3.7193657383418892e-16) 959027438521293.50000 ulp for j1(-3.8317059702075134) = 3.7193657383418892e-16) 303062371422965.38000 ulp for j1(-3.8317059702075205) = 3.1614608775906103e-15) 303062371422965.38000 ulp for j1(3.8317059702075205) = -3.1614608775906103e-15) 156287044687378.31000 ulp for j1(-3.831705970207524) = 4.8351754598444657e-15) 156287044687378.31000 ulp for j1(3.831705970207524) = -4.8351754598444657e-15) 82072526676975.34400 ulp for j1(13.323691936314223) = -5.5770734366331129e-17) 82072526676975.34400 ulp for j1(-13.323691936314223) = 5.5770734366331129e-17) j11024500062037853668740.00000-7.0155866698156188-3.1101236429065784e-15