Normal execution of the
parabolic approximation


This is the normal output produced by the C implementation (the Maple output is identical). 
 1      6.000      6.000      0.000      0.000
 2     11.000     11.000      5.000      0.000
 3     18.000     18.000     -8.000      1.000
 4     27.000     27.000    -10.000      1.000
 5     38.000     38.000    -12.000      1.000
 6     51.000     51.000    -14.000      1.000
 7     66.000     66.000    -16.000      1.000
 8     83.000     83.000    -18.000      1.000
 9    102.000    102.000    -20.000      1.000
10    123.000    123.000    -22.000      1.000