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) 1276498969817541.50000 ulp for j1(29.046828534916855) = -4.9308544767599701e-17) 1276498969817541.50000 ulp for j1(-29.046828534916855) = 4.9308544767599701e-17) 937918057966874.12000 ulp for j1(-22.760084380592772) = 7.0790211477589501e-17) 937918057966874.12000 ulp for j1(22.760084380592772) = -7.0790211477589501e-17) 523092309117987.81000 ulp for j1(60.469457845347492) = -2.6075935887631564e-17) 422992643178276.69000 ulp for j1(38.474766234771614) = -9.7736421293644016e-17) 390374805725733.69000 ulp for j1(41.61709421281445) = 7.5308076821311013e-17) 244279428281981.75000 ulp for j1(25.903672087618382) = -7.3698549660584186e-17) 244279428281981.75000 ulp for j1(-25.903672087618382) = 7.3698549660584186e-17) 113531257386682.34000 ulp for j1(22.760084380592769) = 5.2280311919105032e-16) 113531257386682.34000 ulp for j1(-22.760084380592769) = -5.2280311919105032e-16) 81201963312064.98400 ulp for j1(29.046828534916852) = 4.7639245074946587e-16) 81201963312064.98400 ulp for j1(-29.046828534916852) = -4.7639245074946587e-16) 78804159372664.96900 ulp for j1(63.611356698481231) = -1.7475614613417457e-16) 69196824179454.82000 ulp for j1(35.332307550083868) = -4.3427332352398035e-16) 68619564390854.68700 ulp for j1(44.759318997652819) = -2.8094955732782794e-16) 66375717593882.43700 ulp for j1(32.189679910974405) = 2.1931668859719546e-16) 51701009284245.95300 ulp for j1(38.474766234771621) = 8.1625531756198566e-16) 50008543473976.48400 ulp for j1(41.617094212814457) = -8.0328615276065077e-16) 45289337584552.50000 ulp for j1(47.901460887185451) = -3.9736883684577191e-16) j110241276498969817541.5000029.046828534916855-4.9308544767599701e-17