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) 1289786618968880128.00000 ulp for j1(-7.01558666981561885) = -7.97717714238962506e-15) 1289786618968880128.00000 ulp for j1(7.01558666981561885) = 7.97717714238962506e-15) 135515044774005344.00000 ulp for j1(7.01558666981561796) = 6.44310461500700279e-15) 135515044774005344.00000 ulp for j1(-7.01558666981561796) = -6.44310461500700279e-15) 47610306648420016.00000 ulp for j1(7.01558666981562151) = 5.52266109857743322e-15) 47610306648420016.00000 ulp for j1(-7.01558666981562151) = -5.52266109857743322e-15) 21047863942359412.00000 ulp for j1(3.83170597020751247) = 1.97936879594712492e-16) 21047863942359412.00000 ulp for j1(-3.83170597020751247) = -1.97936879594712492e-16) 19327011183990332.00000 ulp for j1(3.83170597020751202) = 5.93810638784137426e-16) 19327011183990332.00000 ulp for j1(-3.83170597020751202) = -5.93810638784137426e-16) 8502797375605146.00000 ulp for j1(-3.83170597020751336) = -0) 8502797375605146.00000 ulp for j1(3.83170597020751336) = 0) 2716995199020033.50000 ulp for j1(-7.0155866698155922) = 3.68177406571828066e-15) 2716995199020033.50000 ulp for j1(7.0155866698155922) = -3.68177406571828066e-15) 452900336231367.68750 ulp for j1(-3.83170597020752401) = 4.35461135108368635e-15) 452900336231367.68750 ulp for j1(3.83170597020752401) = -4.35461135108368635e-15) 289044194357272.37500 ulp for j1(3.83170597020752046) = -3.16699007351540578e-15) 289044194357272.37500 ulp for j1(-3.83170597020752046) = 3.16699007351540578e-15) 162509535327770.90625 ulp for j1(7.015586669815006) = -1.79793300209234736e-13) 162509535327770.90625 ulp for j1(-7.015586669815006) = 1.79793300209234736e-13) j110241289786618968880128.00000-7.01558666981561885-7.97717714238962506e-15