quad precision
p1 = (-1.00000000000000000000000000000000e+00,0.00000000000000000000000000000000e+00,0.00000000000000000000000000000000e+00,1.00000000000000000000000000000000e+00)
p2 = (-1.00000000000000000000000000000000e+00,0.00000000000000000000000000000000e+00,0.00000000000000000000000000000000e+00,-1.00000000000000000000000000000000e+00)
p3 = (1.52392439451318091965667048882850e-01,-3.95662577808356314351356022854614e-02,9.33358581519928282825168263510354e-03,-1.46870184254580498554078492981181e-01)
p4 = (9.50786237105579294297122901823194e-01,-7.04407739589035804550791054540153e-01,1.59369349899834522262669912529847e-01,6.18389533699395726703130217077556e-01)
p5 = (8.96821323443102613737210049293936e-01,7.43973997369871435985926656825626e-01,-1.68702935715033805090921595164947e-01,-4.71519349444815228149051724096371e-01)
mur^2 = 2.00000000000000000000000000000000e+00
H10 = 6.58837658472508532019247118263321e+00
H20 = 2.72856525392978432211171702472908e+02