double c[]= {+4.00e+000, +1.50e+001 , +4.00e+001 , +8.50e+001 , +1.56e+002 , +2.59e+002, +4.00e+002 , +5.85e+002 , +8.20e+002, +1.11e+003, +1.50e+001, +8.50e+001, +2.59e+002 , +5.85e+002 , +1.11e+003 , +1.89e+003 , +2.96e+003 , +4.37e+003 , +6.18e+003 , +8.42e+003, +4.00e+001, +2.59e+002 , +8.20e+002 , +1.89e+003 , +3.62e+003 , +6.18e+003 , +9.72e+003 , +1.44e+004 , +2.04e+004 , +2.79e+004, +8.50e+001, +5.85e+002 , +1.89e+003 , +4.37e+003 , +8.42e+003 , +1.44e+004 , +2.28e+004 , +3.38e+004 , +4.80e+004 , +6.56e+004, +1.56e+002 , +1.11e+003 , +3.62e+003 , +8.42e+003 , +1.63e+004 , +2.79e+004 , +4.41e+004, +6.56e+004 , +9.32e+004, +1.28e+005, +2.59e+002, +1.89e+003, +6.18e+003 , +1.44e+004 ,+2.79e+004 ,+4.80e+004 , +7.59e+004, +1.13e+005 , +1.60e+005 , +2.20e+005, +4.00e+002 , +2.96e+003 , +9.72e+003 , +2.28e+004 , +4.41e+004 , +7.59e+004 , +1.20e+005 , +1.79e+005 , +2.54e+005 , +3.48e+005, +5.85e+002 , +4.37e+003 , +1.44e+004 , +3.38e+004 , +6.56e+004 , +1.13e+005 , +1.79e+005 , +2.66e+005 , +3.79e+005 , +5.18e+005, +8.20e+002 , +6.18e+003 , +2.04e+004 , +4.80e+004 , +9.32e+004 , +1.60e+005 , +2.54e+005 , +3.79e+005 , +5.38e+005 , +7.37e+005, +1.11e+003, +8.42e+003 , +2.79e+004 ,+6.56e+004 ,+1.28e+005 ,+2.20e+005 , +3.48e+005 , +5.18e+005 , +7.37e+005 , +1.01e+006}; double b[]= {+1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +2.00e+000 +4.00e+000 +8.00e+000 +1.00e+000 +3.00e+000 +9.00e+000 +2.70e+001 +1.00e+000 +4.00e+000 +1.60e+001 +6.40e+001 +1.00e+000 +5.00e+000 +2.50e+001 +1.25e+002 +1.00e+000 +6.00e+000 +3.60e+001 +2.16e+002 +1.00e+000 +7.00e+000 +4.90e+001 +3.43e+002 +1.00e+000 +8.00e+000 +6.40e+001 +5.12e+002 +1.00e+000 +9.00e+000 +8.10e+001 +7.29e+002 +1.00e+000 +1.00e+001 +1.00e+002 +1.00e+003]; a=[+1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +1.00e+000 +2.00e+000 +3.00e+000 +4.00e+000 +5.00e+000 +6.00e+000 +7.00e+000 +8.00e+000 +9.00e+000 +1.00e+001 +1.00e+000 +4.00e+000 +9.00e+000 +1.60e+001 +2.50e+001 +3.60e+001 +4.90e+001 +6.40e+001 +8.10e+001 +1.00e+002 +1.00e+000 +8.00e+000 +2.70e+001 +6.40e+001 +1.25e+002 +2.16e+002 +3.43e+002 +5.12e+002 +7.29e+002 +1.00e+003];