13 #include "lifev/core/array/VectorEpetra.hpp"    31             return vBar * (-0.9969);
    34             ERROR_MSG (
"This entrie is not allowed: ud_functions.hpp");
    42     Real vBar = 0.8 * linearPontdist (t, x, y, z, i) / 0.357434;
    46             return vBar * 0.113589;
    49             return -0.213165 * vBar;
    52             return vBar * (0.970391);
    55             ERROR_MSG (
"This entrie is not allowed: ud_functions.hpp");
    63     Real vBar = 0.2 * linearPontdist (t, x, y, z, i) / 0.357434;
    67             return vBar * -0.127337;
    70             return 0.3022305 * vBar;
    73             return vBar * (0.9446915);
    76             ERROR_MSG (
"This entrie is not allowed: ud_functions.hpp");
    85     Real vBar = linearPontdist (t, 0, 0, 0, 0) / 0.35;
   100     if (t <= 0.02 + 0.01)
   104     if (t <= 0.03 + 0.01)
   108     if (t <= 0.04 + 0.01)
   110         return   -9.1759e-06 * 1e6;
   112     if (t <= 0.05 + 0.01)
   114         return   -3.0930e-05 * 1e6;
   116     if (t <= 0.06 + 0.01)
   118         return   -6.2639e-05 * 1e6;
   120     if (t <= 0.07 + 0.01)
   122         return   -1.0212e-04 * 1e6;
   124     if (t <= 0.08 + 0.01)
   126         return   -1.4760e-04 * 1e6;
   128     if (t <= 0.09 + 0.01)
   130         return   -1.9726e-04 * 1e6;
   132     if (t <= 0.10 + 0.01)
   134         return   -2.4980e-04 * 1e6;
   136     if (t <= 0.11 + 0.01)
   138         return   -2.9526e-04 * 1e6;
   140     if (t <= 0.12 + 0.01)
   142         return   -3.2956e-04 * 1e6;
   144     if (t <= 0.13 + 0.01)
   146         return   -3.5469e-04 * 1e6;
   148     if (t <= 0.14 + 0.01)
   150         return   -3.7250e-04 * 1e6;
   152     if (t <= 0.15 + 0.01)
   154         return   -3.8429e-04 * 1e6;
   156     if (t <= 0.16 + 0.01)
   158         return   -3.9123e-04 * 1e6;
   160     if (t <= 0.17 + 0.01)
   162         return   -3.9431e-04 * 1e6;
   164     if (t <= 0.18 + 0.01)
   166         return   -3.9349e-04 * 1e6;
   168     if (t <= 0.19 + 0.01)
   170         return   -3.8858e-04 * 1e6;
   172     if (t <= 0.20 + 0.01)
   174         return   -3.7985e-04 * 1e6;
   176     if (t <= 0.21 + 0.01)
   178         return   -3.6756e-04 * 1e6;
   180     if (t <= 0.22 + 0.01)
   182         return   -3.5207e-04 * 1e6;
   184     if (t <= 0.23 + 0.01)
   186         return   -3.3408e-04 * 1e6;
   188     if (t <= 0.24 + 0.01)
   190         return   -3.1402e-04 * 1e6;
   192     if (t <= 0.25 + 0.01)
   194         return   -2.9288e-04 * 1e6;
   196     if (t <= 0.26 + 0.01)
   198         return   -2.7154e-04 * 1e6;
   200     if (t <= 0.27 + 0.01)
   202         return   -2.5054e-04 * 1e6;
   204     if (t <= 0.28 + 0.01)
   206         return   -2.2979e-04 * 1e6;
   208     if (t <= 0.29 + 0.01)
   210         return   -2.0904e-04 * 1e6;
   212     if (t <= 0.30 + 0.01)
   214         return   -1.8880e-04 * 1e6;
   216     if (t <= 0.31 + 0.01)
   218         return   -1.6899e-04 * 1e6;
   220     if (t <= 0.32 + 0.01)
   222         return   -1.4864e-04 * 1e6;
   224     if (t <= 0.33 + 0.01)
   226         return   -1.2730e-04 * 1e6;
   228     if (t <= 0.34 + 0.01)
   230         return   -1.0400e-04 * 1e6;
   232     if (t <= 0.35 + 0.01)
   234         return   -7.9755e-05 * 1e6;
   236     if (t <= 0.36 + 0.01)
   238         return   -5.8719e-05 * 1e6;
   240     if (t <= 0.37 + 0.01)
   242         return   -4.0345e-05 * 1e6;
   244     if (t <= 0.38 + 0.01)
   246         return   -2.4596e-05 * 1e6;
   248     if (t <= 0.39 + 0.01)
   250         return   -1.2259e-05 * 1e6;
   252     if (t <= 0.40 + 0.01)
   254         return   -3.8110e-06 * 1e6;
   256     if (t <= 0.41 + 0.01)
   260     if (t <= 0.42 + 0.01)
   264     if (t <= 0.43 + 0.01)
   268     if (t <= 0.44 + 0.01)
   272     if (t <= 0.45 + 0.01)
   276     if (t <= 0.46 + 0.01)
   280     if (t <= 0.47 + 0.01)
   284     if (t <= 0.48 + 0.01)
   288     if (t <= 0.49 + 0.01)
   292     if (t <= 0.50 + 0.01)
   296     if (t <= 0.51 + 0.01)
   300     if (t <= 0.52 + 0.01)
   304     if (t <= 0.53 + 0.01)
   308     if (t <= 0.54 + 0.01)
   312     if (t <= 0.55 + 0.01)
   316     if (t <= 0.56 + 0.01)
   320     if (t <= 0.57 + 0.01)
   324     if (t <= 0.58 + 0.01)
   328     if (t <= 0.59 + 0.01)
   332     if (t <= 0.60 + 0.01)
   336     if (t <= 0.61 + 0.01)
   340     if (t <= 0.62 + 0.01)
   344     if (t <= 0.63 + 0.01)
   348     if (t <= 0.64 + 0.01)
   352     if (t <= 0.65 + 0.01)
   356     if (t <= 0.66 + 0.01)
   360     if (t <= 0.67 + 0.01)
   364     if (t <= 0.68 + 0.01)
   368     if (t <= 0.69 + 0.01)
   372     if (t <= 0.70 + 0.01)
   376     if (t <= 0.71 + 0.01)
   380     if (t <= 0.72 + 0.01)
   384     if (t <= 0.73 + 0.01)
   388     if (t <= 0.74 + 0.01)
   392     if (t <= 0.75 + 0.01)
   396     if (t <= 0.76 + 0.01)
   400     if (t <= 0.77 + 0.01)
   404     if (t <= 0.78 + 0.01)
   408     if (t <= 0.79 + 0.01)
   421     if (t <= 0.00 + 0.01)
   423         return -1.e1 * 11017;
   425     if (t <= 0.01 + 0.01)
   427         return -1.e1 * 10954;
   429     if (t <= 0.02 + 0.01)
   431         return -1.e1 * 10893;
   433     if (t <= 0.03 + 0.01)
   435         return -1.e1 * 10832;
   437     if (t <= 0.04 + 0.01)
   439         return -1.e1 * 10771;
   441     if (t <= 0.05 + 0.01)
   443         return -1.e1 * 10712;
   445     if (t <= 0.06 + 0.01)
   447         return -1.e1 * 10653;
   449     if (t <= 0.07 + 0.01)
   451         return -1.e1 * 11113;
   453     if (t <= 0.08 + 0.01)
   455         return -1.e1 * 11544;
   457     if (t <= 0.09 + 0.01)
   459         return -1.e1 * 11869;
   461     if (t <= 0.10 + 0.01)
   463         return -1.e1 * 12146;
   465     if (t <= 0.11 + 0.01)
   467         return -1.e1 * 12394;
   469     if (t <= 0.12 + 0.01)
   471         return -1.e1 * 12635;
   473     if (t <= 0.13 + 0.01)
   475         return -1.e1 * 12889;
   477     if (t <= 0.14 + 0.01)
   479         return -1.e1 * 13151;
   481     if (t <= 0.15 + 0.01)
   483         return -1.e1 * 13398;
   485     if (t <= 0.16 + 0.01)
   487         return -1.e1 * 13620;
   489     if (t <= 0.17 + 0.01)
   491         return -1.e1 * 13833;
   493     if (t <= 0.18 + 0.01)
   495         return -1.e1 * 14035;
   497     if (t <= 0.19 + 0.01)
   499         return -1.e1 * 14229;
   501     if (t <= 0.20 + 0.01)
   503         return -1.e1 * 14436;
   505     if (t <= 0.21 + 0.01)
   507         return -1.e1 * 14613;
   509     if (t <= 0.22 + 0.01)
   511         return -1.e1 * 14753;
   513     if (t <= 0.23 + 0.01)
   515         return -1.e1 * 14878;
   517     if (t <= 0.24 + 0.01)
   519         return -1.e1 * 14974;
   521     if (t <= 0.25 + 0.01)
   523         return -1.e1 * 15032;
   525     if (t <= 0.26 + 0.01)
   527         return -1.e1 * 15047;
   529     if (t <= 0.27 + 0.01)
   531         return -1.e1 * 15025;
   533     if (t <= 0.28 + 0.01)
   535         return -1.e1 * 14975;
   537     if (t <= 0.29 + 0.01)
   539         return -1.e1 * 14899;
   541     if (t <= 0.30 + 0.01)
   543         return -1.e1 * 14822;
   545     if (t <= 0.31 + 0.01)
   547         return -1.e1 * 14721;
   549     if (t <= 0.32 + 0.01)
   551         return -1.e1 * 14594;
   553     if (t <= 0.33 + 0.01)
   555         return -1.e1 * 14496;
   557     if (t <= 0.34 + 0.01)
   559         return -1.e1 * 14375;
   561     if (t <= 0.35 + 0.01)
   563         return -1.e1 * 14198;
   565     if (t <= 0.36 + 0.01)
   567         return -1.e1 * 13990;
   569     if (t <= 0.37 + 0.01)
   571         return -1.e1 * 13726;
   573     if (t <= 0.38 + 0.01)
   575         return -1.e1 * 13397;
   577     if (t <= 0.39 + 0.01)
   579         return -1.e1 * 13167;
   581     if (t <= 0.40 + 0.01)
   583         return -1.e1 * 13132;
   585     if (t <= 0.41 + 0.01)
   587         return -1.e1 * 13315;
   589     if (t <= 0.42 + 0.01)
   591         return -1.e1 * 13271;
   593     if (t <= 0.43 + 0.01)
   595         return -1.e1 * 13157;
   597     if (t <= 0.44 + 0.01)
   599         return -1.e1 * 13028;
   601     if (t <= 0.45 + 0.01)
   603         return -1.e1 * 12975;
   605     if (t <= 0.46 + 0.01)
   607         return -1.e1 * 12933;
   609     if (t <= 0.47 + 0.01)
   611         return -1.e1 * 12891;
   613     if (t <= 0.48 + 0.01)
   615         return -1.e1 * 12836;
   617     if (t <= 0.49 + 0.01)
   619         return -1.e1 * 12768;
   621     if (t <= 0.50 + 0.01)
   623         return -1.e1 * 12700;
   625     if (t <= 0.51 + 0.01)
   627         return -1.e1 * 12641;
   629     if (t <= 0.52 + 0.01)
   631         return -1.e1 * 12592;
   633     if (t <= 0.53 + 0.01)
   635         return -1.e1 * 12548;
   637     if (t <= 0.54 + 0.01)
   639         return -1.e1 * 12504;
   641     if (t <= 0.55 + 0.01)
   643         return -1.e1 * 12456;
   645     if (t <= 0.56 + 0.01)
   647         return -1.e1 * 12405;
   649     if (t <= 0.57 + 0.01)
   651         return -1.e1 * 12353;
   653     if (t <= 0.58 + 0.01)
   655         return -1.e1 * 12300;
   657     if (t <= 0.59 + 0.01)
   659         return -1.e1 * 12244;
   661     if (t <= 0.60 + 0.01)
   663         return -1.e1 * 12184;
   665     if (t <= 0.61 + 0.01)
   667         return -1.e1 * 12122;
   669     if (t <= 0.62 + 0.01)
   671         return -1.e1 * 12058;
   673     if (t <= 0.63 + 0.01)
   675         return -1.e1 * 11995;
   677     if (t <= 0.64 + 0.01)
   679         return -1.e1 * 11933;
   681     if (t <= 0.65 + 0.01)
   683         return -1.e1 * 11871;
   685     if (t <= 0.66 + 0.01)
   687         return -1.e1 * 11810;
   689     if (t <= 0.67 + 0.01)
   691         return -1.e1 * 11747;
   693     if (t <= 0.68 + 0.01)
   695         return -1.e1 * 11684;
   697     if (t <= 0.69 + 0.01)
   699         return -1.e1 * 11620;
   701     if (t <= 0.70 + 0.01)
   703         return -1.e1 * 11556;
   705     if (t <= 0.71 + 0.01)
   707         return -1.e1 * 11492;
   709     if (t <= 0.72 + 0.01)
   711         return -1.e1 * 11428;
   713     if (t <= 0.73 + 0.01)
   715         return -1.e1 * 11365;
   717     if (t <= 0.74 + 0.01)
   719         return -1.e1 * 11302;
   721     if (t <= 0.75 + 0.01)
   723         return -1.e1 * 11240;
   725     if (t <= 0.76 + 0.01)
   727         return -1.e1 * 11179;
   729     if (t <= 0.77 + 0.01)
   731         return -1.e1 * 11120;
   733     if (t <= 0.78 + 0.01)
   735         return -1.e1 * 11062;
   737     if (t <= 0.79 + 0.01)
   739         return -1.e1 * 11006;
   751     if (t <= 0.00 + 0.01)
   753         return 3.350E-05 * 1e6;
   755     if (t <= 0.01 + 0.01)
   757         return 3.373E-05 * 1e6;
   759     if (t <= 0.02 + 0.01)
   761         return 3.402E-05 * 1e6;
   763     if (t <= 0.03 + 0.01)
   765         return 3.434E-05 * 1e6;
   767     if (t <= 0.04 + 0.01)
   769         return 3.466E-05 * 1e6;
   771     if (t <= 0.05 + 0.01)
   773         return 3.495E-05 * 1e6;
   775     if (t <= 0.06 + 0.01)
   777         return 3.519E-05 * 1e6;
   779     if (t <= 0.07 + 0.01)
   781         return 3.539E-05 * 1e6;
   783     if (t <= 0.08 + 0.01)
   785         return 3.564E-05 * 1e6;
   787     if (t <= 0.09 + 0.01)
   789         return 3.617E-05 * 1e6;
   791     if (t <= 0.10 + 0.01)
   793         return 3.773E-05 * 1e6;
   795     if (t <= 0.11 + 0.01)
   797         return 4.176E-05 * 1e6;
   799     if (t <= 0.12 + 0.01)
   801         return 5.037E-05 * 1e6;
   803     if (t <= 0.13 + 0.01)
   805         return 6.546E-05 * 1e6;
   807     if (t <= 0.14 + 0.01)
   809         return 8.701E-05 * 1e6;
   811     if (t <= 0.15 + 0.01)
   813         return 1.117E-04 * 1e6;
   815     if (t <= 0.16 + 0.01)
   817         return 1.345E-04 * 1e6;
   819     if (t <= 0.17 + 0.01)
   821         return 1.519E-04 * 1e6;
   823     if (t <= 0.18 + 0.01)
   825         return 1.642E-04 * 1e6;
   827     if (t <= 0.19 + 0.01)
   829         return 1.737E-04 * 1e6;
   831     if (t <= 0.20 + 0.01)
   833         return 1.821E-04 * 1e6;
   835     if (t <= 0.21 + 0.01)
   837         return 1.897E-04 * 1e6;
   839     if (t <= 0.22 + 0.01)
   841         return 1.958E-04 * 1e6;
   843     if (t <= 0.23 + 0.01)
   845         return 1.999E-04 * 1e6;
   847     if (t <= 0.24 + 0.01)
   849         return 2.019E-04 * 1e6;
   851     if (t <= 0.25 + 0.01)
   853         return 2.020E-04 * 1e6;
   855     if (t <= 0.26 + 0.01)
   857         return 2.004E-04 * 1e6;
   859     if (t <= 0.27 + 0.01)
   861         return 1.972E-04 * 1e6;
   863     if (t <= 0.28 + 0.01)
   865         return 1.926E-04 * 1e6;
   867     if (t <= 0.29 + 0.01)
   869         return 1.868E-04 * 1e6;
   871     if (t <= 0.30 + 0.01)
   873         return 1.798E-04 * 1e6;
   875     if (t <= 0.31 + 0.01)
   877         return 1.719E-04 * 1e6;
   879     if (t <= 0.32 + 0.01)
   881         return 1.632E-04 * 1e6;
   883     if (t <= 0.33 + 0.01)
   885         return 1.540E-04 * 1e6;
   887     if (t <= 0.34 + 0.01)
   889         return 1.446E-04 * 1e6;
   891     if (t <= 0.35 + 0.01)
   893         return 1.350E-04 * 1e6;
   895     if (t <= 0.36 + 0.01)
   897         return 1.254E-04 * 1e6;
   899     if (t <= 0.37 + 0.01)
   901         return 1.158E-04 * 1e6;
   903     if (t <= 0.38 + 0.01)
   905         return 1.062E-04 * 1e6;
   907     if (t <= 0.39 + 0.01)
   909         return 9.651E-05 * 1e6;
   911     if (t <= 0.40 + 0.01)
   913         return 8.634E-05 * 1e6;
   915     if (t <= 0.41 + 0.01)
   917         return 7.558E-05 * 1e6;
   919     if (t <= 0.42 + 0.01)
   921         return 6.447E-05 * 1e6;
   923     if (t <= 0.43 + 0.01)
   925         return 5.382E-05 * 1e6;
   927     if (t <= 0.44 + 0.01)
   929         return 4.484E-05 * 1e6;
   931     if (t <= 0.45 + 0.01)
   933         return 3.865E-05 * 1e6;
   935     if (t <= 0.46 + 0.01)
   937         return 3.556E-05 * 1e6;
   939     if (t <= 0.47 + 0.01)
   941         return 3.473E-05 * 1e6;
   943     if (t <= 0.48 + 0.01)
   945         return 3.457E-05 * 1e6;
   947     if (t <= 0.49 + 0.01)
   949         return 3.373E-05 * 1e6;
   951     if (t <= 0.50 + 0.01)
   953         return 3.191E-05 * 1e6;
   955     if (t <= 0.51 + 0.01)
   957         return 2.975E-05 * 1e6;
   959     if (t <= 0.52 + 0.01)
   961         return 2.809E-05 * 1e6;
   963     if (t <= 0.53 + 0.01)
   965         return 2.730E-05 * 1e6;
   967     if (t <= 0.54 + 0.01)
   969         return 2.718E-05 * 1e6;
   971     if (t <= 0.55 + 0.01)
   973         return 2.732E-05 * 1e6;
   975     if (t <= 0.56 + 0.01)
   977         return 2.744E-05 * 1e6;
   979     if (t <= 0.57 + 0.01)
   981         return 2.753E-05 * 1e6;
   983     if (t <= 0.58 + 0.01)
   985         return 2.772E-05 * 1e6;
   987     if (t <= 0.59 + 0.01)
   989         return 2.811E-05 * 1e6;
   991     if (t <= 0.60 + 0.01)
   993         return 2.866E-05 * 1e6;
   995     if (t <= 0.61 + 0.01)
   997         return 2.929E-05 * 1e6;
   999     if (t <= 0.62 + 0.01)
  1001         return 2.990E-05 * 1e6;
  1003     if (t <= 0.63 + 0.01)
  1005         return 3.044E-05 * 1e6;
  1007     if (t <= 0.64 + 0.01)
  1009         return 3.091E-05 * 1e6;
  1011     if (t <= 0.65 + 0.01)
  1013         return 3.132E-05 * 1e6;
  1015     if (t <= 0.66 + 0.01)
  1017         return 3.168E-05 * 1e6;
  1019     if (t <= 0.67 + 0.01)
  1021         return 3.199E-05 * 1e6;
  1023     if (t <= 0.68 + 0.01)
  1025         return 3.224E-05 * 1e6;
  1027     if (t <= 0.69 + 0.01)
  1029         return 3.244E-05 * 1e6;
  1031     if (t <= 0.70 + 0.01)
  1033         return 3.259E-05 * 1e6;
  1035     if (t <= 0.71 + 0.01)
  1037         return 3.270E-05 * 1e6;
  1039     if (t <= 0.72 + 0.01)
  1041         return 3.277E-05 * 1e6;
  1043     if (t <= 0.73 + 0.01)
  1045         return 3.281E-05 * 1e6;
  1047     if (t <= 0.74 + 0.01)
  1049         return 3.282E-05 * 1e6;
  1051     if (t <= 0.75 + 0.01)
  1053         return 3.283E-05 * 1e6;
  1055     if (t <= 0.76 + 0.01)
  1057         return 3.283E-05 * 1e6;
  1059     if (t <= 0.77 + 0.01)
  1061         return 3.286E-05 * 1e6;
  1063     if (t <= 0.78 + 0.01)
  1065         return 3.291E-05 * 1e6;
  1067     if (t <= 0.79 + 0.01)
  1069         return 3.300E-05 * 1e6;
  1075     if (t <= 0.00 + 0.01)
  1077         return  3.033E-06 * 1e6;
  1079     if (t <= 0.01 + 0.01)
  1081         return  3.041E-06 * 1e6;
  1083     if (t <= 0.02 + 0.01)
  1085         return  3.062E-06 * 1e6;
  1087     if (t <= 0.03 + 0.01)
  1089         return  3.094E-06 * 1e6;
  1091     if (t <= 0.04 + 0.01)
  1093         return  3.127E-06 * 1e6;
  1095     if (t <= 0.05 + 0.01)
  1097         return  3.150E-06 * 1e6;
  1099     if (t <= 0.06 + 0.01)
  1101         return  3.152E-06 * 1e6;
  1103     if (t <= 0.07 + 0.01)
  1105         return  3.141E-06 * 1e6;
  1107     if (t <= 0.08 + 0.01)
  1109         return  3.196E-06 * 1e6;
  1111     if (t <= 0.09 + 0.01)
  1113         return  3.574E-06 * 1e6;
  1115     if (t <= 0.10 + 0.01)
  1117         return  4.778E-06 * 1e6;
  1119     if (t <= 0.11 + 0.01)
  1121         return  7.387E-06 * 1e6;
  1123     if (t <= 0.12 + 0.01)
  1125         return  1.150E-05 * 1e6;
  1127     if (t <= 0.13 + 0.01)
  1129         return  1.609E-05 * 1e6;
  1131     if (t <= 0.14 + 0.01)
  1133         return  1.933E-05 * 1e6;
  1135     if (t <= 0.15 + 0.01)
  1137         return  2.007E-05 * 1e6;
  1139     if (t <= 0.16 + 0.01)
  1141         return  1.885E-05 * 1e6;
  1143     if (t <= 0.17 + 0.01)
  1145         return  1.706E-05 * 1e6;
  1147     if (t <= 0.18 + 0.01)
  1149         return  1.569E-05 * 1e6;
  1151     if (t <= 0.19 + 0.01)
  1153         return  1.481E-05 * 1e6;
  1155     if (t <= 0.20 + 0.01)
  1157         return  1.401E-05 * 1e6;
  1159     if (t <= 0.21 + 0.01)
  1161         return  1.294E-05 * 1e6;
  1163     if (t <= 0.22 + 0.01)
  1165         return  1.160E-05 * 1e6;
  1167     if (t <= 0.23 + 0.01)
  1169         return  1.018E-05 * 1e6;
  1171     if (t <= 0.24 + 0.01)
  1173         return  8.832E-06 * 1e6;
  1175     if (t <= 0.25 + 0.01)
  1177         return  7.609E-06 * 1e6;
  1179     if (t <= 0.26 + 0.01)
  1181         return  6.578E-06 * 1e6;
  1183     if (t <= 0.27 + 0.01)
  1185         return  5.843E-06 * 1e6;
  1187     if (t <= 0.28 + 0.01)
  1189         return  5.472E-06 * 1e6;
  1191     if (t <= 0.29 + 0.01)
  1193         return  5.412E-06 * 1e6;
  1195     if (t <= 0.30 + 0.01)
  1197         return  5.491E-06 * 1e6;
  1199     if (t <= 0.31 + 0.01)
  1201         return  5.527E-06 * 1e6;
  1203     if (t <= 0.32 + 0.01)
  1205         return  5.420E-06 * 1e6;
  1207     if (t <= 0.33 + 0.01)
  1209         return  5.169E-06 * 1e6;
  1211     if (t <= 0.34 + 0.01)
  1213         return  4.829E-06 * 1e6;
  1215     if (t <= 0.35 + 0.01)
  1217         return  4.465E-06 * 1e6;
  1219     if (t <= 0.36 + 0.01)
  1221         return  4.111E-06 * 1e6;
  1223     if (t <= 0.37 + 0.01)
  1225         return  3.750E-06 * 1e6;
  1227     if (t <= 0.38 + 0.01)
  1229         return  3.304E-06 * 1e6;
  1231     if (t <= 0.39 + 0.01)
  1233         return  2.668E-06 * 1e6;
  1235     if (t <= 0.40 + 0.01)
  1237         return  1.800E-06 * 1e6;
  1239     if (t <= 0.41 + 0.01)
  1241         return  8.269E-07 * 1e6;
  1243     if (t <= 0.42 + 0.01)
  1245         return  9.760E-08 * 1e6;
  1247     if (t <= 0.43 + 0.01)
  1249         return  7.311E-08 * 1e6;
  1251     if (t <= 0.44 + 0.01)
  1253         return  1.041E-06 * 1e6;
  1255     if (t <= 0.45 + 0.01)
  1257         return  2.783E-06 * 1e6;
  1259     if (t <= 0.46 + 0.01)
  1261         return  4.537E-06 * 1e6;
  1263     if (t <= 0.47 + 0.01)
  1265         return  5.488E-06 * 1e6;
  1267     if (t <= 0.48 + 0.01)
  1269         return  5.431E-06 * 1e6;
  1271     if (t <= 0.49 + 0.01)
  1273         return  4.863E-06 * 1e6;
  1275     if (t <= 0.50 + 0.01)
  1277         return  4.452E-06 * 1e6;
  1279     if (t <= 0.51 + 0.01)
  1281         return  4.499E-06 * 1e6;
  1283     if (t <= 0.52 + 0.01)
  1285         return  4.824E-06 * 1e6;
  1287     if (t <= 0.53 + 0.01)
  1289         return  5.059E-06 * 1e6;
  1291     if (t <= 0.54 + 0.01)
  1293         return  4.989E-06 * 1e6;
  1295     if (t <= 0.55 + 0.01)
  1297         return  4.671E-06 * 1e6;
  1299     if (t <= 0.56 + 0.01)
  1301         return  4.292E-06 * 1e6;
  1303     if (t <= 0.57 + 0.01)
  1305         return  3.981E-06 * 1e6;
  1307     if (t <= 0.58 + 0.01)
  1309         return  3.749E-06 * 1e6;
  1311     if (t <= 0.59 + 0.01)
  1313         return  3.553E-06 * 1e6;
  1315     if (t <= 0.60 + 0.01)
  1317         return  3.377E-06 * 1e6;
  1319     if (t <= 0.61 + 0.01)
  1321         return  3.255E-06 * 1e6;
  1323     if (t <= 0.62 + 0.01)
  1325         return  3.224E-06 * 1e6;
  1327     if (t <= 0.63 + 0.01)
  1329         return  3.281E-06 * 1e6;
  1331     if (t <= 0.64 + 0.01)
  1333         return  3.377E-06 * 1e6;
  1335     if (t <= 0.65 + 0.01)
  1337         return  3.452E-06 * 1e6;
  1339     if (t <= 0.66 + 0.01)
  1341         return  3.472E-06 * 1e6;
  1343     if (t <= 0.67 + 0.01)
  1345         return  3.441E-06 * 1e6;
  1347     if (t <= 0.68 + 0.01)
  1349         return  3.389E-06 * 1e6;
  1351     if (t <= 0.69 + 0.01)
  1353         return  3.343E-06 * 1e6;
  1355     if (t <= 0.70 + 0.01)
  1357         return  3.312E-06 * 1e6;
  1359     if (t <= 0.71 + 0.01)
  1361         return  3.289E-06 * 1e6;
  1363     if (t <= 0.72 + 0.01)
  1365         return  3.262E-06 * 1e6;
  1367     if (t <= 0.73 + 0.01)
  1369         return  3.223E-06 * 1e6;
  1371     if (t <= 0.74 + 0.01)
  1373         return  3.177E-06 * 1e6;
  1375     if (t <= 0.75 + 0.01)
  1377         return  3.132E-06 * 1e6;
  1379     if (t <= 0.76 + 0.01)
  1381         return  3.094E-06 * 1e6;
  1383     if (t <= 0.77 + 0.01)
  1385         return  3.065E-06 * 1e6;
  1387     if (t <= 0.78 + 0.01)
  1389         return  3.040E-06 * 1e6;
  1391     if (t <= 0.79 + 0.01)
  1393         return  3.016E-06 * 1e6;
  1399     if (t <= 0.00 + 0.01)
  1401         return  7.817E-07 * 1e6;
  1403     if (t <= 0.01 + 0.01)
  1405         return  7.879E-07 * 1e6;
  1407     if (t <= 0.02 + 0.01)
  1409         return  7.977E-07 * 1e6;
  1411     if (t <= 0.03 + 0.01)
  1413         return  8.077E-07 * 1e6;
  1415     if (t <= 0.04 + 0.01)
  1417         return  8.144E-07 * 1e6;
  1419     if (t <= 0.05 + 0.01)
  1421         return  8.164E-07 * 1e6;
  1423     if (t <= 0.06 + 0.01)
  1425         return  8.136E-07 * 1e6;
  1427     if (t <= 0.07 + 0.01)
  1429         return  8.129E-07 * 1e6;
  1431     if (t <= 0.08 + 0.01)
  1433         return  8.462E-07 * 1e6;
  1435     if (t <= 0.09 + 0.01)
  1437         return  9.858E-07 * 1e6;
  1439     if (t <= 0.10 + 0.01)
  1441         return  1.317E-06 * 1e6;
  1443     if (t <= 0.11 + 0.01)
  1445         return  1.854E-06 * 1e6;
  1447     if (t <= 0.12 + 0.01)
  1449         return  2.476E-06 * 1e6;
  1451     if (t <= 0.13 + 0.01)
  1453         return  3.005E-06 * 1e6;
  1455     if (t <= 0.14 + 0.01)
  1457         return  3.366E-06 * 1e6;
  1459     if (t <= 0.15 + 0.01)
  1461         return  3.626E-06 * 1e6;
  1463     if (t <= 0.16 + 0.01)
  1465         return  3.860E-06 * 1e6;
  1467     if (t <= 0.17 + 0.01)
  1469         return  4.044E-06 * 1e6;
  1471     if (t <= 0.18 + 0.01)
  1473         return  4.065E-06 * 1e6;
  1475     if (t <= 0.19 + 0.01)
  1477         return  3.824E-06 * 1e6;
  1479     if (t <= 0.20 + 0.01)
  1481         return  3.320E-06 * 1e6;
  1483     if (t <= 0.21 + 0.01)
  1485         return  2.659E-06 * 1e6;
  1487     if (t <= 0.22 + 0.01)
  1489         return  2.006E-06 * 1e6;
  1491     if (t <= 0.23 + 0.01)
  1493         return  1.503E-06 * 1e6;
  1495     if (t <= 0.24 + 0.01)
  1497         return  1.215E-06 * 1e6;
  1499     if (t <= 0.25 + 0.01)
  1501         return  1.117E-06 * 1e6;
  1503     if (t <= 0.26 + 0.01)
  1505         return  1.136E-06 * 1e6;
  1507     if (t <= 0.27 + 0.01)
  1509         return  1.190E-06 * 1e6;
  1511     if (t <= 0.28 + 0.01)
  1513         return  1.222E-06 * 1e6;
  1515     if (t <= 0.29 + 0.01)
  1517         return  1.219E-06 * 1e6;
  1519     if (t <= 0.30 + 0.01)
  1521         return  1.197E-06 * 1e6;
  1523     if (t <= 0.31 + 0.01)
  1525         return  1.179E-06 * 1e6;
  1527     if (t <= 0.32 + 0.01)
  1529         return  1.175E-06 * 1e6;
  1531     if (t <= 0.33 + 0.01)
  1533         return  1.176E-06 * 1e6;
  1535     if (t <= 0.34 + 0.01)
  1537         return  1.164E-06 * 1e6;
  1539     if (t <= 0.35 + 0.01)
  1541         return  1.129E-06 * 1e6;
  1543     if (t <= 0.36 + 0.01)
  1545         return  1.063E-06 * 1e6;
  1547     if (t <= 0.37 + 0.01)
  1549         return  9.647E-07 * 1e6;
  1551     if (t <= 0.38 + 0.01)
  1553         return  8.310E-07 * 1e6;
  1555     if (t <= 0.39 + 0.01)
  1557         return  6.635E-07 * 1e6;
  1559     if (t <= 0.40 + 0.01)
  1561         return  4.774E-07 * 1e6;
  1563     if (t <= 0.41 + 0.01)
  1565         return  3.116E-07 * 1e6;
  1567     if (t <= 0.42 + 0.01)
  1569         return  2.251E-07 * 1e6;
  1571     if (t <= 0.43 + 0.01)
  1573         return  2.627E-07 * 1e6;
  1575     if (t <= 0.44 + 0.01)
  1577         return  4.099E-07 * 1e6;
  1579     if (t <= 0.45 + 0.01)
  1581         return  5.913E-07 * 1e6;
  1583     if (t <= 0.46 + 0.01)
  1585         return  7.359E-07 * 1e6;
  1587     if (t <= 0.47 + 0.01)
  1589         return  8.403E-07 * 1e6;
  1591     if (t <= 0.48 + 0.01)
  1593         return  9.515E-07 * 1e6;
  1595     if (t <= 0.49 + 0.01)
  1597         return  1.097E-06 * 1e6;
  1599     if (t <= 0.50 + 0.01)
  1601         return  1.245E-06 * 1e6;
  1603     if (t <= 0.51 + 0.01)
  1605         return  1.331E-06 * 1e6;
  1607     if (t <= 0.52 + 0.01)
  1609         return  1.314E-06 * 1e6;
  1611     if (t <= 0.53 + 0.01)
  1613         return  1.204E-06 * 1e6;
  1615     if (t <= 0.54 + 0.01)
  1617         return  1.048E-06 * 1e6;
  1619     if (t <= 0.55 + 0.01)
  1621         return  9.004E-07 * 1e6;
  1623     if (t <= 0.56 + 0.01)
  1625         return  7.937E-07 * 1e6;
  1627     if (t <= 0.57 + 0.01)
  1629         return  7.358E-07 * 1e6;
  1631     if (t <= 0.58 + 0.01)
  1633         return  7.163E-07 * 1e6;
  1635     if (t <= 0.59 + 0.01)
  1637         return  7.186E-07 * 1e6;
  1639     if (t <= 0.60 + 0.01)
  1641         return  7.280E-07 * 1e6;
  1643     if (t <= 0.61 + 0.01)
  1645         return  7.359E-07 * 1e6;
  1647     if (t <= 0.62 + 0.01)
  1649         return  7.400E-07 * 1e6;
  1651     if (t <= 0.63 + 0.01)
  1653         return  7.423E-07 * 1e6;
  1655     if (t <= 0.64 + 0.01)
  1657         return  7.465E-07 * 1e6;
  1659     if (t <= 0.65 + 0.01)
  1661         return  7.556E-07 * 1e6;
  1663     if (t <= 0.66 + 0.01)
  1665         return  7.700E-07 * 1e6;
  1667     if (t <= 0.67 + 0.01)
  1669         return  7.871E-07 * 1e6;
  1671     if (t <= 0.68 + 0.01)
  1673         return  8.028E-07 * 1e6;
  1675     if (t <= 0.69 + 0.01)
  1677         return  8.129E-07 * 1e6;
  1679     if (t <= 0.70 + 0.01)
  1681         return  8.151E-07 * 1e6;
  1683     if (t <= 0.71 + 0.01)
  1685         return  8.101E-07 * 1e6;
  1687     if (t <= 0.72 + 0.01)
  1689         return  8.006E-07 * 1e6;
  1691     if (t <= 0.73 + 0.01)
  1693         return  7.905E-07 * 1e6;
  1695     if (t <= 0.74 + 0.01)
  1697         return  7.830E-07 * 1e6;
  1699     if (t <= 0.75 + 0.01)
  1701         return  7.792E-07 * 1e6;
  1703     if (t <= 0.76 + 0.01)
  1705         return  7.782E-07 * 1e6;
  1707     if (t <= 0.77 + 0.01)
  1709         return  7.782E-07 * 1e6;
  1711     if (t <= 0.78 + 0.01)
  1713         return  7.775E-07 * 1e6;
  1715     if (t <= 0.79 + 0.01)
  1717         return  7.754E-07 * 1e6;
  1723     if (t <= 0.00 + 0.01)
  1725         return  1.930E-06 * 1e6;
  1727     if (t <= 0.01 + 0.01)
  1729         return  1.710E-06 * 1e6;
  1731     if (t <= 0.02 + 0.01)
  1733         return  1.495E-06 * 1e6;
  1735     if (t <= 0.03 + 0.01)
  1737         return  1.289E-06 * 1e6;
  1739     if (t <= 0.04 + 0.01)
  1741         return  1.096E-06 * 1e6;
  1743     if (t <= 0.05 + 0.01)
  1745         return  9.184E-07 * 1e6;
  1747     if (t <= 0.06 + 0.01)
  1749         return  7.568E-07 * 1e6;
  1751     if (t <= 0.07 + 0.01)
  1753         return  6.240E-07 * 1e6;
  1755     if (t <= 0.08 + 0.01)
  1757         return  5.854E-07 * 1e6;
  1759     if (t <= 0.09 + 0.01)
  1761         return  8.113E-07 * 1e6;
  1763     if (t <= 0.10 + 0.01)
  1765         return  1.559E-06 * 1e6;
  1767     if (t <= 0.11 + 0.01)
  1769         return  3.032E-06 * 1e6;
  1771     if (t <= 0.12 + 0.01)
  1773         return  5.221E-06 * 1e6;
  1775     if (t <= 0.13 + 0.01)
  1777         return  7.903E-06 * 1e6;
  1779     if (t <= 0.14 + 0.01)
  1781         return  1.081E-05 * 1e6;
  1783     if (t <= 0.15 + 0.01)
  1785         return  1.376E-05 * 1e6;
  1787     if (t <= 0.16 + 0.01)
  1789         return  1.661E-05 * 1e6;
  1791     if (t <= 0.17 + 0.01)
  1793         return  1.919E-05 * 1e6;
  1795     if (t <= 0.18 + 0.01)
  1797         return  2.134E-05 * 1e6;
  1799     if (t <= 0.19 + 0.01)
  1801         return  2.292E-05 * 1e6;
  1803     if (t <= 0.20 + 0.01)
  1805         return  2.389E-05 * 1e6;
  1807     if (t <= 0.21 + 0.01)
  1809         return  2.433E-05 * 1e6;
  1811     if (t <= 0.22 + 0.01)
  1813         return  2.435E-05 * 1e6;
  1815     if (t <= 0.23 + 0.01)
  1817         return  2.406E-05 * 1e6;
  1819     if (t <= 0.24 + 0.01)
  1821         return  2.353E-05 * 1e6;
  1823     if (t <= 0.25 + 0.01)
  1825         return  2.280E-05 * 1e6;
  1827     if (t <= 0.26 + 0.01)
  1829         return  2.188E-05 * 1e6;
  1831     if (t <= 0.27 + 0.01)
  1833         return  2.078E-05 * 1e6;
  1835     if (t <= 0.28 + 0.01)
  1837         return  1.947E-05 * 1e6;
  1839     if (t <= 0.29 + 0.01)
  1841         return  1.797E-05 * 1e6;
  1843     if (t <= 0.30 + 0.01)
  1845         return  1.628E-05 * 1e6;
  1847     if (t <= 0.31 + 0.01)
  1849         return  1.445E-05 * 1e6;
  1851     if (t <= 0.32 + 0.01)
  1853         return  1.254E-05 * 1e6;
  1855     if (t <= 0.33 + 0.01)
  1857         return  1.060E-05 * 1e6;
  1859     if (t <= 0.34 + 0.01)
  1861         return  8.684E-06 * 1e6;
  1863     if (t <= 0.35 + 0.01)
  1865         return  6.838E-06 * 1e6;
  1867     if (t <= 0.36 + 0.01)
  1869         return  5.084E-06 * 1e6;
  1871     if (t <= 0.37 + 0.01)
  1873         return  3.412E-06 * 1e6;
  1875     if (t <= 0.38 + 0.01)
  1877         return  1.784E-06 * 1e6;
  1879     if (t <= 0.39 + 0.01)
  1881         return  1.534E-07 * 1e6;
  1883     if (t <= 0.40 + 0.01)
  1885         return  -1.494E-06 * 1e6;
  1887     if (t <= 0.41 + 0.01)
  1889         return  -3.093E-06 * 1e6;
  1891     if (t <= 0.42 + 0.01)
  1893         return  -4.495E-06 * 1e6;
  1895     if (t <= 0.43 + 0.01)
  1897         return  -5.521E-06 * 1e6;
  1899     if (t <= 0.44 + 0.01)
  1901         return  -6.086E-06 * 1e6;
  1903     if (t <= 0.45 + 0.01)
  1905         return  -6.252E-06 * 1e6;
  1907     if (t <= 0.46 + 0.01)
  1909         return  -6.160E-06 * 1e6;
  1911     if (t <= 0.47 + 0.01)
  1913         return  -5.908E-06 * 1e6;
  1915     if (t <= 0.48 + 0.01)
  1917         return  -5.508E-06 * 1e6;
  1919     if (t <= 0.49 + 0.01)
  1921         return  -4.944E-06 * 1e6;
  1923     if (t <= 0.50 + 0.01)
  1925         return  -4.234E-06 * 1e6;
  1927     if (t <= 0.51 + 0.01)
  1929         return  -3.441E-06 * 1e6;
  1931     if (t <= 0.52 + 0.01)
  1933         return  -2.639E-06 * 1e6;
  1935     if (t <= 0.53 + 0.01)
  1937         return  -1.873E-06 * 1e6;
  1939     if (t <= 0.54 + 0.01)
  1941         return  -1.160E-06 * 1e6;
  1943     if (t <= 0.55 + 0.01)
  1945         return  -5.019E-07 * 1e6;
  1947     if (t <= 0.56 + 0.01)
  1949         return  1.024E-07 * 1e6;
  1951     if (t <= 0.57 + 0.01)
  1953         return  6.543E-07 * 1e6;
  1955     if (t <= 0.58 + 0.01)
  1957         return  1.156E-06 * 1e6;
  1959     if (t <= 0.59 + 0.01)
  1961         return  1.614E-06 * 1e6;
  1963     if (t <= 0.60 + 0.01)
  1965         return  2.034E-06 * 1e6;
  1967     if (t <= 0.61 + 0.01)
  1969         return  2.419E-06 * 1e6;
  1971     if (t <= 0.62 + 0.01)
  1973         return  2.765E-06 * 1e6;
  1975     if (t <= 0.63 + 0.01)
  1977         return  3.065E-06 * 1e6;
  1979     if (t <= 0.64 + 0.01)
  1981         return  3.307E-06 * 1e6;
  1983     if (t <= 0.65 + 0.01)
  1985         return  3.486E-06 * 1e6;
  1987     if (t <= 0.66 + 0.01)
  1989         return  3.600E-06 * 1e6;
  1991     if (t <= 0.67 + 0.01)
  1993         return  3.654E-06 * 1e6;
  1995     if (t <= 0.68 + 0.01)
  1997         return  3.657E-06 * 1e6;
  1999     if (t <= 0.69 + 0.01)
  2001         return  3.621E-06 * 1e6;
  2003     if (t <= 0.70 + 0.01)
  2005         return  3.554E-06 * 1e6;
  2007     if (t <= 0.71 + 0.01)
  2009         return  3.465E-06 * 1e6;
  2011     if (t <= 0.72 + 0.01)
  2013         return  3.357E-06 * 1e6;
  2015     if (t <= 0.73 + 0.01)
  2017         return  3.233E-06 * 1e6;
  2019     if (t <= 0.74 + 0.01)
  2021         return  3.094E-06 * 1e6;
  2023     if (t <= 0.75 + 0.01)
  2025         return  2.941E-06 * 1e6;
  2027     if (t <= 0.76 + 0.01)
  2029         return  2.774E-06 * 1e6;
  2031     if (t <= 0.77 + 0.01)
  2033         return  2.591E-06 * 1e6;
  2035     if (t <= 0.78 + 0.01)
  2037         return  2.395E-06 * 1e6;
  2039     if (t <= 0.79 + 0.01)
  2041         return  2.185E-06 * 1e6;
  2047     if (t <= 0.00 + 0.01)
  2049         return  2.445E-06 * 1e6;
  2051     if (t <= 0.01 + 0.01)
  2053         return  2.470E-06 * 1e6;
  2055     if (t <= 0.02 + 0.01)
  2057         return  2.519E-06 * 1e6;
  2059     if (t <= 0.03 + 0.01)
  2061         return  2.578E-06 * 1e6;
  2063     if (t <= 0.04 + 0.01)
  2065         return  2.631E-06 * 1e6;
  2067     if (t <= 0.05 + 0.01)
  2069         return  2.665E-06 * 1e6;
  2071     if (t <= 0.06 + 0.01)
  2073         return  2.674E-06 * 1e6;
  2075     if (t <= 0.07 + 0.01)
  2077         return  2.674E-06 * 1e6;
  2079     if (t <= 0.08 + 0.01)
  2081         return  2.777E-06 * 1e6;
  2083     if (t <= 0.09 + 0.01)
  2085         return  3.260E-06 * 1e6;
  2087     if (t <= 0.10 + 0.01)
  2089         return  4.539E-06 * 1e6;
  2091     if (t <= 0.11 + 0.01)
  2093         return  6.956E-06 * 1e6;
  2095     if (t <= 0.12 + 0.01)
  2097         return  1.045E-05 * 1e6;
  2099     if (t <= 0.13 + 0.01)
  2101         return  1.438E-05 * 1e6;
  2103     if (t <= 0.14 + 0.01)
  2105         return  1.769E-05 * 1e6;
  2107     if (t <= 0.15 + 0.01)
  2109         return  1.958E-05 * 1e6;
  2111     if (t <= 0.16 + 0.01)
  2113         return  1.991E-05 * 1e6;
  2115     if (t <= 0.17 + 0.01)
  2117         return  1.910E-05 * 1e6;
  2119     if (t <= 0.18 + 0.01)
  2121         return  1.775E-05 * 1e6;
  2123     if (t <= 0.19 + 0.01)
  2125         return  1.628E-05 * 1e6;
  2127     if (t <= 0.20 + 0.01)
  2129         return  1.479E-05 * 1e6;
  2131     if (t <= 0.21 + 0.01)
  2133         return  1.326E-05 * 1e6;
  2135     if (t <= 0.22 + 0.01)
  2137         return  1.165E-05 * 1e6;
  2139     if (t <= 0.23 + 0.01)
  2141         return  1.001E-05 * 1e6;
  2143     if (t <= 0.24 + 0.01)
  2145         return  8.448E-06 * 1e6;
  2147     if (t <= 0.25 + 0.01)
  2149         return  7.065E-06 * 1e6;
  2151     if (t <= 0.26 + 0.01)
  2153         return  5.930E-06 * 1e6;
  2155     if (t <= 0.27 + 0.01)
  2157         return  5.089E-06 * 1e6;
  2159     if (t <= 0.28 + 0.01)
  2161         return  4.566E-06 * 1e6;
  2163     if (t <= 0.29 + 0.01)
  2165         return  4.350E-06 * 1e6;
  2167     if (t <= 0.30 + 0.01)
  2169         return  4.382E-06 * 1e6;
  2171     if (t <= 0.31 + 0.01)
  2173         return  4.552E-06 * 1e6;
  2175     if (t <= 0.32 + 0.01)
  2177         return  4.729E-06 * 1e6;
  2179     if (t <= 0.33 + 0.01)
  2181         return  4.798E-06 * 1e6;
  2183     if (t <= 0.34 + 0.01)
  2185         return  4.698E-06 * 1e6;
  2187     if (t <= 0.35 + 0.01)
  2189         return  4.428E-06 * 1e6;
  2191     if (t <= 0.36 + 0.01)
  2193         return  4.026E-06 * 1e6;
  2195     if (t <= 0.37 + 0.01)
  2197         return  3.513E-06 * 1e6;
  2199     if (t <= 0.38 + 0.01)
  2201         return  2.875E-06 * 1e6;
  2203     if (t <= 0.39 + 0.01)
  2205         return  2.076E-06 * 1e6;
  2207     if (t <= 0.40 + 0.01)
  2209         return  1.131E-06 * 1e6;
  2211     if (t <= 0.41 + 0.01)
  2213         return  1.861E-07 * 1e6;
  2215     if (t <= 0.42 + 0.01)
  2217         return  -4.577E-07 * 1e6;
  2219     if (t <= 0.43 + 0.01)
  2221         return  -4.657E-07 * 1e6;
  2223     if (t <= 0.44 + 0.01)
  2225         return  3.124E-07 * 1e6;
  2227     if (t <= 0.45 + 0.01)
  2229         return  1.684E-06 * 1e6;
  2231     if (t <= 0.46 + 0.01)
  2233         return  3.174E-06 * 1e6;
  2235     if (t <= 0.47 + 0.01)
  2237         return  4.306E-06 * 1e6;
  2239     if (t <= 0.48 + 0.01)
  2241         return  4.873E-06 * 1e6;
  2243     if (t <= 0.49 + 0.01)
  2245         return  4.980E-06 * 1e6;
  2247     if (t <= 0.50 + 0.01)
  2249         return  4.876E-06 * 1e6;
  2251     if (t <= 0.51 + 0.01)
  2253         return  4.757E-06 * 1e6;
  2255     if (t <= 0.52 + 0.01)
  2257         return  4.681E-06 * 1e6;
  2259     if (t <= 0.53 + 0.01)
  2261         return  4.603E-06 * 1e6;
  2263     if (t <= 0.54 + 0.01)
  2265         return  4.458E-06 * 1e6;
  2267     if (t <= 0.55 + 0.01)
  2269         return  4.224E-06 * 1e6;
  2271     if (t <= 0.56 + 0.01)
  2273         return  3.928E-06 * 1e6;
  2275     if (t <= 0.57 + 0.01)
  2277         return  3.618E-06 * 1e6;
  2279     if (t <= 0.58 + 0.01)
  2281         return  3.335E-06 * 1e6;
  2283     if (t <= 0.59 + 0.01)
  2285         return  3.098E-06 * 1e6;
  2287     if (t <= 0.60 + 0.01)
  2289         return  2.913E-06 * 1e6;
  2291     if (t <= 0.61 + 0.01)
  2293         return  2.778E-06 * 1e6;
  2295     if (t <= 0.62 + 0.01)
  2297         return  2.697E-06 * 1e6;
  2299     if (t <= 0.63 + 0.01)
  2301         return  2.670E-06 * 1e6;
  2303     if (t <= 0.64 + 0.01)
  2305         return  2.692E-06 * 1e6;
  2307     if (t <= 0.65 + 0.01)
  2309         return  2.746E-06 * 1e6;
  2311     if (t <= 0.66 + 0.01)
  2313         return  2.811E-06 * 1e6;
  2315     if (t <= 0.67 + 0.01)
  2317         return  2.868E-06 * 1e6;
  2319     if (t <= 0.68 + 0.01)
  2321         return  2.904E-06 * 1e6;
  2323     if (t <= 0.69 + 0.01)
  2325         return  2.915E-06 * 1e6;
  2327     if (t <= 0.70 + 0.01)
  2329         return  2.903E-06 * 1e6;
  2331     if (t <= 0.71 + 0.01)
  2333         return  2.872E-06 * 1e6;
  2335     if (t <= 0.72 + 0.01)
  2337         return  2.826E-06 * 1e6;
  2339     if (t <= 0.73 + 0.01)
  2341         return  2.769E-06 * 1e6;
  2343     if (t <= 0.74 + 0.01)
  2345         return  2.704E-06 * 1e6;
  2347     if (t <= 0.75 + 0.01)
  2349         return  2.636E-06 * 1e6;
  2351     if (t <= 0.76 + 0.01)
  2353         return  2.569E-06 * 1e6;
  2355     if (t <= 0.77 + 0.01)
  2357         return  2.511E-06 * 1e6;
  2359     if (t <= 0.78 + 0.01)
  2361         return  2.465E-06 * 1e6;
  2363     if (t <= 0.79 + 0.01)
  2365         return  2.433E-06 * 1e6;
  2371     if (t <= 0.00 + 0.01)
  2373         return  7.854E-07 * 1e6;
  2375     if (t <= 0.01 + 0.01)
  2377         return  7.900E-07 * 1e6;
  2379     if (t <= 0.02 + 0.01)
  2381         return  7.988E-07 * 1e6;
  2383     if (t <= 0.03 + 0.01)
  2385         return  8.077E-07 * 1e6;
  2387     if (t <= 0.04 + 0.01)
  2389         return  8.131E-07 * 1e6;
  2391     if (t <= 0.05 + 0.01)
  2393         return  8.141E-07 * 1e6;
  2395     if (t <= 0.06 + 0.01)
  2397         return  8.120E-07 * 1e6;
  2399     if (t <= 0.07 + 0.01)
  2401         return  8.085E-07 * 1e6;
  2403     if (t <= 0.08 + 0.01)
  2405         return  8.082E-07 * 1e6;
  2407     if (t <= 0.09 + 0.01)
  2409         return  8.255E-07 * 1e6;
  2411     if (t <= 0.10 + 0.01)
  2413         return  8.927E-07 * 1e6;
  2415     if (t <= 0.11 + 0.01)
  2417         return  1.055E-06 * 1e6;
  2419     if (t <= 0.12 + 0.01)
  2421         return  1.343E-06 * 1e6;
  2423     if (t <= 0.13 + 0.01)
  2425         return  1.728E-06 * 1e6;
  2427     if (t <= 0.14 + 0.01)
  2429         return  2.118E-06 * 1e6;
  2431     if (t <= 0.15 + 0.01)
  2433         return  2.419E-06 * 1e6;
  2435     if (t <= 0.16 + 0.01)
  2437         return  2.598E-06 * 1e6;
  2439     if (t <= 0.17 + 0.01)
  2441         return  2.673E-06 * 1e6;
  2443     if (t <= 0.18 + 0.01)
  2445         return  2.642E-06 * 1e6;
  2447     if (t <= 0.19 + 0.01)
  2449         return  2.465E-06 * 1e6;
  2451     if (t <= 0.20 + 0.01)
  2453         return  2.118E-06 * 1e6;
  2455     if (t <= 0.21 + 0.01)
  2457         return  1.660E-06 * 1e6;
  2459     if (t <= 0.22 + 0.01)
  2461         return  1.210E-06 * 1e6;
  2463     if (t <= 0.23 + 0.01)
  2465         return  8.723E-07 * 1e6;
  2467     if (t <= 0.24 + 0.01)
  2469         return  6.812E-07 * 1e6;
  2471     if (t <= 0.25 + 0.01)
  2473         return  6.155E-07 * 1e6;
  2475     if (t <= 0.26 + 0.01)
  2477         return  6.323E-07 * 1e6;
  2479     if (t <= 0.27 + 0.01)
  2481         return  6.927E-07 * 1e6;
  2483     if (t <= 0.28 + 0.01)
  2485         return  7.708E-07 * 1e6;
  2487     if (t <= 0.29 + 0.01)
  2489         return  8.538E-07 * 1e6;
  2491     if (t <= 0.30 + 0.01)
  2493         return  9.397E-07 * 1e6;
  2495     if (t <= 0.31 + 0.01)
  2497         return  1.030E-06 * 1e6;
  2499     if (t <= 0.32 + 0.01)
  2501         return  1.120E-06 * 1e6;
  2503     if (t <= 0.33 + 0.01)
  2505         return  1.203E-06 * 1e6;
  2507     if (t <= 0.34 + 0.01)
  2509         return  1.263E-06 * 1e6;
  2511     if (t <= 0.35 + 0.01)
  2513         return  1.292E-06 * 1e6;
  2515     if (t <= 0.36 + 0.01)
  2517         return  1.288E-06 * 1e6;
  2519     if (t <= 0.37 + 0.01)
  2521         return  1.258E-06 * 1e6;
  2523     if (t <= 0.38 + 0.01)
  2525         return  1.207E-06 * 1e6;
  2527     if (t <= 0.39 + 0.01)
  2529         return  1.142E-06 * 1e6;
  2531     if (t <= 0.40 + 0.01)
  2533         return  1.064E-06 * 1e6;
  2535     if (t <= 0.41 + 0.01)
  2537         return  9.816E-07 * 1e6;
  2539     if (t <= 0.42 + 0.01)
  2541         return  9.133E-07 * 1e6;
  2543     if (t <= 0.43 + 0.01)
  2545         return  8.870E-07 * 1e6;
  2547     if (t <= 0.44 + 0.01)
  2549         return  9.268E-07 * 1e6;
  2551     if (t <= 0.45 + 0.01)
  2553         return  1.035E-06 * 1e6;
  2555     if (t <= 0.46 + 0.01)
  2557         return  1.183E-06 * 1e6;
  2559     if (t <= 0.47 + 0.01)
  2561         return  1.329E-06 * 1e6;
  2563     if (t <= 0.48 + 0.01)
  2565         return  1.445E-06 * 1e6;
  2567     if (t <= 0.49 + 0.01)
  2569         return  1.529E-06 * 1e6;
  2571     if (t <= 0.50 + 0.01)
  2573         return  1.587E-06 * 1e6;
  2575     if (t <= 0.51 + 0.01)
  2577         return  1.612E-06 * 1e6;
  2579     if (t <= 0.52 + 0.01)
  2581         return  1.587E-06 * 1e6;
  2583     if (t <= 0.53 + 0.01)
  2585         return  1.501E-06 * 1e6;
  2587     if (t <= 0.54 + 0.01)
  2589         return  1.370E-06 * 1e6;
  2591     if (t <= 0.55 + 0.01)
  2593         return  1.230E-06 * 1e6;
  2595     if (t <= 0.56 + 0.01)
  2597         return  1.112E-06 * 1e6;
  2599     if (t <= 0.57 + 0.01)
  2601         return  1.029E-06 * 1e6;
  2603     if (t <= 0.58 + 0.01)
  2605         return  9.776E-07 * 1e6;
  2607     if (t <= 0.59 + 0.01)
  2609         return  9.460E-07 * 1e6;
  2611     if (t <= 0.60 + 0.01)
  2613         return  9.239E-07 * 1e6;
  2615     if (t <= 0.61 + 0.01)
  2617         return  9.061E-07 * 1e6;
  2619     if (t <= 0.62 + 0.01)
  2621         return  8.917E-07 * 1e6;
  2623     if (t <= 0.63 + 0.01)
  2625         return  8.821E-07 * 1e6;
  2627     if (t <= 0.64 + 0.01)
  2629         return  8.786E-07 * 1e6;
  2631     if (t <= 0.65 + 0.01)
  2633         return  8.809E-07 * 1e6;
  2635     if (t <= 0.66 + 0.01)
  2637         return  8.871E-07 * 1e6;
  2639     if (t <= 0.67 + 0.01)
  2641         return  8.940E-07 * 1e6;
  2643     if (t <= 0.68 + 0.01)
  2645         return  8.978E-07 * 1e6;
  2647     if (t <= 0.69 + 0.01)
  2649         return  8.959E-07 * 1e6;
  2651     if (t <= 0.70 + 0.01)
  2653         return  8.874E-07 * 1e6;
  2655     if (t <= 0.71 + 0.01)
  2657         return  8.734E-07 * 1e6;
  2659     if (t <= 0.72 + 0.01)
  2661         return  8.561E-07 * 1e6;
  2663     if (t <= 0.73 + 0.01)
  2665         return  8.383E-07 * 1e6;
  2667     if (t <= 0.74 + 0.01)
  2669         return  8.221E-07 * 1e6;
  2671     if (t <= 0.75 + 0.01)
  2673         return  8.088E-07 * 1e6;
  2675     if (t <= 0.76 + 0.01)
  2677         return  7.986E-07 * 1e6;
  2679     if (t <= 0.77 + 0.01)
  2681         return  7.909E-07 * 1e6;
  2683     if (t <= 0.78 + 0.01)
  2685         return  7.852E-07 * 1e6;
  2687     if (t <= 0.79 + 0.01)
  2689         return  7.807E-07 * 1e6;
  2695     if (t <= 0.00 + 0.01)
  2697         return  2.581E-06 * 1e6;
  2699     if (t <= 0.01 + 0.01)
  2701         return  2.499E-06 * 1e6;
  2703     if (t <= 0.02 + 0.01)
  2705         return  2.399E-06 * 1e6;
  2707     if (t <= 0.03 + 0.01)
  2709         return  2.281E-06 * 1e6;
  2711     if (t <= 0.04 + 0.01)
  2713         return  2.147E-06 * 1e6;
  2715     if (t <= 0.05 + 0.01)
  2717         return  2.002E-06 * 1e6;
  2719     if (t <= 0.06 + 0.01)
  2721         return  1.849E-06 * 1e6;
  2723     if (t <= 0.07 + 0.01)
  2725         return  1.695E-06 * 1e6;
  2727     if (t <= 0.08 + 0.01)
  2729         return  1.565E-06 * 1e6;
  2731     if (t <= 0.09 + 0.01)
  2733         return  1.537E-06 * 1e6;
  2735     if (t <= 0.10 + 0.01)
  2737         return  1.766E-06 * 1e6;
  2739     if (t <= 0.11 + 0.01)
  2741         return  2.460E-06 * 1e6;
  2743     if (t <= 0.12 + 0.01)
  2745         return  3.773E-06 * 1e6;
  2747     if (t <= 0.13 + 0.01)
  2749         return  5.709E-06 * 1e6;
  2751     if (t <= 0.14 + 0.01)
  2753         return  8.131E-06 * 1e6;
  2755     if (t <= 0.15 + 0.01)
  2757         return  1.087E-05 * 1e6;
  2759     if (t <= 0.16 + 0.01)
  2761         return  1.379E-05 * 1e6;
  2763     if (t <= 0.17 + 0.01)
  2765         return  1.675E-05 * 1e6;
  2767     if (t <= 0.18 + 0.01)
  2769         return  1.957E-05 * 1e6;
  2771     if (t <= 0.19 + 0.01)
  2773         return  2.205E-05 * 1e6;
  2775     if (t <= 0.20 + 0.01)
  2777         return  2.403E-05 * 1e6;
  2779     if (t <= 0.21 + 0.01)
  2781         return  2.538E-05 * 1e6;
  2783     if (t <= 0.22 + 0.01)
  2785         return  2.606E-05 * 1e6;
  2787     if (t <= 0.23 + 0.01)
  2789         return  2.611E-05 * 1e6;
  2791     if (t <= 0.24 + 0.01)
  2793         return  2.562E-05 * 1e6;
  2795     if (t <= 0.25 + 0.01)
  2797         return  2.474E-05 * 1e6;
  2799     if (t <= 0.26 + 0.01)
  2801         return  2.362E-05 * 1e6;
  2803     if (t <= 0.27 + 0.01)
  2805         return  2.238E-05 * 1e6;
  2807     if (t <= 0.28 + 0.01)
  2809         return  2.111E-05 * 1e6;
  2811     if (t <= 0.29 + 0.01)
  2813         return  1.981E-05 * 1e6;
  2815     if (t <= 0.30 + 0.01)
  2817         return  1.850E-05 * 1e6;
  2819     if (t <= 0.31 + 0.01)
  2821         return  1.715E-05 * 1e6;
  2823     if (t <= 0.32 + 0.01)
  2825         return  1.576E-05 * 1e6;
  2827     if (t <= 0.33 + 0.01)
  2829         return  1.432E-05 * 1e6;
  2831     if (t <= 0.34 + 0.01)
  2833         return  1.284E-05 * 1e6;
  2835     if (t <= 0.35 + 0.01)
  2837         return  1.132E-05 * 1e6;
  2839     if (t <= 0.36 + 0.01)
  2841         return  9.768E-06 * 1e6;
  2843     if (t <= 0.37 + 0.01)
  2845         return  8.180E-06 * 1e6;
  2847     if (t <= 0.38 + 0.01)
  2849         return  6.543E-06 * 1e6;
  2851     if (t <= 0.39 + 0.01)
  2853         return  4.831E-06 * 1e6;
  2855     if (t <= 0.40 + 0.01)
  2857         return  3.030E-06 * 1e6;
  2859     if (t <= 0.41 + 0.01)
  2861         return  1.163E-06 * 1e6;
  2863     if (t <= 0.42 + 0.01)
  2865         return  -6.817E-07 * 1e6;
  2867     if (t <= 0.43 + 0.01)
  2869         return  -2.362E-06 * 1e6;
  2871     if (t <= 0.44 + 0.01)
  2873         return  -3.738E-06 * 1e6;
  2875     if (t <= 0.45 + 0.01)
  2877         return  -4.742E-06 * 1e6;
  2879     if (t <= 0.46 + 0.01)
  2881         return  -5.400E-06 * 1e6;
  2883     if (t <= 0.47 + 0.01)
  2885         return  -5.785E-06 * 1e6;
  2887     if (t <= 0.48 + 0.01)
  2889         return  -5.956E-06 * 1e6;
  2891     if (t <= 0.49 + 0.01)
  2893         return  -5.932E-06 * 1e6;
  2895     if (t <= 0.50 + 0.01)
  2897         return  -5.723E-06 * 1e6;
  2899     if (t <= 0.51 + 0.01)
  2901         return  -5.358E-06 * 1e6;
  2903     if (t <= 0.52 + 0.01)
  2905         return  -4.889E-06 * 1e6;
  2907     if (t <= 0.53 + 0.01)
  2909         return  -4.370E-06 * 1e6;
  2911     if (t <= 0.54 + 0.01)
  2913         return  -3.846E-06 * 1e6;
  2915     if (t <= 0.55 + 0.01)
  2917         return  -3.341E-06 * 1e6;
  2919     if (t <= 0.56 + 0.01)
  2921         return  -2.866E-06 * 1e6;
  2923     if (t <= 0.57 + 0.01)
  2925         return  -2.422E-06 * 1e6;
  2927     if (t <= 0.58 + 0.01)
  2929         return  -2.004E-06 * 1e6;
  2931     if (t <= 0.59 + 0.01)
  2933         return  -1.601E-06 * 1e6;
  2935     if (t <= 0.60 + 0.01)
  2937         return  -1.206E-06 * 1e6;
  2939     if (t <= 0.61 + 0.01)
  2941         return  -8.123E-07 * 1e6;
  2943     if (t <= 0.62 + 0.01)
  2945         return  -4.195E-07 * 1e6;
  2947     if (t <= 0.63 + 0.01)
  2949         return  -3.198E-08 * 1e6;
  2951     if (t <= 0.64 + 0.01)
  2953         return  3.429E-07 * 1e6;
  2955     if (t <= 0.65 + 0.01)
  2957         return  6.970E-07 * 1e6;
  2959     if (t <= 0.66 + 0.01)
  2961         return  1.024E-06 * 1e6;
  2963     if (t <= 0.67 + 0.01)
  2965         return  1.320E-06 * 1e6;
  2967     if (t <= 0.68 + 0.01)
  2969         return  1.583E-06 * 1e6;
  2971     if (t <= 0.69 + 0.01)
  2973         return  1.813E-06 * 1e6;
  2975     if (t <= 0.70 + 0.01)
  2977         return  2.011E-06 * 1e6;
  2979     if (t <= 0.71 + 0.01)
  2981         return  2.180E-06 * 1e6;
  2983     if (t <= 0.72 + 0.01)
  2985         return  2.323E-06 * 1e6;
  2987     if (t <= 0.73 + 0.01)
  2989         return  2.442E-06 * 1e6;
  2991     if (t <= 0.74 + 0.01)
  2993         return  2.539E-06 * 1e6;
  2995     if (t <= 0.75 + 0.01)
  2997         return  2.613E-06 * 1e6;
  2999     if (t <= 0.76 + 0.01)
  3001         return  2.665E-06 * 1e6;
  3003     if (t <= 0.77 + 0.01)
  3005         return  2.691E-06 * 1e6;
  3007     if (t <= 0.78 + 0.01)
  3009         return  2.690E-06 * 1e6;
  3011     if (t <= 0.79 + 0.01)
  3013         return  2.661E-06 * 1e6;
  3052     if (z < 5 && z >= 0)
  3059         return  60 * 1e3 + 30 * (z - 5) * 1e3;    
  3071     if ( z < 0 && z >= -3 )
  3097             ERROR_MSG (
"This entrie is not allowed: ud_functions.hpp");
  3117             ERROR_MSG (
"This entrie is not allowed: ud_functions.hpp");
  3128         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3131     Real ti = floor (t * 100.0) / 100.0;
  3132     Real tii = ti + 0.01;
  3141         t = ( ( (
int) floor (t * 1000.0) ) % 800) / 1000.0;
  3144     Real ti = floor (t * 100.0) / 100.0;
  3145     Real tii = ti + 0.01;
  3154         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3157     Real ti = floor (t * 100.0) / 100.0;
  3158     Real tii = ti + 0.01;
  3168         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3171     Real ti = floor (t * 100.0) / 100.0;
  3172     Real tii = ti + 0.01;
  3182         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3185     Real ti = floor (t * 100) / 100.0;
  3186     Real tii = ti + 0.01;
  3196         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3199     Real ti = floor (t * 100) / 100.0;
  3200     Real tii = ti + 0.01;
  3209         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3212     Real ti = floor (t * 100) / 100.0;
  3213     Real tii = ti + 0.01;
  3222         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3225     Real ti = floor (t * 100) / 100.0;
  3226     Real tii = ti + 0.01;
  3234         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3237     Real ti = floor (t * 100) / 100.0;
  3238     Real tii = ti + 0.01;
  3247         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3250     Real ti = floor (t * 100) / 100.0;
  3251     Real tii = ti + 0.01;
  3260         t = ( ( (
int) floor (t * 1000) ) % 800) / 1000.0;
  3263     Real ti = floor (t * 100) / 100.0;
  3264     Real tii = ti + 0.01;
  3270     if (t <= 0.00 + 0.01)
  3274     if (t <= 0.01 + 0.01)
  3278     if (t <= 0.02 + 0.01)
  3282     if (t <= 0.03 + 0.01)
  3286     if (t <= 0.04 + 0.01)
  3290     if (t <= 0.05 + 0.01)
  3294     if (t <= 0.06 + 0.01)
  3298     if (t <= 0.07 + 0.01)
  3302     if (t <= 0.08 + 0.01)
  3306     if (t <= 0.09 + 0.01)
  3310     if (t <= 0.10 + 0.01)
  3314     if (t <= 0.11 + 0.01)
  3318     if (t <= 0.12 + 0.01)
  3322     if (t <= 0.13 + 0.01)
  3324         return  65.55999999999999;
  3326     if (t <= 0.14 + 0.01)
  3328         return  87.10999999999999;
  3330     if (t <= 0.15 + 0.01)
  3334     if (t <= 0.16 + 0.01)
  3338     if (t <= 0.17 + 0.01)
  3342     if (t <= 0.18 + 0.01)
  3346     if (t <= 0.19 + 0.01)
  3350     if (t <= 0.20 + 0.01)
  3354     if (t <= 0.21 + 0.01)
  3358     if (t <= 0.22 + 0.01)
  3362     if (t <= 0.23 + 0.01)
  3366     if (t <= 0.24 + 0.01)
  3370     if (t <= 0.25 + 0.01)
  3374     if (t <= 0.26 + 0.01)
  3378     if (t <= 0.27 + 0.01)
  3382     if (t <= 0.28 + 0.01)
  3386     if (t <= 0.29 + 0.01)
  3390     if (t <= 0.30 + 0.01)
  3394     if (t <= 0.31 + 0.01)
  3398     if (t <= 0.32 + 0.01)
  3402     if (t <= 0.33 + 0.01)
  3406     if (t <= 0.34 + 0.01)
  3410     if (t <= 0.35 + 0.01)
  3414     if (t <= 0.36 + 0.01)
  3418     if (t <= 0.37 + 0.01)
  3422     if (t <= 0.38 + 0.01)
  3426     if (t <= 0.39 + 0.01)
  3428         return  96.60999999999999;
  3430     if (t <= 0.40 + 0.01)
  3434     if (t <= 0.41 + 0.01)
  3436         return  75.67999999999999;
  3438     if (t <= 0.42 + 0.01)
  3440         return  64.56999999999999;
  3442     if (t <= 0.43 + 0.01)
  3446     if (t <= 0.44 + 0.01)
  3450     if (t <= 0.45 + 0.01)
  3454     if (t <= 0.46 + 0.01)
  3458     if (t <= 0.47 + 0.01)
  3460         return  34.83000000000001;
  3462     if (t <= 0.48 + 0.01)
  3466     if (t <= 0.49 + 0.01)
  3470     if (t <= 0.50 + 0.01)
  3474     if (t <= 0.51 + 0.01)
  3478     if (t <= 0.52 + 0.01)
  3482     if (t <= 0.53 + 0.01)
  3486     if (t <= 0.54 + 0.01)
  3490     if (t <= 0.55 + 0.01)
  3494     if (t <= 0.56 + 0.01)
  3498     if (t <= 0.57 + 0.01)
  3502     if (t <= 0.58 + 0.01)
  3506     if (t <= 0.59 + 0.01)
  3510     if (t <= 0.60 + 0.01)
  3514     if (t <= 0.61 + 0.01)
  3518     if (t <= 0.62 + 0.01)
  3522     if (t <= 0.63 + 0.01)
  3526     if (t <= 0.64 + 0.01)
  3530     if (t <= 0.65 + 0.01)
  3534     if (t <= 0.66 + 0.01)
  3538     if (t <= 0.67 + 0.01)
  3542     if (t <= 0.68 + 0.01)
  3546     if (t <= 0.69 + 0.01)
  3550     if (t <= 0.70 + 0.01)
  3554     if (t <= 0.71 + 0.01)
  3558     if (t <= 0.72 + 0.01)
  3562     if (t <= 0.73 + 0.01)
  3566     if (t <= 0.74 + 0.01)
  3570     if (t <= 0.75 + 0.01)
  3572         return  32.93000000000001;
  3574     if (t <= 0.76 + 0.01)
  3576         return  32.93000000000001;
  3578     if (t <= 0.77 + 0.01)
  3582     if (t <= 0.78 + 0.01)
  3584         return  33.01000000000001;
  3586     if (t <= 0.79 + 0.01)
  3594     if (t <= 0.00 + 0.01)
  3596         return    0.6737071250000001;
  3598     if (t <= 0.01 + 0.01)
  3600         return    0.6799071250000001;
  3602     if (t <= 0.02 + 0.01)
  3604         return    0.6897071250000001;
  3606     if (t <= 0.03 + 0.01)
  3608         return    0.6997071250000001;
  3610     if (t <= 0.04 + 0.01)
  3612         return    0.706407125       ;
  3614     if (t <= 0.05 + 0.01)
  3616         return    0.7084071250000001;
  3618     if (t <= 0.06 + 0.01)
  3620         return    0.7056071250000001;
  3622     if (t <= 0.07 + 0.01)
  3624         return    0.7049071250000001;
  3626     if (t <= 0.08 + 0.01)
  3628         return    0.7382071250000001;
  3630     if (t <= 0.09 + 0.01)
  3632         return    0.8778071249999999;
  3634     if (t <= 0.10 + 0.01)
  3636         return    1.209007125       ;
  3638     if (t <= 0.11 + 0.01)
  3640         return    1.746007125       ;
  3642     if (t <= 0.12 + 0.01)
  3644         return    2.368007125       ;
  3646     if (t <= 0.13 + 0.01)
  3648         return    2.897007125       ;
  3650     if (t <= 0.14 + 0.01)
  3652         return    3.258007125       ;
  3654     if (t <= 0.15 + 0.01)
  3656         return    3.518007125       ;
  3658     if (t <= 0.16 + 0.01)
  3660         return    3.752007125       ;
  3662     if (t <= 0.17 + 0.01)
  3664         return    3.936007124999999 ;
  3666     if (t <= 0.18 + 0.01)
  3668         return    3.957007125       ;
  3670     if (t <= 0.19 + 0.01)
  3672         return    3.716007125       ;
  3674     if (t <= 0.20 + 0.01)
  3676         return    3.212007125       ;
  3678     if (t <= 0.21 + 0.01)
  3680         return    2.551007125       ;
  3682     if (t <= 0.22 + 0.01)
  3684         return    1.898007125       ;
  3686     if (t <= 0.23 + 0.01)
  3688         return    1.395007125       ;
  3690     if (t <= 0.24 + 0.01)
  3692         return    1.107007125       ;
  3694     if (t <= 0.25 + 0.01)
  3696         return    1.009007125       ;
  3698     if (t <= 0.26 + 0.01)
  3700         return    1.028007125       ;
  3702     if (t <= 0.27 + 0.01)
  3704         return    1.082007125       ;
  3706     if (t <= 0.28 + 0.01)
  3708         return    1.114007125       ;
  3710     if (t <= 0.29 + 0.01)
  3712         return    1.111007125       ;
  3714     if (t <= 0.30 + 0.01)
  3716         return    1.089007125       ;
  3718     if (t <= 0.31 + 0.01)
  3720         return    1.071007125       ;
  3722     if (t <= 0.32 + 0.01)
  3724         return    1.067007125       ;
  3726     if (t <= 0.33 + 0.01)
  3728         return    1.068007125       ;
  3730     if (t <= 0.34 + 0.01)
  3732         return    1.056007125       ;
  3734     if (t <= 0.35 + 0.01)
  3736         return    1.021007125       ;
  3738     if (t <= 0.36 + 0.01)
  3740         return    0.955007125       ;
  3742     if (t <= 0.37 + 0.01)
  3744         return    0.856707125       ;
  3746     if (t <= 0.38 + 0.01)
  3748         return    0.7230071250000001;
  3750     if (t <= 0.39 + 0.01)
  3752         return    0.5555071250000001;
  3754     if (t <= 0.40 + 0.01)
  3756         return    0.3694071249999999;
  3758     if (t <= 0.41 + 0.01)
  3760         return    0.203607125       ;
  3762     if (t <= 0.42 + 0.01)
  3764         return    0.117107125       ;
  3766     if (t <= 0.43 + 0.01)
  3768         return    0.154707125       ;
  3770     if (t <= 0.44 + 0.01)
  3772         return    0.3019071249999999;
  3774     if (t <= 0.45 + 0.01)
  3776         return    0.4833071250000001;
  3778     if (t <= 0.46 + 0.01)
  3780         return    0.6279071250000001;
  3782     if (t <= 0.47 + 0.01)
  3784         return    0.7323071250000001;
  3786     if (t <= 0.48 + 0.01)
  3788         return    0.8435071250000001;
  3790     if (t <= 0.49 + 0.01)
  3792         return    0.989007125       ;
  3794     if (t <= 0.50 + 0.01)
  3796         return    1.137007125       ;
  3798     if (t <= 0.51 + 0.01)
  3800         return    1.223007125       ;
  3802     if (t <= 0.52 + 0.01)
  3804         return    1.206007125       ;
  3806     if (t <= 0.53 + 0.01)
  3808         return    1.096007125       ;
  3810     if (t <= 0.54 + 0.01)
  3812         return    0.9400071249999998;
  3814     if (t <= 0.55 + 0.01)
  3816         return    0.7924071250000001;
  3818     if (t <= 0.56 + 0.01)
  3820         return    0.6857071250000001;
  3822     if (t <= 0.57 + 0.01)
  3824         return    0.6278071250000001;
  3826     if (t <= 0.58 + 0.01)
  3828         return    0.6083071250000001;
  3830     if (t <= 0.59 + 0.01)
  3832         return    0.610607125       ;
  3834     if (t <= 0.60 + 0.01)
  3836         return    0.620007125       ;
  3838     if (t <= 0.61 + 0.01)
  3840         return    0.6279071250000001;
  3842     if (t <= 0.62 + 0.01)
  3844         return    0.6320071250000001;
  3846     if (t <= 0.63 + 0.01)
  3848         return    0.6343071250000001;
  3850     if (t <= 0.64 + 0.01)
  3852         return    0.6385071250000001;
  3854     if (t <= 0.65 + 0.01)
  3856         return    0.6476071250000001;
  3858     if (t <= 0.66 + 0.01)
  3860         return    0.6620071250000001;
  3862     if (t <= 0.67 + 0.01)
  3864         return    0.679107125       ;
  3866     if (t <= 0.68 + 0.01)
  3868         return    0.6948071250000001;
  3870     if (t <= 0.69 + 0.01)
  3872         return    0.7049071250000001;
  3874     if (t <= 0.70 + 0.01)
  3876         return    0.7071071250000001;
  3878     if (t <= 0.71 + 0.01)
  3880         return    0.7021071250000001;
  3882     if (t <= 0.72 + 0.01)
  3884         return    0.6926071250000001;
  3886     if (t <= 0.73 + 0.01)
  3888         return    0.6825071250000001;
  3890     if (t <= 0.74 + 0.01)
  3892         return    0.675007125       ;
  3894     if (t <= 0.75 + 0.01)
  3896         return    0.671207125       ;
  3898     if (t <= 0.76 + 0.01)
  3900         return    0.6702071250000001;
  3902     if (t <= 0.77 + 0.01)
  3904         return    0.6702071250000001;
  3906     if (t <= 0.78 + 0.01)
  3908         return    0.669507125       ;
  3910     if (t <= 0.79 + 0.01)
  3912         return    0.6674071250000001;
  3927     if (t <=  0.0050 + 0.005)
  3931     if (t <=  0.0100 + 0.005)
  3935     if (t <=  0.0150 + 0.005)
  3939     if (t <=  0.0200 + 0.005)
  3943     if (t <=  0.0250 + 0.005)
  3947     if (t <=  0.0300 + 0.005)
  3951     if (t <=  0.0350 + 0.005)
  3955     if (t <=  0.0400 + 0.005)
  3959     if (t <=  0.0450 + 0.005)
  3963     if (t <=  0.0500 + 0.005)
  3967     if (t <=  0.0550 + 0.005)
  3971     if (t <=  0.0600 + 0.005)
  3975     if (t <=  0.0650 + 0.005)
  3979     if (t <=  0.0700 + 0.005)
  3983     if (t <=  0.0750 + 0.005)
  3987     if (t <=  0.0800 + 0.005)
  3991     if (t <=  0.0850 + 0.005)
  3995     if (t <=  0.0900 + 0.005)
  3999     if (t <=  0.0950 + 0.005)
  4003     if (t <=  0.1000 + 0.005)
  4007     if (t <=  0.1050 + 0.005)
  4011     if (t <=  0.1100 + 0.005)
  4015     if (t <=  0.1150 + 0.005)
  4019     if (t <=  0.1200 + 0.005)
  4023     if (t <=  0.1250 + 0.005)
  4027     if (t <=  0.1300 + 0.005)
  4031     if (t <=  0.1350 + 0.005)
  4035     if (t <=  0.1400 + 0.005)
  4039     if (t <=  0.1450 + 0.005)
  4043     if (t <=  0.1500 + 0.005)
  4047     if (t <=  0.1550 + 0.005)
  4051     if (t <=  0.1600 + 0.005)
  4055     if (t <=  0.1650 + 0.005)
  4059     if (t <=  0.1700 + 0.005)
  4063     if (t <=  0.1750 + 0.005)
  4067     if (t <=  0.1800 + 0.005)
  4071     if (t <=  0.1850 + 0.005)
  4075     if (t <=  0.1900 + 0.005)
  4079     if (t <=  0.1950 + 0.005)
  4083     if (t <=  0.2000 + 0.005)
  4087     if (t <=  0.2050 + 0.005)
  4091     if (t <=  0.2100 + 0.005)
  4095     if (t <=  0.2150 + 0.005)
  4099     if (t <=  0.2200 + 0.005)
  4103     if (t <=  0.2250 + 0.005)
  4107     if (t <=  0.2300 + 0.005)
  4111     if (t <=  0.2350 + 0.005)
  4115     if (t <=  0.2400 + 0.005)
  4119     if (t <=  0.2450 + 0.005)
  4123     if (t <=  0.2500 + 0.005)
  4127     if (t <=  0.2550 + 0.005)
  4131     if (t <=  0.2600 + 0.005)
  4135     if (t <=  0.2650 + 0.005)
  4139     if (t <=  0.2700 + 0.005)
  4143     if (t <=  0.2750 + 0.005)
  4147     if (t <=  0.2800 + 0.005)
  4151     if (t <=  0.2850 + 0.005)
  4155     if (t <=  0.2900 + 0.005)
  4159     if (t <=  0.2950 + 0.005)
  4163     if (t <=  0.3000 + 0.005)
  4167     if (t <=  0.3050 + 0.005)
  4171     if (t <=  0.3100 + 0.005)
  4175     if (t <=  0.3150 + 0.005)
  4179     if (t <=  0.3200 + 0.005)
  4183     if (t <=  0.3250 + 0.005)
  4187     if (t <=  0.3300 + 0.005)
  4191     if (t <=  0.3350 + 0.005)
  4195     if (t <=  0.3400 + 0.005)
  4199     if (t <=  0.3450 + 0.005)
  4203     if (t <=  0.3500 + 0.005)
  4207     if (t <=  0.3550 + 0.005)
  4211     if (t <=  0.3600 + 0.005)
  4215     if (t <=  0.3650 + 0.005)
  4219     if (t <=  0.3700 + 0.005)
  4223     if (t <=  0.3750 + 0.005)
  4227     if (t <=  0.3800 + 0.005)
  4231     if (t <=  0.3850 + 0.005)
  4235     if (t <=  0.3900 + 0.005)
  4239     if (t <=  0.3950 + 0.005)
  4243     if (t <=  0.4000 + 0.005)
  4247     if (t <=  0.4050 + 0.005)
  4251     if (t <=  0.4100 + 0.005)
  4255     if (t <=  0.4150 + 0.005)
  4259     if (t <=  0.4200 + 0.005)
  4263     if (t <=  0.4250 + 0.005)
  4267     if (t <=  0.4300 + 0.005)
  4271     if (t <=  0.4350 + 0.005)
  4275     if (t <=  0.4400 + 0.005)
  4279     if (t <=  0.4450 + 0.005)
  4283     if (t <=  0.4500 + 0.005)
  4287     if (t <=  0.4550 + 0.005)
  4291     if (t <=  0.4600 + 0.005)
  4295     if (t <=  0.4650 + 0.005)
  4299     if (t <=  0.4700 + 0.005)
  4303     if (t <=  0.4750 + 0.005)
  4307     if (t <=  0.4800 + 0.005)
  4311     if (t <=  0.4850 + 0.005)
  4315     if (t <=  0.4900 + 0.005)
  4319     if (t <=  0.4950 + 0.005)
  4323     if (t <=  0.5000 + 0.005)
  4327     if (t <=  0.5050 + 0.005)
  4331     if (t <=  0.5100 + 0.005)
  4335     if (t <=  0.5150 + 0.005)
  4339     if (t <=  0.5200 + 0.005)
  4343     if (t <=  0.5250 + 0.005)
  4347     if (t <=  0.5300 + 0.005)
  4351     if (t <=  0.5350 + 0.005)
  4355     if (t <=  0.5400 + 0.005)
  4359     if (t <=  0.5450 + 0.005)
  4363     if (t <=  0.5500 + 0.005)
  4367     if (t <=  0.5550 + 0.005)
  4371     if (t <=  0.5600 + 0.005)
  4375     if (t <=  0.5650 + 0.005)
  4379     if (t <=  0.5700 + 0.005)
  4383     if (t <=  0.5750 + 0.005)
  4387     if (t <=  0.5800 + 0.005)
  4391     if (t <=  0.5850 + 0.005)
  4395     if (t <=  0.5900 + 0.005)
  4399     if (t <=  0.5950 + 0.005)
  4403     if (t <=  0.6000 + 0.005)
  4407     if (t <=  0.6050 + 0.005)
  4411     if (t <=  0.6100 + 0.005)
  4415     if (t <=  0.6150 + 0.005)
  4419     if (t <=  0.6200 + 0.005)
  4423     if (t <=  0.6250 + 0.005)
  4427     if (t <=  0.6300 + 0.005)
  4431     if (t <=  0.6350 + 0.005)
  4435     if (t <=  0.6400 + 0.005)
  4439     if (t <=  0.6450 + 0.005)
  4443     if (t <=  0.6500 + 0.005)
  4447     if (t <=  0.6550 + 0.005)
  4451     if (t <=  0.6600 + 0.005)
  4455     if (t <=  0.6650 + 0.005)
  4459     if (t <=  0.6700 + 0.005)
  4463     if (t <=  0.6750 + 0.005)
  4467     if (t <=  0.6800 + 0.005)
  4471     if (t <=  0.6850 + 0.005)
  4475     if (t <=  0.6900 + 0.005)
  4479     if (t <=  0.6950 + 0.005)
  4483     if (t <=  0.7000 + 0.005)
  4487     if (t <=  0.7050 + 0.005)
  4491     if (t <=  0.7100 + 0.005)
  4495     if (t <=  0.7150 + 0.005)
  4499     if (t <=  0.7200 + 0.005)
  4503     if (t <=  0.7250 + 0.005)
  4507     if (t <=  0.7300 + 0.005)
  4511     if (t <=  0.7350 + 0.005)
  4515     if (t <=  0.7400 + 0.005)
  4519     if (t <=  0.7450 + 0.005)
  4523     if (t <=  0.7500 + 0.005)
  4527     if (t <=  0.7550 + 0.005)
  4531     if (t <=  0.7600 + 0.005)
  4535     if (t <=  0.7650 + 0.005)
  4539     if (t <=  0.7700 + 0.005)
  4543     if (t <=  0.7750 + 0.005)
  4547     if (t <=  0.7800 + 0.005)
  4551     if (t <=  0.7850 + 0.005)
  4555     if (t <=  0.7900 + 0.005)
  4559     if (t <=  0.7950 + 0.005)
  4563     if (t <=  0.8000 + 0.005)
  4567     if (t <=  0.8050 + 0.005)
  4571     if (t <=  0.8100 + 0.005)
  4575     if (t <=  0.8150 + 0.005)
  4579     if (t <=  0.8200 + 0.005)
  4583     if (t <=  0.8250 + 0.005)
  4587     if (t <=  0.8300 + 0.005)
  4591     if (t <=  0.8350 + 0.005)
  4595     if (t <=  0.8400 + 0.005)
  4599     if (t <=  0.8450 + 0.005)
  4603     if (t <=  0.8500 + 0.005)
  4607     if (t <=  0.8550 + 0.005)
  4611     if (t <=  0.8600 + 0.005)
  4615     if (t <=  0.8650 + 0.005)
  4619     if (t <=  0.8700 + 0.005)
  4623     if (t <=  0.8750 + 0.005)
  4627     if (t <=  0.8800 + 0.005)
  4631     if (t <=  0.8850 + 0.005)
  4635     if (t <=  0.8900 + 0.005)
  4639     if (t <=  0.8950 + 0.005)
  4643     if (t <=  0.9000 + 0.005)
  4647     if (t <=  0.9050 + 0.005)
  4651     if (t <=  0.9100 + 0.005)
  4655     if (t <=  0.9150 + 0.005)
  4659     if (t <=  0.9200 + 0.005)
  4663     if (t <=  0.9250 + 0.005)
  4667     if (t <=  0.9300 + 0.005)
  4671     if (t <=  0.9350 + 0.005)
  4675     if (t <=  0.9400 + 0.005)
  4679     if (t <=  0.9450 + 0.005)
  4683     if (t <=  0.9500 + 0.005)
  4687     if (t <=  0.9550 + 0.005)
  4691     if (t <=  0.9600 + 0.005)
  4695     if (t <=  0.9650 + 0.005)
  4699     if (t <=  0.9700 + 0.005)
  4703     if (t <=  0.9750 + 0.005)
  4707     if (t <=  0.9800 + 0.005)
  4711     if (t <=  0.9850 + 0.005)
  4715     if (t <=  0.9900 + 0.005)
  4719     if (t <=  0.9950 + 0.005)
  4723     if (t <=  1.0000 + 0.005)
  4727     if (t <=  1.0050 + 0.005)
  4731     if (t <=  1.0100 + 0.005)
  4735     if (t <=  1.0150 + 0.005)
  4739     if (t <=  1.0200 + 0.005)
  4743     if (t <=  1.0250 + 0.005)
  4747     if (t <=  1.0300 + 0.005)
  4751     if (t <=  1.0350 + 0.005)
  4755     if (t <=  1.0400 + 0.005)
  4759     if (t <=  1.0450 + 0.005)
  4763     if (t <=  1.0500 + 0.005)
  4767     if (t <=  1.0550 + 0.005)
  4771     if (t <=  1.0600 + 0.005)
  4775     if (t <=  1.0650 + 0.005)
  4779     if (t <=  1.0700 + 0.005)
  4783     if (t <=  1.0750 + 0.005)
  4787     if (t <=  1.0800 + 0.005)
  4791     if (t <=  1.0850 + 0.005)
  4795     if (t <=  1.0900 + 0.005)
  4799     if (t <=  1.0950 + 0.005)
  4803     if (t <=  1.1000 + 0.005)
  4807     if (t <=  1.1050 + 0.005)
  4811     if (t <=  1.1100 + 0.005)
  4824         tNew = ( ( (
int) floor (t * 10000) ) % 11000) / 10000.0;
  4827     Real dt =  
static_cast<
int> ( std::floor ( std::ceil ( tNew * 1000000.0 ) / 1000.0 ) ) % 5 / 1000.0;
  4828     Real ti = floor (tNew * 1000) / 1000.0 - dt;
  4829     Real tii = ti + 0.005;
  4830     return - ( (inletCylinder (tii, 0, 0, 0, 0) - inletCylinder (ti, 0, 0, 0, 0) ) / (0.005) * (tNew - (ti) ) + inletCylinder (ti, 0, 0, 0, 0) );
  4846     Real area (0.5 * 0.5 * 3.141592653589793); 
  4848     Real radiusSquared (0);
  4849     radiusSquared = 0.5 * 0.5; 
  4852     peak = ( 2.0 * flux ) / ( area );
  4857             return n1 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4859             return n2 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4861             return n3 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4880     Real area (0.5 * 0.5 * 3.141592653589793); 
  4882     Real radiusSquared (0);
  4883     radiusSquared = 0.5 * 0.5; 
  4886     peak = ( 2.0 * flux ) / ( area );
  4891             return n1 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4893             return n2 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4895             return n3 * ( peak * std::max ( 0.0, ( (radiusSquared - ( (x - x0) * (x - x0) + (y - y0) * (y - y0) ) ) / radiusSquared) ) );
  4909     Real vBar ( -15.0 * std::sin ( t * 2 * 3.14 / 0.1  ) );
  4926     Real newt = ( (
Real) ( ( (
int) (t * 1000) ) % 792) ) / 1000;
  4927     if ( newt <= 0.    + 0.004 )
  4931     else if ( newt <= 0.004 + 0.004 )
  4935     else if ( newt <= 0.008 + 0.004 )
  4939     else if ( newt <= 0.012 + 0.004 )
  4943     else if ( newt <= 0.016 + 0.004 )
  4947     else if ( newt <= 0.02  + 0.004 )
  4951     else if ( newt <= 0.024 + 0.004 )
  4955     else if ( newt <= 0.028 + 0.004 )
  4959     else if ( newt <= 0.032 + 0.004 )
  4963     else if ( newt <= 0.036 + 0.004 )
  4967     else if ( newt <= 0.04  + 0.004 )
  4971     else if ( newt <= 0.044 + 0.004 )
  4975     else if ( newt <= 0.048 + 0.004 )
  4979     else if ( newt <= 0.052 + 0.004 )
  4983     else if ( newt <= 0.056 + 0.004 )
  4987     else if ( newt <= 0.06  + 0.004 )
  4991     else if ( newt <= 0.064 + 0.004 )
  4995     else if ( newt <= 0.068 + 0.004 )
  4999     else if ( newt <= 0.072 + 0.004 )
  5003     else if ( newt <= 0.076 + 0.004 )
  5007     else if ( newt <= 0.08  + 0.004 )
  5011     else if ( newt <= 0.084 + 0.004 )
  5015     else if ( newt <= 0.088 + 0.004 )
  5019     else if ( newt <= 0.092 + 0.004 )
  5023     else if ( newt <= 0.096 + 0.004 )
  5027     else if ( newt <= 0.1   + 0.004 )
  5031     else if ( newt <= 0.104 + 0.004 )
  5035     else if ( newt <= 0.108 + 0.004 )
  5039     else if ( newt <= 0.112 + 0.004 )
  5043     else if ( newt <= 0.116 + 0.004 )
  5047     else if ( newt <= 0.12  + 0.004 )
  5051     else if ( newt <= 0.124 + 0.004 )
  5055     else if ( newt <= 0.128 + 0.004 )
  5059     else if ( newt <= 0.132 + 0.004 )
  5063     else if ( newt <= 0.136 + 0.004 )
  5067     else if ( newt <= 0.14  + 0.004 )
  5071     else if ( newt <= 0.144 + 0.004 )
  5075     else if ( newt <= 0.148 + 0.004 )
  5079     else if ( newt <= 0.152 + 0.004 )
  5083     else if ( newt <= 0.156 + 0.004 )
  5087     else if ( newt <= 0.16  + 0.004 )
  5091     else if ( newt <= 0.164 + 0.004 )
  5095     else if ( newt <= 0.168 + 0.004 )
  5099     else if ( newt <= 0.172 + 0.004 )
  5103     else if ( newt <= 0.176 + 0.004 )
  5107     else if ( newt <= 0.18  + 0.004 )
  5111     else if ( newt <= 0.184 + 0.004 )
  5115     else if ( newt <= 0.188 + 0.004 )
  5119     else if ( newt <= 0.192 + 0.004 )
  5123     else if ( newt <= 0.196 + 0.004 )
  5127     else if ( newt <= 0.2   + 0.004 )
  5131     else if ( newt <= 0.204 + 0.004 )
  5135     else if ( newt <= 0.208 + 0.004 )
  5139     else if ( newt <= 0.212 + 0.004 )
  5143     else if ( newt <= 0.216 + 0.004 )
  5147     else if ( newt <= 0.22  + 0.004 )
  5151     else if ( newt <= 0.224 + 0.004 )
  5155     else if ( newt <= 0.228 + 0.004 )
  5159     else if ( newt <= 0.232 + 0.004 )
  5163     else if ( newt <= 0.236 + 0.004 )
  5167     else if ( newt <= 0.24  + 0.004 )
  5171     else if ( newt <= 0.244 + 0.004 )
  5175     else if ( newt <= 0.248 + 0.004 )
  5179     else if ( newt <= 0.252 + 0.004 )
  5183     else if ( newt <= 0.256 + 0.004 )
  5187     else if ( newt <= 0.26  + 0.004 )
  5191     else if ( newt <= 0.264 + 0.004 )
  5195     else if ( newt <= 0.268 + 0.004 )
  5199     else if ( newt <= 0.272 + 0.004 )
  5203     else if ( newt <= 0.276 + 0.004 )
  5207     else if ( newt <= 0.28  + 0.004 )
  5211     else if ( newt <= 0.284 + 0.004 )
  5215     else if ( newt <= 0.288 + 0.004 )
  5219     else if ( newt <= 0.292 + 0.004 )
  5223     else if ( newt <= 0.296 + 0.004 )
  5227     else if ( newt <= 0.3   + 0.004 )
  5231     else if ( newt <= 0.304 + 0.004 )
  5235     else if ( newt <= 0.308 + 0.004 )
  5239     else if ( newt <= 0.312 + 0.004 )
  5243     else if ( newt <= 0.316 + 0.004 )
  5247     else if ( newt <= 0.32  + 0.004 )
  5251     else if ( newt <= 0.324 + 0.004 )
  5255     else if ( newt <= 0.328 + 0.004 )
  5259     else if ( newt <= 0.332 + 0.004 )
  5263     else if ( newt <= 0.336 + 0.004 )
  5267     else if ( newt <= 0.34  + 0.004 )
  5271     else if ( newt <= 0.344 + 0.004 )
  5275     else if ( newt <= 0.348 + 0.004 )
  5279     else if ( newt <= 0.352 + 0.004 )
  5283     else if ( newt <= 0.356 + 0.004 )
  5287     else if ( newt <= 0.36  + 0.004 )
  5291     else if ( newt <= 0.364 + 0.004 )
  5295     else if ( newt <= 0.368 + 0.004 )
  5299     else if ( newt <= 0.372 + 0.004 )
  5303     else if ( newt <= 0.376 + 0.004 )
  5307     else if ( newt <= 0.38  + 0.004 )
  5311     else if ( newt <= 0.384 + 0.004 )
  5315     else if ( newt <= 0.388 + 0.004 )
  5319     else if ( newt <= 0.392 + 0.004 )
  5323     else if ( newt <= 0.396 + 0.004 )
  5327     else if ( newt <= 0.4   + 0.004 )
  5331     else if ( newt <= 0.404 + 0.004 )
  5335     else if ( newt <= 0.408 + 0.004 )
  5339     else if ( newt <= 0.412 + 0.004 )
  5343     else if ( newt <= 0.416 + 0.004 )
  5347     else if ( newt <= 0.42  + 0.004 )
  5351     else if ( newt <= 0.424 + 0.004 )
  5355     else if ( newt <= 0.428 + 0.004 )
  5359     else if ( newt <= 0.432 + 0.004 )
  5363     else if ( newt <= 0.436 + 0.004 )
  5367     else if ( newt <= 0.44  + 0.004 )
  5371     else if ( newt <= 0.444 + 0.004 )
  5375     else if ( newt <= 0.448 + 0.004 )
  5379     else if ( newt <= 0.452 + 0.004 )
  5383     else if ( newt <= 0.456 + 0.004 )
  5387     else if ( newt <= 0.46  + 0.004 )
  5391     else if ( newt <= 0.464 + 0.004 )
  5395     else if ( newt <= 0.468 + 0.004 )
  5399     else if ( newt <= 0.472 + 0.004 )
  5403     else if ( newt <= 0.476 + 0.004 )
  5407     else if ( newt <= 0.48  + 0.004 )
  5411     else if ( newt <= 0.484 + 0.004 )
  5415     else if ( newt <= 0.488 + 0.004 )
  5419     else if ( newt <= 0.492 + 0.004 )
  5423     else if ( newt <= 0.496 + 0.004 )
  5427     else if ( newt <= 0.5   + 0.004 )
  5431     else if ( newt <= 0.504 + 0.004 )
  5435     else if ( newt <= 0.508 + 0.004 )
  5439     else if ( newt <= 0.512 + 0.004 )
  5443     else if ( newt <= 0.516 + 0.004 )
  5447     else if ( newt <= 0.52  + 0.004 )
  5451     else if ( newt <= 0.524 + 0.004 )
  5455     else if ( newt <= 0.528 + 0.004 )
  5459     else if ( newt <= 0.532 + 0.004 )
  5463     else if ( newt <= 0.536 + 0.004 )
  5467     else if ( newt <= 0.54  + 0.004 )
  5471     else if ( newt <= 0.544 + 0.004 )
  5475     else if ( newt <= 0.548 + 0.004 )
  5479     else if ( newt <= 0.552 + 0.004 )
  5483     else if ( newt <= 0.556 + 0.004 )
  5487     else if ( newt <= 0.56  + 0.004 )
  5491     else if ( newt <= 0.564 + 0.004 )
  5495     else if ( newt <= 0.568 + 0.004 )
  5499     else if ( newt <= 0.572 + 0.004 )
  5503     else if ( newt <= 0.576 + 0.004 )
  5507     else if ( newt <= 0.58  + 0.004 )
  5511     else if ( newt <= 0.584 + 0.004 )
  5515     else if ( newt <= 0.588 + 0.004 )
  5519     else if ( newt <= 0.592 + 0.004 )
  5523     else if ( newt <= 0.596 + 0.004 )
  5527     else if ( newt <= 0.6   + 0.004 )
  5531     else if ( newt <= 0.604 + 0.004 )
  5535     else if ( newt <= 0.608 + 0.004 )
  5539     else if ( newt <= 0.612 + 0.004 )
  5543     else if ( newt <= 0.616 + 0.004 )
  5547     else if ( newt <= 0.62  + 0.004 )
  5551     else if ( newt <= 0.624 + 0.004 )
  5555     else if ( newt <= 0.628 + 0.004 )
  5559     else if ( newt <= 0.632 + 0.004 )
  5563     else if ( newt <= 0.636 + 0.004 )
  5567     else if ( newt <= 0.64  + 0.004 )
  5571     else if ( newt <= 0.644 + 0.004 )
  5575     else if ( newt <= 0.648 + 0.004 )
  5579     else if ( newt <= 0.652 + 0.004 )
  5583     else if ( newt <= 0.656 + 0.004 )
  5587     else if ( newt <= 0.66  + 0.004 )
  5591     else if ( newt <= 0.664 + 0.004 )
  5595     else if ( newt <= 0.668 + 0.004 )
  5599     else if ( newt <= 0.672 + 0.004 )
  5603     else if ( newt <= 0.676 + 0.004 )
  5607     else if ( newt <= 0.68  + 0.004 )
  5611     else if ( newt <= 0.684 + 0.004 )
  5615     else if ( newt <= 0.688 + 0.004 )
  5619     else if ( newt <= 0.692 + 0.004 )
  5623     else if ( newt <= 0.696 + 0.004 )
  5627     else if ( newt <= 0.7   + 0.004 )
  5631     else if ( newt <= 0.704 + 0.004 )
  5635     else if ( newt <= 0.708 + 0.004 )
  5639     else if ( newt <= 0.712 + 0.004 )
  5643     else if ( newt <= 0.716 + 0.004 )
  5647     else if ( newt <= 0.72  + 0.004 )
  5651     else if ( newt <= 0.724 + 0.004 )
  5655     else if ( newt <= 0.728 + 0.004 )
  5659     else if ( newt <= 0.732 + 0.004 )
  5663     else if ( newt <= 0.736 + 0.004 )
  5667     else if ( newt <= 0.74  + 0.004 )
  5671     else if ( newt <= 0.744 + 0.004 )
  5675     else if ( newt <= 0.748 + 0.004 )
  5679     else if ( newt <= 0.752 + 0.004 )
  5683     else if ( newt <= 0.756 + 0.004 )
  5687     else if ( newt <= 0.76  + 0.004 )
  5691     else if ( newt <= 0.764 + 0.004 )
  5695     else if ( newt <= 0.768 + 0.004 )
  5699     else if ( newt <= 0.772 + 0.004 )
  5703     else if ( newt <= 0.776 + 0.004 )
  5707     else if ( newt <= 0.78  + 0.004 )
  5711     else if ( newt <= 0.784 + 0.004 )
  5715     else if ( newt <= 0.788 + 0.004 )
  5719     else if ( newt <= 0.792 + 0.004 )
  5728     Real newt = ( (
Real) ( ( (
int) (t * 1000) ) % 792) ) / 1000;
  5729     if ( newt <= 0.    + 0.004 )
  5731         return  1.96970000e-02   * (-10);
  5733     else if ( newt <= 0.004 + 0.004 )
  5735         return  1.93150000e-02   * (-10);
  5737     else if ( newt <= 0.008 + 0.004 )
  5739         return  2.12310000e-02   * (-10);
  5741     else if ( newt <= 0.012 + 0.004 )
  5743         return  3.04300000e-02   * (-10);
  5745     else if ( newt <= 0.016 + 0.004 )
  5747         return  4.50000000e-02   * (-10);
  5749     else if ( newt <= 0.02  + 0.004 )
  5751         return  6.18670000e-02   * (-10);
  5753     else if ( newt <= 0.024 + 0.004 )
  5755         return  8.48730000e-02   * (-10);
  5757     else if ( newt <= 0.028 + 0.004 )
  5759         return  1.07874000e-01   * (-10);
  5761     else if ( newt <= 0.032 + 0.004 )
  5763         return  1.32410000e-01   * (-10);
  5765     else if ( newt <= 0.036 + 0.004 )
  5767         return  1.57714000e-01   * (-10);
  5769     else if ( newt <= 0.04  + 0.004 )
  5771         return  1.86849000e-01   * (-10);
  5773     else if ( newt <= 0.044 + 0.004 )
  5775         return  2.19821000e-01   * (-10);
  5777     else if ( newt <= 0.048 + 0.004 )
  5779         return  2.52788000e-01   * (-10);
  5781     else if ( newt <= 0.052 + 0.004 )
  5783         return  2.85760000e-01   * (-10);
  5785     else if ( newt <= 0.056 + 0.004 )
  5787         return  3.17198000e-01   * (-10);
  5789     else if ( newt <= 0.06  + 0.004 )
  5791         return  3.44798000e-01   * (-10);
  5793     else if ( newt <= 0.064 + 0.004 )
  5795         return  3.70869000e-01   * (-10);
  5797     else if ( newt <= 0.068 + 0.004 )
  5799         return  3.96959000e-01   * (-10);
  5801     else if ( newt <= 0.072 + 0.004 )
  5803         return  4.26080000e-01   * (-10);
  5805     else if ( newt <= 0.076 + 0.004 )
  5807         return  4.53675000e-01   * (-10);
  5809     else if ( newt <= 0.08  + 0.004 )
  5811         return  4.82034000e-01   * (-10);
  5813     else if ( newt <= 0.084 + 0.004 )
  5815         return  5.09630000e-01   * (-10);
  5817     else if ( newt <= 0.088 + 0.004 )
  5819         return  5.27264000e-01   * (-10);
  5821     else if ( newt <= 0.092 + 0.004 )
  5823         return  5.38798000e-01   * (-10);
  5825     else if ( newt <= 0.096 + 0.004 )
  5827         return  5.48759000e-01   * (-10);
  5829     else if ( newt <= 0.1   + 0.004 )
  5831         return  5.52572000e-01   * (-10);
  5833     else if ( newt <= 0.104 + 0.004 )
  5835         return  5.63344000e-01   * (-10);
  5837     else if ( newt <= 0.108 + 0.004 )
  5839         return  5.80978000e-01   * (-10);
  5841     else if ( newt <= 0.112 + 0.004 )
  5843         return  6.00138000e-01   * (-10);
  5845     else if ( newt <= 0.116 + 0.004 )
  5847         return  6.19298000e-01   * (-10);
  5849     else if ( newt <= 0.12  + 0.004 )
  5851         return  6.32310000e-01   * (-10);
  5853     else if ( newt <= 0.124 + 0.004 )
  5855         return  6.40746000e-01   * (-10);
  5857     else if ( newt <= 0.128 + 0.004 )
  5859         return  6.41508000e-01   * (-10);
  5861     else if ( newt <= 0.132 + 0.004 )
  5863         return  6.43081000e-01   * (-10);
  5865     else if ( newt <= 0.136 + 0.004 )
  5867         return  6.49944000e-01   * (-10);
  5869     else if ( newt <= 0.14  + 0.004 )
  5871         return  6.59190000e-01   * (-10);
  5873     else if ( newt <= 0.144 + 0.004 )
  5875         return  6.64528000e-01   * (-10);
  5877     else if ( newt <= 0.148 + 0.004 )
  5879         return  6.60716000e-01   * (-10);
  5881     else if ( newt <= 0.152 + 0.004 )
  5883         return  6.56092000e-01   * (-10);
  5885     else if ( newt <= 0.156 + 0.004 )
  5887         return  6.59190000e-01   * (-10);
  5889     else if ( newt <= 0.16  + 0.004 )
  5891         return  6.62241000e-01   * (-10);
  5893     else if ( newt <= 0.164 + 0.004 )
  5895         return  6.62241000e-01   * (-10);
  5897     else if ( newt <= 0.168 + 0.004 )
  5899         return  6.55330000e-01   * (-10);
  5901     else if ( newt <= 0.172 + 0.004 )
  5903         return  6.46131000e-01   * (-10);
  5905     else if ( newt <= 0.176 + 0.004 )
  5907         return  6.35408000e-01   * (-10);
  5909     else if ( newt <= 0.18  + 0.004 )
  5911         return  6.23111000e-01   * (-10);
  5913     else if ( newt <= 0.184 + 0.004 )
  5915         return  6.17010000e-01   * (-10);
  5917     else if ( newt <= 0.188 + 0.004 )
  5919         return  6.10862000e-01   * (-10);
  5921     else if ( newt <= 0.192 + 0.004 )
  5923         return  5.99376000e-01   * (-10);
  5925     else if ( newt <= 0.196 + 0.004 )
  5927         return  5.86316000e-01   * (-10);
  5929     else if ( newt <= 0.2   + 0.004 )
  5931         return  5.70970000e-01   * (-10);
  5933     else if ( newt <= 0.204 + 0.004 )
  5935         return  5.54145000e-01   * (-10);
  5937     else if ( newt <= 0.208 + 0.004 )
  5939         return  5.43374000e-01   * (-10);
  5941     else if ( newt <= 0.212 + 0.004 )
  5943         return  5.35700000e-01   * (-10);
  5945     else if ( newt <= 0.216 + 0.004 )
  5947         return  5.26502000e-01   * (-10);
  5949     else if ( newt <= 0.22  + 0.004 )
  5951         return  5.11965000e-01   * (-10);
  5953     else if ( newt <= 0.224 + 0.004 )
  5955         return  4.91232000e-01   * (-10);
  5957     else if ( newt <= 0.228 + 0.004 )
  5959         return  4.69022000e-01   * (-10);
  5961     else if ( newt <= 0.232 + 0.004 )
  5963         return  4.46002000e-01   * (-10);
  5965     else if ( newt <= 0.236 + 0.004 )
  5967         return  4.27605000e-01   * (-10);
  5969     else if ( newt <= 0.24  + 0.004 )
  5971         return  4.10733000e-01   * (-10);
  5973     else if ( newt <= 0.244 + 0.004 )
  5975         return  3.86950000e-01   * (-10);
  5977     else if ( newt <= 0.248 + 0.004 )
  5979         return  3.65502000e-01   * (-10);
  5981     else if ( newt <= 0.252 + 0.004 )
  5983         return  3.42501000e-01   * (-10);
  5985     else if ( newt <= 0.256 + 0.004 )
  5987         return  3.24862000e-01   * (-10);
  5989     else if ( newt <= 0.26  + 0.004 )
  5991         return  3.07994000e-01   * (-10);
  5993     else if ( newt <= 0.264 + 0.004 )
  5995         return  2.88058000e-01   * (-10);
  5997     else if ( newt <= 0.268 + 0.004 )
  5999         return  2.69656000e-01   * (-10);
  6001     else if ( newt <= 0.272 + 0.004 )
  6003         return  2.47422000e-01   * (-10);
  6005     else if ( newt <= 0.276 + 0.004 )
  6007         return  2.23653000e-01   * (-10);
  6009     else if ( newt <= 0.28  + 0.004 )
  6011         return  2.04484000e-01   * (-10);
  6013     else if ( newt <= 0.284 + 0.004 )
  6015         return  1.84547000e-01   * (-10);
  6017     else if ( newt <= 0.288 + 0.004 )
  6019         return  1.67680000e-01   * (-10);
  6021     else if ( newt <= 0.292 + 0.004 )
  6023         return  1.54644000e-01   * (-10);
  6025     else if ( newt <= 0.296 + 0.004 )
  6027         return  1.42376000e-01   * (-10);
  6029     else if ( newt <= 0.3   + 0.004 )
  6031         return  1.32410000e-01   * (-10);
  6033     else if ( newt <= 0.304 + 0.004 )
  6035         return  1.20142000e-01   * (-10);
  6037     else if ( newt <= 0.308 + 0.004 )
  6039         return  1.06340000e-01   * (-10);
  6041     else if ( newt <= 0.312 + 0.004 )
  6043         return  9.25370000e-02   * (-10);
  6045     else if ( newt <= 0.316 + 0.004 )
  6047         return  7.64370000e-02   * (-10);
  6049     else if ( newt <= 0.32  + 0.004 )
  6051         return  5.88030000e-02   * (-10);
  6053     else if ( newt <= 0.324 + 0.004 )
  6055         return  9.73060000e-03   * (-10);
  6057     else if ( newt <= 0.328 + 0.004 )
  6059         return  -3.09101000e-02  * (-10);
  6061     else if ( newt <= 0.332 + 0.004 )
  6063         return  -6.84790000e-02  * (-10);
  6065     else if ( newt <= 0.336 + 0.004 )
  6067         return  -1.58956600e-01  * (-10);
  6069     else if ( newt <= 0.34  + 0.004 )
  6071         return  -2.08029000e-01  * (-10);
  6073     else if ( newt <= 0.344 + 0.004 )
  6075         return  -2.37931000e-01  * (-10);
  6077     else if ( newt <= 0.348 + 0.004 )
  6079         return  -2.68601000e-01  * (-10);
  6081     else if ( newt <= 0.352 + 0.004 )
  6083         return  -2.71666000e-01  * (-10);
  6085     else if ( newt <= 0.356 + 0.004 )
  6087         return  -2.74735000e-01  * (-10);
  6089     else if ( newt <= 0.36  + 0.004 )
  6091         return  -2.70136000e-01  * (-10);
  6093     else if ( newt <= 0.364 + 0.004 )
  6095         return  -2.65532000e-01  * (-10);
  6097     else if ( newt <= 0.368 + 0.004 )
  6099         return  -2.59398000e-01  * (-10);
  6101     else if ( newt <= 0.372 + 0.004 )
  6103         return  -2.59398000e-01  * (-10);
  6105     else if ( newt <= 0.376 + 0.004 )
  6107         return  -2.60165000e-01  * (-10);
  6109     else if ( newt <= 0.38  + 0.004 )
  6111         return  -2.60932000e-01  * (-10);
  6113     else if ( newt <= 0.384 + 0.004 )
  6115         return  -2.60165000e-01  * (-10);
  6117     else if ( newt <= 0.388 + 0.004 )
  6119         return  -2.57101000e-01  * (-10);
  6121     else if ( newt <= 0.392 + 0.004 )
  6123         return  -2.58635000e-01  * (-10);
  6125     else if ( newt <= 0.396 + 0.004 )
  6127         return  -2.56333000e-01  * (-10);
  6129     else if ( newt <= 0.4   + 0.004 )
  6131         return  -2.49432000e-01  * (-10);
  6133     else if ( newt <= 0.404 + 0.004 )
  6135         return  -2.43298000e-01  * (-10);
  6137     else if ( newt <= 0.408 + 0.004 )
  6139         return  -2.23361000e-01  * (-10);
  6141     else if ( newt <= 0.412 + 0.004 )
  6143         return  -2.09563000e-01  * (-10);
  6145     else if ( newt <= 0.416 + 0.004 )
  6147         return  -1.71992000e-01  * (-10);
  6149     else if ( newt <= 0.42  + 0.004 )
  6151         return  -1.22151700e-01  * (-10);
  6153     else if ( newt <= 0.424 + 0.004 )
  6155         return  -9.14815000e-02  * (-10);
  6157     else if ( newt <= 0.428 + 0.004 )
  6159         return  -6.08113000e-02  * (-10);
  6161     else if ( newt <= 0.432 + 0.004 )
  6163         return  -5.69779000e-02  * (-10);
  6165     else if ( newt <= 0.436 + 0.004 )
  6167         return  -5.16103000e-02  * (-10);
  6169     else if ( newt <= 0.44  + 0.004 )
  6171         return  -3.93414000e-02  * (-10);
  6173     else if ( newt <= 0.444 + 0.004 )
  6175         return  -1.78748000e-02  * (-10);
  6177     else if ( newt <= 0.448 + 0.004 )
  6179         return  2.42960000e-02   * (-10);
  6181     else if ( newt <= 0.452 + 0.004 )
  6183         return  3.80990000e-02   * (-10);
  6185     else if ( newt <= 0.456 + 0.004 )
  6187         return  4.26980000e-02   * (-10);
  6189     else if ( newt <= 0.46  + 0.004 )
  6191         return  4.26980000e-02   * (-10);
  6193     else if ( newt <= 0.464 + 0.004 )
  6195         return  4.50000000e-02   * (-10);
  6197     else if ( newt <= 0.468 + 0.004 )
  6199         return  4.57670000e-02   * (-10);
  6201     else if ( newt <= 0.472 + 0.004 )
  6203         return  4.73020000e-02   * (-10);
  6205     else if ( newt <= 0.476 + 0.004 )
  6207         return  4.88320000e-02   * (-10);
  6209     else if ( newt <= 0.48  + 0.004 )
  6211         return  5.26690000e-02   * (-10);
  6213     else if ( newt <= 0.484 + 0.004 )
  6215         return  5.49660000e-02   * (-10);
  6217     else if ( newt <= 0.488 + 0.004 )
  6219         return  5.49660000e-02   * (-10);
  6221     else if ( newt <= 0.492 + 0.004 )
  6223         return  5.72680000e-02   * (-10);
  6225     else if ( newt <= 0.496 + 0.004 )
  6227         return  5.95700000e-02   * (-10);
  6229     else if ( newt <= 0.5   + 0.004 )
  6231         return  6.11000000e-02   * (-10);
  6233     else if ( newt <= 0.504 + 0.004 )
  6235         return  6.18670000e-02   * (-10);
  6237     else if ( newt <= 0.508 + 0.004 )
  6239         return  6.49370000e-02   * (-10);
  6241     else if ( newt <= 0.512 + 0.004 )
  6243         return  6.87680000e-02   * (-10);
  6245     else if ( newt <= 0.516 + 0.004 )
  6247         return  6.64710000e-02   * (-10);
  6249     else if ( newt <= 0.52  + 0.004 )
  6251         return  6.41690000e-02   * (-10);
  6253     else if ( newt <= 0.524 + 0.004 )
  6255         return  6.57040000e-02   * (-10);
  6257     else if ( newt <= 0.528 + 0.004 )
  6259         return  6.87680000e-02   * (-10);
  6261     else if ( newt <= 0.532 + 0.004 )
  6263         return  7.18380000e-02   * (-10);
  6265     else if ( newt <= 0.536 + 0.004 )
  6267         return  7.64370000e-02   * (-10);
  6269     else if ( newt <= 0.54  + 0.004 )
  6271         return  7.64370000e-02   * (-10);
  6273     else if ( newt <= 0.544 + 0.004 )
  6275         return  7.41350000e-02   * (-10);
  6277     else if ( newt <= 0.548 + 0.004 )
  6279         return  7.49020000e-02   * (-10);
  6281     else if ( newt <= 0.552 + 0.004 )
  6283         return  7.56700000e-02   * (-10);
  6285     else if ( newt <= 0.556 + 0.004 )
  6287         return  7.33680000e-02   * (-10);
  6289     else if ( newt <= 0.56  + 0.004 )
  6291         return  6.87680000e-02   * (-10);
  6293     else if ( newt <= 0.564 + 0.004 )
  6295         return  7.03030000e-02   * (-10);
  6297     else if ( newt <= 0.568 + 0.004 )
  6299         return  6.72340000e-02   * (-10);
  6301     else if ( newt <= 0.572 + 0.004 )
  6303         return  6.41690000e-02   * (-10);
  6305     else if ( newt <= 0.576 + 0.004 )
  6307         return  6.11000000e-02   * (-10);
  6309     else if ( newt <= 0.58  + 0.004 )
  6311         return  6.11000000e-02   * (-10);
  6313     else if ( newt <= 0.584 + 0.004 )
  6315         return  6.41690000e-02   * (-10);
  6317     else if ( newt <= 0.588 + 0.004 )
  6319         return  7.18380000e-02   * (-10);
  6321     else if ( newt <= 0.592 + 0.004 )
  6323         return  7.95020000e-02   * (-10);
  6325     else if ( newt <= 0.596 + 0.004 )
  6327         return  8.56360000e-02   * (-10);
  6329     else if ( newt <= 0.6   + 0.004 )
  6331         return  8.25710000e-02   * (-10);
  6333     else if ( newt <= 0.604 + 0.004 )
  6335         return  7.64370000e-02   * (-10);
  6337     else if ( newt <= 0.608 + 0.004 )
  6339         return  7.10710000e-02   * (-10);
  6341     else if ( newt <= 0.612 + 0.004 )
  6343         return  7.33680000e-02   * (-10);
  6345     else if ( newt <= 0.616 + 0.004 )
  6347         return  7.18380000e-02   * (-10);
  6349     else if ( newt <= 0.62  + 0.004 )
  6351         return  7.10710000e-02   * (-10);
  6353     else if ( newt <= 0.624 + 0.004 )
  6355         return  6.80010000e-02   * (-10);
  6357     else if ( newt <= 0.628 + 0.004 )
  6359         return  6.95360000e-02   * (-10);
  6361     else if ( newt <= 0.632 + 0.004 )
  6363         return  7.03030000e-02   * (-10);
  6365     else if ( newt <= 0.636 + 0.004 )
  6367         return  6.72340000e-02   * (-10);
  6369     else if ( newt <= 0.64  + 0.004 )
  6371         return  6.18670000e-02   * (-10);
  6373     else if ( newt <= 0.644 + 0.004 )
  6375         return  5.65010000e-02   * (-10);
  6377     else if ( newt <= 0.648 + 0.004 )
  6379         return  5.65010000e-02   * (-10);
  6381     else if ( newt <= 0.652 + 0.004 )
  6383         return  5.49660000e-02   * (-10);
  6385     else if ( newt <= 0.656 + 0.004 )
  6387         return  5.65010000e-02   * (-10);
  6389     else if ( newt <= 0.66  + 0.004 )
  6391         return  5.88030000e-02   * (-10);
  6393     else if ( newt <= 0.664 + 0.004 )
  6395         return  5.95700000e-02   * (-10);
  6397     else if ( newt <= 0.668 + 0.004 )
  6399         return  6.03370000e-02   * (-10);
  6401     else if ( newt <= 0.672 + 0.004 )
  6403         return  6.03370000e-02   * (-10);
  6405     else if ( newt <= 0.676 + 0.004 )
  6407         return  5.95700000e-02   * (-10);
  6409     else if ( newt <= 0.68  + 0.004 )
  6411         return  6.11000000e-02   * (-10);
  6413     else if ( newt <= 0.684 + 0.004 )
  6415         return  6.11000000e-02   * (-10);
  6417     else if ( newt <= 0.688 + 0.004 )
  6419         return  5.95700000e-02   * (-10);
  6421     else if ( newt <= 0.692 + 0.004 )
  6423         return  5.88030000e-02   * (-10);
  6425     else if ( newt <= 0.696 + 0.004 )
  6427         return  5.80350000e-02   * (-10);
  6429     else if ( newt <= 0.7   + 0.004 )
  6431         return  6.03370000e-02   * (-10);
  6433     else if ( newt <= 0.704 + 0.004 )
  6435         return  5.95700000e-02   * (-10);
  6437     else if ( newt <= 0.708 + 0.004 )
  6439         return  5.80350000e-02   * (-10);
  6441     else if ( newt <= 0.712 + 0.004 )
  6443         return  5.42030000e-02   * (-10);
  6445     else if ( newt <= 0.716 + 0.004 )
  6447         return  5.19010000e-02   * (-10);
  6449     else if ( newt <= 0.72  + 0.004 )
  6451         return  4.65350000e-02   * (-10);
  6453     else if ( newt <= 0.724 + 0.004 )
  6455         return  4.50000000e-02   * (-10);
  6457     else if ( newt <= 0.728 + 0.004 )
  6459         return  4.19350000e-02   * (-10);
  6461     else if ( newt <= 0.732 + 0.004 )
  6463         return  4.26980000e-02   * (-10);
  6465     else if ( newt <= 0.736 + 0.004 )
  6467         return  4.42330000e-02   * (-10);
  6469     else if ( newt <= 0.74  + 0.004 )
  6471         return  4.80690000e-02   * (-10);
  6473     else if ( newt <= 0.744 + 0.004 )
  6475         return  4.42330000e-02   * (-10);
  6477     else if ( newt <= 0.748 + 0.004 )
  6479         return  4.11680000e-02   * (-10);
  6481     else if ( newt <= 0.752 + 0.004 )
  6483         return  3.65640000e-02   * (-10);
  6485     else if ( newt <= 0.756 + 0.004 )
  6487         return  3.58010000e-02   * (-10);
  6489     else if ( newt <= 0.76  + 0.004 )
  6491         return  3.42670000e-02   * (-10);
  6493     else if ( newt <= 0.764 + 0.004 )
  6495         return  3.27320000e-02   * (-10);
  6497     else if ( newt <= 0.768 + 0.004 )
  6499         return  3.11970000e-02   * (-10);
  6501     else if ( newt <= 0.772 + 0.004 )
  6503         return  3.04300000e-02   * (-10);
  6505     else if ( newt <= 0.776 + 0.004 )
  6507         return  2.96670000e-02   * (-10);
  6509     else if ( newt <= 0.78  + 0.004 )
  6511         return  3.04300000e-02   * (-10);
  6513     else if ( newt <= 0.784 + 0.004 )
  6515         return  2.96670000e-02   * (-10);
  6517     else if ( newt <= 0.788 + 0.004 )
  6519         return  2.50630000e-02   * (-10);
  6521     else if ( newt <= 0.792 + 0.004 )
  6523         return  2.27660000e-02   * (-10);
  6525     else if ( newt <= 0.796 + 0.004 )
  6527         return  2.12310000e-02   * (-10);
  6529     else if ( newt <= 0.8   + 0.004 )
  6531         return  2.23800000e-02   * (-10);
  6539         t = ( ( (
int) floor (t * 1000.0) ) % 800) / 1000.0;
  6542     Real dt = ( (
int) floor (t * 1000.0) ) % 4 / 1000.0;
  6543     Real ti = floor (t * 1000.0) / 1000.0 - dt;
  6544     Real tii = ti + 0.004;
  6545     return (popliteal (tii) - popliteal (ti) ) / (0.004) * (t - (ti) ) + popliteal (ti);
  6552         t = ( ( (
int) floor (t * 1000.0) ) % 800) / 1000.0;
  6555     Real dt = ( (
int) floor (t * 1000.0) ) % 4 / 1000.0;
  6556     Real ti = floor (t * 1000.0) / 1000.0 - dt;
  6557     Real tii = ti + 0.004;
  6558     return (pont_dist (tii) - pont_dist (ti) ) / (0.004) * (t - (ti) ) + pont_dist (ti);
  6564     Real newt = ( (
Real) ( ( (
int) (t * 1000) ) % 800) ) / 1000;
  6565     if ( newt <=          0 + 0.001 )
  6569     else if ( newt <=     0.0010 + 0.001 )
  6571         return    -0.2595   * (1.0e4);
  6573     else if ( newt <=     0.0020 + 0.001 )
  6575         return    -0.1949   * (1.0e4);
  6577     else if ( newt <=     0.0030 + 0.001 )
  6579         return    -0.1639   * (1.0e4);
  6581     else if ( newt <=     0.0040 + 0.001 )
  6583         return    -0.1380   * (1.0e4);
  6585     else if ( newt <=     0.0050 + 0.001 )
  6587         return    -0.1340   * (1.0e4);
  6589     else if ( newt <=     0.0060 + 0.001 )
  6591         return    -0.1268   * (1.0e4);
  6593     else if ( newt <=     0.0070 + 0.001 )
  6595         return    -0.1215   * (1.0e4);
  6597     else if ( newt <=     0.0080 + 0.001 )
  6599         return    -0.1210   * (1.0e4);
  6601     else if ( newt <=     0.0090 + 0.001 )
  6603         return    -0.1411   * (1.0e4);
  6605     else if ( newt <=     0.0100 + 0.001 )
  6607         return    -0.1653   * (1.0e4);
  6609     else if ( newt <=     0.0110 + 0.001 )
  6611         return    -0.1953   * (1.0e4);
  6613     else if ( newt <=     0.0120 + 0.001 )
  6615         return    -0.2301   * (1.0e4);
  6617     else if ( newt <=     0.0130 + 0.001 )
  6619         return    -0.2657   * (1.0e4);
  6621     else if ( newt <=     0.0140 + 0.001 )
  6623         return    -0.3023   * (1.0e4);
  6625     else if ( newt <=     0.0150 + 0.001 )
  6627         return    -0.3377   * (1.0e4);
  6629     else if ( newt <=     0.0160 + 0.001 )
  6631         return    -0.3705   * (1.0e4);
  6633     else if ( newt <=     0.0170 + 0.001 )
  6635         return    -0.3932   * (1.0e4);
  6637     else if ( newt <=     0.0180 + 0.001 )
  6639         return    -0.4152   * (1.0e4);
  6641     else if ( newt <=     0.0190 + 0.001 )
  6643         return    -0.4363   * (1.0e4);
  6645     else if ( newt <=     0.0200 + 0.001 )
  6647         return    -0.4580   * (1.0e4);
  6649     else if ( newt <=     0.0210 + 0.001 )
  6651         return    -0.4829   * (1.0e4);
  6653     else if ( newt <=     0.0220 + 0.001 )
  6655         return    -0.5100   * (1.0e4);
  6657     else if ( newt <=     0.0230 + 0.001 )
  6659         return    -0.5399   * (1.0e4);
  6661     else if ( newt <=     0.0240 + 0.001 )
  6663         return    -0.5725   * (1.0e4);
  6665     else if ( newt <=     0.0250 + 0.001 )
  6667         return    -0.6086   * (1.0e4);
  6669     else if ( newt <=     0.0260 + 0.001 )
  6671         return    -0.6456   * (1.0e4);
  6673     else if ( newt <=     0.0270 + 0.001 )
  6675         return    -0.6830   * (1.0e4);
  6677     else if ( newt <=     0.0280 + 0.001 )
  6679         return    -0.7199   * (1.0e4);
  6681     else if ( newt <=     0.0290 + 0.001 )
  6683         return    -0.7515   * (1.0e4);
  6685     else if ( newt <=     0.0300 + 0.001 )
  6687         return    -0.7830   * (1.0e4);
  6689     else if ( newt <=     0.0310 + 0.001 )
  6691         return    -0.8136   * (1.0e4);
  6693     else if ( newt <=     0.0320 + 0.001 )
  6695         return    -0.8436   * (1.0e4);
  6697     else if ( newt <=     0.0330 + 0.001 )
  6699         return    -0.8720   * (1.0e4);
  6701     else if ( newt <=     0.0340 + 0.001 )
  6703         return    -0.9012   * (1.0e4);
  6705     else if ( newt <=     0.0350 + 0.001 )
  6707         return    -0.9312   * (1.0e4);
  6709     else if ( newt <=     0.0360 + 0.001 )
  6711         return    -0.9621   * (1.0e4);
  6713     else if ( newt <=     0.0370 + 0.001 )
  6715         return    -0.9926   * (1.0e4);
  6717     else if ( newt <=     0.0380 + 0.001 )
  6719         return    -1.0240   * (1.0e4);
  6721     else if ( newt <=     0.0390 + 0.001 )
  6723         return    -1.0559   * (1.0e4);
  6725     else if ( newt <=     0.0400 + 0.001 )
  6727         return    -1.2131   * (1.0e4);
  6729     else if ( newt <=     0.0410 + 0.001 )
  6731         return    -1.1907   * (1.0e4);
  6733     else if ( newt <=     0.0420 + 0.001 )
  6735         return    -1.1903   * (1.0e4);
  6737     else if ( newt <=     0.0430 + 0.001 )
  6739         return    -1.1610   * (1.0e4);
  6741     else if ( newt <=     0.0440 + 0.001 )
  6743         return    -1.2035   * (1.0e4);
  6745     else if ( newt <=     0.0450 + 0.001 )
  6747         return    -1.2322   * (1.0e4);
  6749     else if ( newt <=     0.0460 + 0.001 )
  6751         return    -1.2613   * (1.0e4);
  6753     else if ( newt <=     0.0470 + 0.001 )
  6755         return    -1.2932   * (1.0e4);
  6757     else if ( newt <=     0.0480 + 0.001 )
  6759         return    -1.3292   * (1.0e4);
  6761     else if ( newt <=     0.0490 + 0.001 )
  6763         return    -1.3769   * (1.0e4);
  6765     else if ( newt <=     0.0500 + 0.001 )
  6767         return    -1.4237   * (1.0e4);
  6769     else if ( newt <=     0.0510 + 0.001 )
  6771         return    -1.4686   * (1.0e4);
  6773     else if ( newt <=     0.0520 + 0.001 )
  6775         return    -1.5097   * (1.0e4);
  6777     else if ( newt <=     0.0530 + 0.001 )
  6779         return    -1.5436   * (1.0e4);
  6781     else if ( newt <=     0.0540 + 0.001 )
  6783         return    -1.5738   * (1.0e4);
  6785     else if ( newt <=     0.0550 + 0.001 )
  6787         return    -1.6007   * (1.0e4);
  6789     else if ( newt <=     0.0560 + 0.001 )
  6791         return    -1.7456   * (1.0e4);
  6793     else if ( newt <=     0.0570 + 0.001 )
  6795         return    -1.7042   * (1.0e4);
  6797     else if ( newt <=     0.0580 + 0.001 )
  6799         return    -1.6879   * (1.0e4);
  6801     else if ( newt <=     0.0590 + 0.001 )
  6803         return    -1.6474   * (1.0e4);
  6805     else if ( newt <=     0.0600 + 0.001 )
  6807         return    -1.6796   * (1.0e4);
  6809     else if ( newt <=     0.0610 + 0.001 )
  6811         return    -1.6954   * (1.0e4);
  6813     else if ( newt <=     0.0620 + 0.001 )
  6815         return    -1.7171   * (1.0e4);
  6817     else if ( newt <=     0.0630 + 0.001 )
  6819         return    -1.7454   * (1.0e4);
  6821     else if ( newt <=     0.0640 + 0.001 )
  6823         return    -1.7803   * (1.0e4);
  6825     else if ( newt <=     0.0650 + 0.001 )
  6827         return    -1.8083   * (1.0e4);
  6829     else if ( newt <=     0.0660 + 0.001 )
  6831         return    -1.8407   * (1.0e4);
  6833     else if ( newt <=     0.0670 + 0.001 )
  6835         return    -1.8726   * (1.0e4);
  6837     else if ( newt <=     0.0680 + 0.001 )
  6839         return    -1.9017   * (1.0e4);
  6841     else if ( newt <=     0.0690 + 0.001 )
  6843         return    -1.9309   * (1.0e4);
  6845     else if ( newt <=     0.0700 + 0.001 )
  6847         return    -1.9547   * (1.0e4);
  6849     else if ( newt <=     0.0710 + 0.001 )
  6851         return    -1.9744   * (1.0e4);
  6853     else if ( newt <=     0.0720 + 0.001 )
  6855         return    -1.9910   * (1.0e4);
  6857     else if ( newt <=     0.0730 + 0.001 )
  6859         return    -2.0102   * (1.0e4);
  6861     else if ( newt <=     0.0740 + 0.001 )
  6863         return    -2.0277   * (1.0e4);
  6865     else if ( newt <=     0.0750 + 0.001 )
  6867         return    -2.0449   * (1.0e4);
  6869     else if ( newt <=     0.0760 + 0.001 )
  6871         return    -2.0622   * (1.0e4);
  6873     else if ( newt <=     0.0770 + 0.001 )
  6875         return    -2.0840   * (1.0e4);
  6877     else if ( newt <=     0.0780 + 0.001 )
  6879         return    -2.1044   * (1.0e4);
  6881     else if ( newt <=     0.0790 + 0.001 )
  6883         return    -2.1239   * (1.0e4);
  6885     else if ( newt <=     0.0800 + 0.001 )
  6887         return    -2.1421   * (1.0e4);
  6889     else if ( newt <=     0.0810 + 0.001 )
  6891         return    -2.1545   * (1.0e4);
  6893     else if ( newt <=     0.0820 + 0.001 )
  6895         return    -2.16662  * (1.0e4);
  6897     else if ( newt <=     0.0830 + 0.001 )
  6899         return    -2.17782  * (1.0e4);
  6901     else if ( newt <=     0.0840 + 0.001 )
  6903         return    -2.18851  * (1.0e4);
  6905     else if ( newt <=     0.0850 + 0.001 )
  6907         return    -2.19451  * (1.0e4);
  6909     else if ( newt <=     0.0860 + 0.001 )
  6911         return    -2.20121  * (1.0e4);
  6913     else if ( newt <=     0.0870 + 0.001 )
  6915         return    -2.20831  * (1.0e4);
  6917     else if ( newt <=     0.0880 + 0.001 )
  6919         return    -2.21571  * (1.0e4);
  6921     else if ( newt <=     0.0890 + 0.001 )
  6923         return    -2.22621  * (1.0e4);
  6925     else if ( newt <=     0.0900 + 0.001 )
  6927         return    -2.23631  * (1.0e4);
  6929     else if ( newt <=     0.0910 + 0.001 )
  6931         return    -2.24611  * (1.0e4);
  6933     else if ( newt <=     0.0920 + 0.001 )
  6935         return    -2.25551  * (1.0e4);
  6937     else if ( newt <=     0.0930 + 0.001 )
  6939         return    -2.25891  * (1.0e4);
  6941     else if ( newt <=     0.0940 + 0.001 )
  6943         return    -2.26351  * (1.0e4);
  6945     else if ( newt <=     0.0950 + 0.001 )
  6947         return    -2.26861  * (1.0e4);
  6949     else if ( newt <=     0.0960 + 0.001 )
  6951         return    -2.27411  * (1.0e4);
  6953     else if ( newt <=     0.0970 + 0.001 )
  6955         return    -2.27451  * (1.0e4);
  6957     else if ( newt <=     0.0980 + 0.001 )
  6959         return    -2.27681  * (1.0e4);
  6961     else if ( newt <=     0.0990 + 0.001 )
  6963         return    -2.28041  * (1.0e4);
  6965     else if ( newt <=     0.1000 + 0.001 )
  6967         return    -2.28481  * (1.0e4);
  6969     else if ( newt <=     0.1010 + 0.001 )
  6971         return    -2.29021  * (1.0e4);
  6973     else if ( newt <=     0.1020 + 0.001 )
  6975         return    -2.29621  * (1.0e4);
  6977     else if ( newt <=     0.1030 + 0.001 )
  6979         return    -2.30261  * (1.0e4);
  6981     else if ( newt <=     0.1040 + 0.001 )
  6983         return    -2.30941  * (1.0e4);
  6985     else if ( newt <=     0.1050 + 0.001 )
  6987         return    -2.32491  * (1.0e4);
  6989     else if ( newt <=     0.1060 + 0.001 )
  6991         return    -2.33822  * (1.0e4);
  6993     else if ( newt <=     0.1070 + 0.001 )
  6995         return    -2.35052  * (1.0e4);
  6997     else if ( newt <=     0.1080 + 0.001 )
  6999         return    -2.36182  * (1.0e4);
  7001     else if ( newt <=     0.1090 + 0.001 )
  7003         return    -2.37172  * (1.0e4);
  7005     else if ( newt <=     0.1100 + 0.001 )
  7007         return    -2.38012  * (1.0e4);
  7009     else if ( newt <=     0.1110 + 0.001 )
  7011         return    -2.38712  * (1.0e4);
  7013     else if ( newt <=     0.1120 + 0.001 )
  7015         return    -2.3926   * (1.0e4);
  7017     else if ( newt <=     0.1130 + 0.001 )
  7019         return    -2.3984   * (1.0e4);
  7021     else if ( newt <=     0.1140 + 0.001 )
  7023         return    -2.4031   * (1.0e4);
  7025     else if ( newt <=     0.1150 + 0.001 )
  7027         return    -2.4070   * (1.0e4);
  7029     else if ( newt <=     0.1160 + 0.001 )
  7031         return    -2.4106   * (1.0e4);
  7033     else if ( newt <=     0.1170 + 0.001 )
  7035         return    -2.3952   * (1.0e4);
  7037     else if ( newt <=     0.1180 + 0.001 )
  7039         return    -2.3841   * (1.0e4);
  7041     else if ( newt <=     0.1190 + 0.001 )
  7043         return    -2.3751   * (1.0e4);
  7045     else if ( newt <=     0.1200 + 0.001 )
  7047         return    -2.3675   * (1.0e4);
  7049     else if ( newt <=     0.1210 + 0.001 )
  7051         return    -2.3569   * (1.0e4);
  7053     else if ( newt <=     0.1220 + 0.001 )
  7055         return    -2.3482   * (1.0e4);
  7057     else if ( newt <=     0.1230 + 0.001 )
  7059         return    -2.3403   * (1.0e4);
  7061     else if ( newt <=     0.1240 + 0.001 )
  7063         return    -2.3329   * (1.0e4);
  7065     else if ( newt <=     0.1250 + 0.001 )
  7067         return    -2.3284   * (1.0e4);
  7069     else if ( newt <=     0.1260 + 0.001 )
  7071         return    -2.3228   * (1.0e4);
  7073     else if ( newt <=     0.1270 + 0.001 )
  7075         return    -2.3163   * (1.0e4);
  7077     else if ( newt <=     0.1280 + 0.001 )
  7079         return    -2.3086   * (1.0e4);
  7081     else if ( newt <=     0.1290 + 0.001 )
  7083         return    -2.3155   * (1.0e4);
  7085     else if ( newt <=     0.1300 + 0.001 )
  7087         return    -2.3172   * (1.0e4);
  7089     else if ( newt <=     0.1310 + 0.001 )
  7091         return    -2.3160   * (1.0e4);
  7093     else if ( newt <=     0.1320 + 0.001 )
  7095         return    -2.3125   * (1.0e4);
  7097     else if ( newt <=     0.1330 + 0.001 )
  7099         return    -2.3087   * (1.0e4);
  7101     else if ( newt <=     0.1340 + 0.001 )
  7103         return    -2.3033   * (1.0e4);
  7105     else if ( newt <=     0.1350 + 0.001 )
  7107         return    -2.2973   * (1.0e4);
  7109     else if ( newt <=     0.1360 + 0.001 )
  7111         return    -2.2912   * (1.0e4);
  7113     else if ( newt <=     0.1370 + 0.001 )
  7115         return    -2.2752   * (1.0e4);
  7117     else if ( newt <=     0.1380 + 0.001 )
  7119         return    -2.2620   * (1.0e4);
  7121     else if ( newt <=     0.1390 + 0.001 )
  7123         return    -2.2505   * (1.0e4);
  7125     else if ( newt <=     0.1400 + 0.001 )
  7127         return    -2.2405   * (1.0e4);
  7129     else if ( newt <=     0.1410 + 0.001 )
  7131         return    -2.2216   * (1.0e4);
  7133     else if ( newt <=     0.1420 + 0.001 )
  7135         return    -2.2060   * (1.0e4);
  7137     else if ( newt <=     0.1430 + 0.001 )
  7139         return    -2.1920   * (1.0e4);
  7141     else if ( newt <=     0.1440 + 0.001 )
  7143         return    -2.1790   * (1.0e4);
  7145     else if ( newt <=     0.1450 + 0.001 )
  7147         return    -2.1681   * (1.0e4);
  7149     else if ( newt <=     0.1460 + 0.001 )
  7151         return    -2.1569   * (1.0e4);
  7153     else if ( newt <=     0.1470 + 0.001 )
  7155         return    -2.1452   * (1.0e4);
  7157     else if ( newt <=     0.1480 + 0.001 )
  7159         return    -2.1326   * (1.0e4);
  7161     else if ( newt <=     0.1490 + 0.001 )
  7163         return    -2.1381   * (1.0e4);
  7165     else if ( newt <=     0.1500 + 0.001 )
  7167         return    -2.1382   * (1.0e4);
  7169     else if ( newt <=     0.1510 + 0.001 )
  7171         return    -2.1357   * (1.0e4);
  7173     else if ( newt <=     0.1520 + 0.001 )
  7175         return    -2.1310   * (1.0e4);
  7177     else if ( newt <=     0.1530 + 0.001 )
  7179         return    -2.1493   * (1.0e4);
  7181     else if ( newt <=     0.1540 + 0.001 )
  7183         return    -2.1608   * (1.0e4);
  7185     else if ( newt <=     0.1550 + 0.001 )
  7187         return    -2.1695   * (1.0e4);
  7189     else if ( newt <=     0.1560 + 0.001 )
  7191         return    -2.1763   * (1.0e4);
  7193     else if ( newt <=     0.1570 + 0.001 )
  7195         return    -2.2154   * (1.0e4);
  7197     else if ( newt <=     0.1580 + 0.001 )
  7199         return    -2.2463   * (1.0e4);
  7201     else if ( newt <=     0.1590 + 0.001 )
  7203         return    -2.2742   * (1.0e4);
  7205     else if ( newt <=     0.1600 + 0.001 )
  7207         return    -2.3002   * (1.0e4);
  7209     else if ( newt <=     0.1610 + 0.001 )
  7211         return    -2.3619   * (1.0e4);
  7213     else if ( newt <=     0.1620 + 0.001 )
  7215         return    -2.4153   * (1.0e4);
  7217     else if ( newt <=     0.1630 + 0.001 )
  7219         return    -2.4664   * (1.0e4);
  7221     else if ( newt <=     0.1640 + 0.001 )
  7223         return    -2.5167   * (1.0e4);
  7225     else if ( newt <=     0.1650 + 0.001 )
  7227         return    -2.5589   * (1.0e4);
  7229     else if ( newt <=     0.1660 + 0.001 )
  7231         return    -2.6054   * (1.0e4);
  7233     else if ( newt <=     0.1670 + 0.001 )
  7235         return    -2.6563   * (1.0e4);
  7237     else if ( newt <=     0.1680 + 0.001 )
  7239         return    -2.7130   * (1.0e4);
  7241     else if ( newt <=     0.1690 + 0.001 )
  7243         return    -2.7059   * (1.0e4);
  7245     else if ( newt <=     0.1700 + 0.001 )
  7247         return    -2.7218   * (1.0e4);
  7249     else if ( newt <=     0.1710 + 0.001 )
  7251         return    -2.7521   * (1.0e4);
  7253     else if ( newt <=     0.1720 + 0.001 )
  7255         return    -2.7960   * (1.0e4);
  7257     else if ( newt <=     0.1730 + 0.001 )
  7259         return    -2.7344   * (1.0e4);
  7261     else if ( newt <=     0.1740 + 0.001 )
  7263         return    -2.7099   * (1.0e4);
  7265     else if ( newt <=     0.1750 + 0.001 )
  7267         return    -2.7057   * (1.0e4);
  7269     else if ( newt <=     0.1760 + 0.001 )
  7271         return    -2.7178   * (1.0e4);
  7273     else if ( newt <=     0.1770 + 0.001 )
  7275         return    -2.7231   * (1.0e4);
  7277     else if ( newt <=     0.1780 + 0.001 )
  7279         return    -2.7414   * (1.0e4);
  7281     else if ( newt <=     0.1790 + 0.001 )
  7283         return    -2.7661   * (1.0e4);
  7285     else if ( newt <=     0.1800 + 0.001 )
  7287         return    -2.7930   * (1.0e4);
  7289     else if ( newt <=     0.1810 + 0.001 )
  7291         return    -2.8229   * (1.0e4);
  7293     else if ( newt <=     0.1820 + 0.001 )
  7295         return    -2.8472   * (1.0e4);
  7297     else if ( newt <=     0.1830 + 0.001 )
  7299         return    -2.8642   * (1.0e4);
  7301     else if ( newt <=     0.1840 + 0.001 )
  7303         return    -2.8726   * (1.0e4);
  7305     else if ( newt <=     0.1850 + 0.001 )
  7307         return    -2.9007   * (1.0e4);
  7309     else if ( newt <=     0.1860 + 0.001 )
  7311         return    -2.9128   * (1.0e4);
  7313     else if ( newt <=     0.1870 + 0.001 )
  7315         return    -2.9136   * (1.0e4);
  7317     else if ( newt <=     0.1880 + 0.001 )
  7319         return    -2.9048   * (1.0e4);
  7321     else if ( newt <=     0.1890 + 0.001 )
  7323         return    -2.9170   * (1.0e4);
  7325     else if ( newt <=     0.1900 + 0.001 )
  7327         return    -2.9159   * (1.0e4);
  7329     else if ( newt <=     0.1910 + 0.001 )
  7331         return    -2.9069   * (1.0e4);
  7333     else if ( newt <=     0.1920 + 0.001 )
  7335         return    -2.8916   * (1.0e4);
  7337     else if ( newt <=     0.1930 + 0.001 )
  7339         return    -2.8918   * (1.0e4);
  7341     else if ( newt <=     0.1940 + 0.001 )
  7343         return    -2.8832   * (1.0e4);
  7345     else if ( newt <=     0.1950 + 0.001 )
  7347         return    -2.8695   * (1.0e4);
  7349     else if ( newt <=     0.1960 + 0.001 )
  7351         return    -2.8520   * (1.0e4);
  7353     else if ( newt <=     0.1970 + 0.001 )
  7355         return    -2.8502   * (1.0e4);
  7357     else if ( newt <=     0.1980 + 0.001 )
  7359         return    -2.8418   * (1.0e4);
  7361     else if ( newt <=     0.1990 + 0.001 )
  7363         return    -2.8301   * (1.0e4);
  7365     else if ( newt <=     0.2000 + 0.001 )
  7367         return    -2.8160   * (1.0e4);
  7369     else if ( newt <=     0.2010 + 0.001 )
  7371         return    -2.8076   * (1.0e4);
  7373     else if ( newt <=     0.2020 + 0.001 )
  7375         return    -2.7965   * (1.0e4);
  7377     else if ( newt <=     0.2030 + 0.001 )
  7379         return    -2.7843   * (1.0e4);
  7381     else if ( newt <=     0.2040 + 0.001 )
  7383         return    -2.7715   * (1.0e4);
  7385     else if ( newt <=     0.2050 + 0.001 )
  7387         return    -2.7587   * (1.0e4);
  7389     else if ( newt <=     0.2060 + 0.001 )
  7391         return    -2.7465   * (1.0e4);
  7393     else if ( newt <=     0.2070 + 0.001 )
  7395         return    -2.7351   * (1.0e4);
  7397     else if ( newt <=     0.2080 + 0.001 )
  7399         return    -2.7250   * (1.0e4);
  7401     else if ( newt <=     0.2090 + 0.001 )
  7403         return    -2.7147   * (1.0e4);
  7405     else if ( newt <=     0.2100 + 0.001 )
  7407         return    -2.7061   * (1.0e4);
  7409     else if ( newt <=     0.2110 + 0.001 )
  7411         return    -2.6990   * (1.0e4);
  7413     else if ( newt <=     0.2120 + 0.001 )
  7415         return    -2.6933   * (1.0e4);
  7417     else if ( newt <=     0.2130 + 0.001 )
  7419         return    -2.6668   * (1.0e4);
  7421     else if ( newt <=     0.2140 + 0.001 )
  7423         return    -2.6461   * (1.0e4);
  7425     else if ( newt <=     0.2150 + 0.001 )
  7427         return    -2.6279   * (1.0e4);
  7429     else if ( newt <=     0.2160 + 0.001 )
  7431         return    -2.6115   * (1.0e4);
  7433     else if ( newt <=     0.2170 + 0.001 )
  7435         return    -2.5895   * (1.0e4);
  7437     else if ( newt <=     0.2180 + 0.001 )
  7439         return    -2.5705   * (1.0e4);
  7441     else if ( newt <=     0.2190 + 0.001 )
  7443         return    -2.5531   * (1.0e4);
  7445     else if ( newt <=     0.2200 + 0.001 )
  7447         return    -2.5369   * (1.0e4);
  7449     else if ( newt <=     0.2210 + 0.001 )
  7451         return    -2.5157   * (1.0e4);
  7453     else if ( newt <=     0.2220 + 0.001 )
  7455         return    -2.4961   * (1.0e4);
  7457     else if ( newt <=     0.2230 + 0.001 )
  7459         return    -2.4770   * (1.0e4);
  7461     else if ( newt <=     0.2240 + 0.001 )
  7463         return    -2.4576   * (1.0e4);
  7465     else if ( newt <=     0.2250 + 0.001 )
  7467         return    -2.4391   * (1.0e4);
  7469     else if ( newt <=     0.2260 + 0.001 )
  7471         return    -2.4195   * (1.0e4);
  7473     else if ( newt <=     0.2270 + 0.001 )
  7475         return    -2.3990   * (1.0e4);
  7477     else if ( newt <=     0.2280 + 0.001 )
  7479         return    -2.3776   * (1.0e4);
  7481     else if ( newt <=     0.2290 + 0.001 )
  7483         return    -2.3639   * (1.0e4);
  7485     else if ( newt <=     0.2300 + 0.001 )
  7487         return    -2.3474   * (1.0e4);
  7489     else if ( newt <=     0.2310 + 0.001 )
  7491         return    -2.3297   * (1.0e4);
  7493     else if ( newt <=     0.2320 + 0.001 )
  7495         return    -2.3107   * (1.0e4);
  7497     else if ( newt <=     0.2330 + 0.001 )
  7499         return    -2.2894   * (1.0e4);
  7501     else if ( newt <=     0.2340 + 0.001 )
  7503         return    -2.2676   * (1.0e4);
  7505     else if ( newt <=     0.2350 + 0.001 )
  7507         return    -2.2453   * (1.0e4);
  7509     else if ( newt <=     0.2360 + 0.001 )
  7511         return    -2.2225   * (1.0e4);
  7513     else if ( newt <=     0.2370 + 0.001 )
  7515         return    -2.1966   * (1.0e4);
  7517     else if ( newt <=     0.2380 + 0.001 )
  7519         return    -2.1713   * (1.0e4);
  7521     else if ( newt <=     0.2390 + 0.001 )
  7523         return    -2.1463   * (1.0e4);
  7525     else if ( newt <=     0.2400 + 0.001 )
  7527         return    -2.1216   * (1.0e4);
  7529     else if ( newt <=     0.2410 + 0.001 )
  7531         return    -2.0871   * (1.0e4);
  7533     else if ( newt <=     0.2420 + 0.001 )
  7535         return    -2.0551   * (1.0e4);
  7537     else if ( newt <=     0.2430 + 0.001 )
  7539         return    -2.0239   * (1.0e4);
  7541     else if ( newt <=     0.2440 + 0.001 )
  7543         return    -1.9930   * (1.0e4);
  7545     else if ( newt <=     0.2450 + 0.001 )
  7547         return    -1.9580   * (1.0e4);
  7549     else if ( newt <=     0.2460 + 0.001 )
  7551         return    -1.9240   * (1.0e4);
  7553     else if ( newt <=     0.2470 + 0.001 )
  7555         return    -1.8902   * (1.0e4);
  7557     else if ( newt <=     0.2480 + 0.001 )
  7559         return    -1.8563   * (1.0e4);
  7561     else if ( newt <=     0.2490 + 0.001 )
  7563         return    -1.7995   * (1.0e4);
  7565     else if ( newt <=     0.2500 + 0.001 )
  7567         return    -1.7480   * (1.0e4);
  7569     else if ( newt <=     0.2510 + 0.001 )
  7571         return    -1.6984   * (1.0e4);
  7573     else if ( newt <=     0.2520 + 0.001 )
  7575         return    -1.6502   * (1.0e4);
  7577     else if ( newt <=     0.2530 + 0.001 )
  7579         return    -1.6212   * (1.0e4);
  7581     else if ( newt <=     0.2540 + 0.001 )
  7583         return    -1.5877   * (1.0e4);
  7585     else if ( newt <=     0.2550 + 0.001 )
  7587         return    -1.5517   * (1.0e4);
  7589     else if ( newt <=     0.2560 + 0.001 )
  7591         return    -1.5129   * (1.0e4);
  7593     else if ( newt <=     0.2570 + 0.001 )
  7595         return    -1.4965   * (1.0e4);
  7597     else if ( newt <=     0.2580 + 0.001 )
  7599         return    -1.4704   * (1.0e4);
  7601     else if ( newt <=     0.2590 + 0.001 )
  7603         return    -1.4382   * (1.0e4);
  7605     else if ( newt <=     0.2600 + 0.001 )
  7607         return    -1.4005   * (1.0e4);
  7609     else if ( newt <=     0.2610 + 0.001 )
  7611         return    -1.3592   * (1.0e4);
  7613     else if ( newt <=     0.2620 + 0.001 )
  7615         return    -1.3134   * (1.0e4);
  7617     else if ( newt <=     0.2630 + 0.001 )
  7619         return    -1.2641   * (1.0e4);
  7621     else if ( newt <=     0.2640 + 0.001 )
  7623         return    -1.2126   * (1.0e4);
  7625     else if ( newt <=     0.2650 + 0.001 )
  7627         return    -1.1473   * (1.0e4);
  7629     else if ( newt <=     0.2660 + 0.001 )
  7631         return    -1.0855   * (1.0e4);
  7633     else if ( newt <=     0.2670 + 0.001 )
  7635         return    -1.0260   * (1.0e4);
  7637     else if ( newt <=     0.2680 + 0.001 )
  7639         return    -0.9694   * (1.0e4);
  7641     else if ( newt <=     0.2690 + 0.001 )
  7643         return    -0.9197   * (1.0e4);
  7645     else if ( newt <=     0.2700 + 0.001 )
  7647         return    -0.8715   * (1.0e4);
  7649     else if ( newt <=     0.2710 + 0.001 )
  7651         return    -0.8246   * (1.0e4);
  7653     else if ( newt <=     0.2720 + 0.001 )
  7655         return    -0.7783   * (1.0e4);
  7657     else if ( newt <=     0.2730 + 0.001 )
  7659         return    -0.7117   * (1.0e4);
  7661     else if ( newt <=     0.2740 + 0.001 )
  7663         return    -0.6489   * (1.0e4);
  7665     else if ( newt <=     0.2750 + 0.001 )
  7667         return    -0.5866   * (1.0e4);
  7669     else if ( newt <=     0.2760 + 0.001 )
  7671         return    -0.5240   * (1.0e4);
  7673     else if ( newt <=     0.2770 + 0.001 )
  7675         return    -0.4540   * (1.0e4);
  7677     else if ( newt <=     0.2780 + 0.001 )
  7679         return    -0.3857   * (1.0e4);
  7681     else if ( newt <=     0.2790 + 0.001 )
  7683         return    -0.3182   * (1.0e4);
  7685     else if ( newt <=     0.2800 + 0.001 )
  7687         return    -0.2513   * (1.0e4);
  7689     else if ( newt <=     0.2810 + 0.001 )
  7691         return    -0.2075   * (1.0e4);
  7693     else if ( newt <=     0.2820 + 0.001 )
  7695         return    -0.15802  * (1.0e4);
  7697     else if ( newt <=     0.2830 + 0.001 )
  7699         return    -0.10562  * (1.0e4);
  7701     else if ( newt <=     0.2840 + 0.001 )
  7703         return    -0.05011  * (1.0e4);
  7705     else if ( newt <=     0.2850 + 0.001 )
  7707         return     0.00771  * (1.0e4);
  7709     else if ( newt <=     0.2860 + 0.001 )
  7711         return     0.06971  * (1.0e4);
  7713     else if ( newt <=     0.2870 + 0.001 )
  7715         return     0.13571  * (1.0e4);
  7717     else if ( newt <=     0.2880 + 0.001 )
  7719         return     0.20521  * (1.0e4);
  7721     else if ( newt <=     0.2890 + 0.001 )
  7723         return     0.29581  * (1.0e4);
  7725     else if ( newt <=     0.2900 + 0.001 )
  7727         return     0.38341  * (1.0e4);
  7729     else if ( newt <=     0.2910 + 0.001 )
  7731         return     0.46981  * (1.0e4);
  7733     else if ( newt <=     0.2920 + 0.001 )
  7735         return     0.55461  * (1.0e4);
  7737     else if ( newt <=     0.2930 + 0.001 )
  7739         return     0.60371  * (1.0e4);
  7741     else if ( newt <=     0.2940 + 0.001 )
  7743         return     0.66001  * (1.0e4);
  7745     else if ( newt <=     0.2950 + 0.001 )
  7747         return     0.71861  * (1.0e4);
  7749     else if ( newt <=     0.2960 + 0.001 )
  7751         return     0.77981  * (1.0e4);
  7753     else if ( newt <=     0.2970 + 0.001 )
  7755         return     0.83161  * (1.0e4);
  7757     else if ( newt <=     0.2980 + 0.001 )
  7759         return     0.89061  * (1.0e4);
  7761     else if ( newt <=     0.2990 + 0.001 )
  7763         return     0.95491  * (1.0e4);
  7765     else if ( newt <=     0.3000 + 0.001 )
  7767         return     1.02421  * (1.0e4);
  7769     else if ( newt <=     0.3010 + 0.001 )
  7771         return     1.08901  * (1.0e4);
  7773     else if ( newt <=     0.3020 + 0.001 )
  7775         return     1.15901  * (1.0e4);
  7777     else if ( newt <=     0.3030 + 0.001 )
  7779         return     1.23141  * (1.0e4);
  7781     else if ( newt <=     0.3040 + 0.001 )
  7783         return     1.30471  * (1.0e4);
  7785     else if ( newt <=     0.3050 + 0.001 )
  7787         return     1.38741  * (1.0e4);
  7789     else if ( newt <=     0.3060 + 0.001 )
  7791         return     1.46602  * (1.0e4);
  7793     else if ( newt <=     0.3070 + 0.001 )
  7795         return     1.54152  * (1.0e4);
  7797     else if ( newt <=     0.3080 + 0.001 )
  7799         return     1.61382  * (1.0e4);
  7801     else if ( newt <=     0.3090 + 0.001 )
  7803         return     1.68312  * (1.0e4);
  7805     else if ( newt <=     0.3100 + 0.001 )
  7807         return     1.74982  * (1.0e4);
  7809     else if ( newt <=     0.3110 + 0.001 )
  7811         return     1.81452  * (1.0e4);
  7813     else if ( newt <=     0.3120 + 0.001 )
  7815         return     1.8778   * (1.0e4);
  7817     else if ( newt <=     0.3130 + 0.001 )
  7819         return     1.9319   * (1.0e4);
  7821     else if ( newt <=     0.3140 + 0.001 )
  7823         return     1.9880   * (1.0e4);
  7825     else if ( newt <=     0.3150 + 0.001 )
  7827         return     2.0449   * (1.0e4);
  7829     else if ( newt <=     0.3160 + 0.001 )
  7831         return     2.1028   * (1.0e4);
  7833     else if ( newt <=     0.3170 + 0.001 )
  7835         return     2.1629   * (1.0e4);
  7837     else if ( newt <=     0.3180 + 0.001 )
  7839         return     2.2232   * (1.0e4);
  7841     else if ( newt <=     0.3190 + 0.001 )
  7843         return     2.2834   * (1.0e4);
  7845     else if ( newt <=     0.3200 + 0.001 )
  7847         return     2.3434   * (1.0e4);
  7849     else if ( newt <=     0.3210 + 0.001 )
  7851         return     2.4023   * (1.0e4);
  7853     else if ( newt <=     0.3220 + 0.001 )
  7855         return     2.4598   * (1.0e4);
  7857     else if ( newt <=     0.3230 + 0.001 )
  7859         return     2.5156   * (1.0e4);
  7861     else if ( newt <=     0.3240 + 0.001 )
  7863         return     2.5693   * (1.0e4);
  7865     else if ( newt <=     0.3250 + 0.001 )
  7867         return     2.6180   * (1.0e4);
  7869     else if ( newt <=     0.3260 + 0.001 )
  7871         return     2.6651   * (1.0e4);
  7873     else if ( newt <=     0.3270 + 0.001 )
  7875         return     2.7099   * (1.0e4);
  7877     else if ( newt <=     0.3280 + 0.001 )
  7879         return     2.7522   * (1.0e4);
  7881     else if ( newt <=     0.3290 + 0.001 )
  7883         return     2.7945   * (1.0e4);
  7885     else if ( newt <=     0.3300 + 0.001 )
  7887         return     2.8325   * (1.0e4);
  7889     else if ( newt <=     0.3310 + 0.001 )
  7891         return     2.8664   * (1.0e4);
  7893     else if ( newt <=     0.3320 + 0.001 )
  7895         return     2.8961   * (1.0e4);
  7897     else if ( newt <=     0.3330 + 0.001 )
  7899         return     2.9207   * (1.0e4);
  7901     else if ( newt <=     0.3340 + 0.001 )
  7903         return     2.9428   * (1.0e4);
  7905     else if ( newt <=     0.3350 + 0.001 )
  7907         return     2.9629   * (1.0e4);
  7909     else if ( newt <=     0.3360 + 0.001 )
  7911         return     2.9814   * (1.0e4);
  7913     else if ( newt <=     0.3370 + 0.001 )
  7915         return     3.0024   * (1.0e4);
  7917     else if ( newt <=     0.3380 + 0.001 )
  7919         return     3.0201   * (1.0e4);
  7921     else if ( newt <=     0.3390 + 0.001 )
  7923         return     3.0343   * (1.0e4);
  7925     else if ( newt <=     0.3400 + 0.001 )
  7927         return     3.0438   * (1.0e4);
  7929     else if ( newt <=     0.3410 + 0.001 )
  7931         return     3.0450   * (1.0e4);
  7933     else if ( newt <=     0.3420 + 0.001 )
  7935         return     3.0401   * (1.0e4);
  7937     else if ( newt <=     0.3430 + 0.001 )
  7939         return     3.0280   * (1.0e4);
  7941     else if ( newt <=     0.3440 + 0.001 )
  7943         return     3.0092   * (1.0e4);
  7945     else if ( newt <=     0.3450 + 0.001 )
  7947         return     2.9883   * (1.0e4);
  7949     else if ( newt <=     0.3460 + 0.001 )
  7951         return     2.9642   * (1.0e4);
  7953     else if ( newt <=     0.3470 + 0.001 )
  7955         return     2.9395   * (1.0e4);
  7957     else if ( newt <=     0.3480 + 0.001 )
  7959         return     2.9157   * (1.0e4);
  7961     else if ( newt <=     0.3490 + 0.001 )
  7963         return     2.8961   * (1.0e4);
  7965     else if ( newt <=     0.3500 + 0.001 )
  7967         return     2.8770   * (1.0e4);
  7969     else if ( newt <=     0.3510 + 0.001 )
  7971         return     2.8581   * (1.0e4);
  7973     else if ( newt <=     0.3520 + 0.001 )
  7975         return     2.8384   * (1.0e4);
  7977     else if ( newt <=     0.3530 + 0.001 )
  7979         return     2.8114   * (1.0e4);
  7981     else if ( newt <=     0.3540 + 0.001 )
  7983         return     2.7836   * (1.0e4);
  7985     else if ( newt <=     0.3550 + 0.001 )
  7987         return     2.7540   * (1.0e4);
  7989     else if ( newt <=     0.3560 + 0.001 )
  7991         return     2.7225   * (1.0e4);
  7993     else if ( newt <=     0.3570 + 0.001 )
  7995         return     2.6901   * (1.0e4);
  7997     else if ( newt <=     0.3580 + 0.001 )
  7999         return     2.6578   * (1.0e4);
  8001     else if ( newt <=     0.3590 + 0.001 )
  8003         return     2.6268   * (1.0e4);
  8005     else if ( newt <=     0.3600 + 0.001 )
  8007         return     2.5974   * (1.0e4);
  8009     else if ( newt <=     0.3610 + 0.001 )
  8011         return     2.5780   * (1.0e4);
  8013     else if ( newt <=     0.3620 + 0.001 )
  8015         return     2.5576   * (1.0e4);
  8017     else if ( newt <=     0.3630 + 0.001 )
  8019         return     2.5370   * (1.0e4);
  8021     else if ( newt <=     0.3640 + 0.001 )
  8023         return     2.5156   * (1.0e4);
  8025     else if ( newt <=     0.3650 + 0.001 )
  8027         return     2.4947   * (1.0e4);
  8029     else if ( newt <=     0.3660 + 0.001 )
  8031         return     2.4724   * (1.0e4);
  8033     else if ( newt <=     0.3670 + 0.001 )
  8035         return     2.4493   * (1.0e4);
  8037     else if ( newt <=     0.3680 + 0.001 )
  8039         return     2.4258   * (1.0e4);
  8041     else if ( newt <=     0.3690 + 0.001 )
  8043         return     2.3965   * (1.0e4);
  8045     else if ( newt <=     0.3700 + 0.001 )
  8047         return     2.3691   * (1.0e4);
  8049     else if ( newt <=     0.3710 + 0.001 )
  8051         return     2.3429   * (1.0e4);
  8053     else if ( newt <=     0.3720 + 0.001 )
  8055         return     2.3181   * (1.0e4);
  8057     else if ( newt <=     0.3730 + 0.001 )
  8059         return     2.2823   * (1.0e4);
  8061     else if ( newt <=     0.3740 + 0.001 )
  8063         return     2.2514   * (1.0e4);
  8065     else if ( newt <=     0.3750 + 0.001 )
  8067         return     2.2231   * (1.0e4);
  8069     else if ( newt <=     0.3760 + 0.001 )
  8071         return     2.1969   * (1.0e4);
  8073     else if ( newt <=     0.3770 + 0.001 )
  8075         return     2.1722   * (1.0e4);
  8077     else if ( newt <=     0.3780 + 0.001 )
  8079         return     2.1484   * (1.0e4);
  8081     else if ( newt <=     0.3790 + 0.001 )
  8083         return     2.1245   * (1.0e4);
  8085     else if ( newt <=     0.3800 + 0.001 )
  8087         return     2.1001   * (1.0e4);
  8089     else if ( newt <=     0.3810 + 0.001 )
  8091         return     2.0776   * (1.0e4);
  8093     else if ( newt <=     0.3820 + 0.001 )
  8095         return     2.0531   * (1.0e4);
  8097     else if ( newt <=     0.3830 + 0.001 )
  8099         return     2.0271   * (1.0e4);
  8101     else if ( newt <=     0.3840 + 0.001 )
  8103         return     1.9996   * (1.0e4);
  8105     else if ( newt <=     0.3850 + 0.001 )
  8107         return     1.9708   * (1.0e4);
  8109     else if ( newt <=     0.3860 + 0.001 )
  8111         return     1.9408   * (1.0e4);
  8113     else if ( newt <=     0.3870 + 0.001 )
  8115         return     1.9099   * (1.0e4);
  8117     else if ( newt <=     0.3880 + 0.001 )
  8119         return     1.8783   * (1.0e4);
  8121     else if ( newt <=     0.3890 + 0.001 )
  8123         return     1.8422   * (1.0e4);
  8125     else if ( newt <=     0.3900 + 0.001 )
  8127         return     1.8070   * (1.0e4);
  8129     else if ( newt <=     0.3910 + 0.001 )
  8131         return     1.7723   * (1.0e4);
  8133     else if ( newt <=     0.3920 + 0.001 )
  8135         return     1.7380   * (1.0e4);
  8137     else if ( newt <=     0.3930 + 0.001 )
  8139         return     1.7041   * (1.0e4);
  8141     else if ( newt <=     0.3940 + 0.001 )
  8143         return     1.6704   * (1.0e4);
  8145     else if ( newt <=     0.3950 + 0.001 )
  8147         return     1.6368   * (1.0e4);
  8149     else if ( newt <=     0.3960 + 0.001 )
  8151         return     1.6027   * (1.0e4);
  8153     else if ( newt <=     0.3970 + 0.001 )
  8155         return     1.5765   * (1.0e4);
  8157     else if ( newt <=     0.3980 + 0.001 )
  8159         return     1.5471   * (1.0e4);
  8161     else if ( newt <=     0.3990 + 0.001 )
  8163         return     1.5157   * (1.0e4);
  8165     else if ( newt <=     0.4000 + 0.001 )
  8167         return     1.4824   * (1.0e4);
  8169     else if ( newt <=     0.4010 + 0.001 )
  8171         return     1.4560   * (1.0e4);
  8173     else if ( newt <=     0.4020 + 0.001 )
  8175         return     1.4260   * (1.0e4);
  8177     else if ( newt <=     0.4030 + 0.001 )
  8179         return     1.3944   * (1.0e4);
  8181     else if ( newt <=     0.4040 + 0.001 )
  8183         return     1.3617   * (1.0e4);
  8185     else if ( newt <=     0.4050 + 0.001 )
  8187         return     1.3271   * (1.0e4);
  8189     else if ( newt <=     0.4060 + 0.001 )
  8191         return     1.2930   * (1.0e4);
  8193     else if ( newt <=     0.4070 + 0.001 )
  8195         return     1.2593   * (1.0e4);
  8197     else if ( newt <=     0.4080 + 0.001 )
  8199         return     1.2264   * (1.0e4);
  8201     else if ( newt <=     0.4090 + 0.001 )
  8203         return     1.1862   * (1.0e4);
  8205     else if ( newt <=     0.4100 + 0.001 )
  8207         return     1.1493   * (1.0e4);
  8209     else if ( newt <=     0.4110 + 0.001 )
  8211         return     1.1147   * (1.0e4);
  8213     else if ( newt <=     0.4120 + 0.001 )
  8215         return     1.0823   * (1.0e4);
  8217     else if ( newt <=     0.4130 + 0.001 )
  8219         return     1.0478   * (1.0e4);
  8221     else if ( newt <=     0.4140 + 0.001 )
  8223         return     1.0168   * (1.0e4);
  8225     else if ( newt <=     0.4150 + 0.001 )
  8227         return     0.9883   * (1.0e4);
  8229     else if ( newt <=     0.4160 + 0.001 )
  8231         return     0.9621   * (1.0e4);
  8233     else if ( newt <=     0.4170 + 0.001 )
  8235         return     0.9347   * (1.0e4);
  8237     else if ( newt <=     0.4180 + 0.001 )
  8239         return     0.9092   * (1.0e4);
  8241     else if ( newt <=     0.4190 + 0.001 )
  8243         return     0.8848   * (1.0e4);
  8245     else if ( newt <=     0.4200 + 0.001 )
  8247         return     0.8611   * (1.0e4);
  8249     else if ( newt <=     0.4210 + 0.001 )
  8251         return     0.8372   * (1.0e4);
  8253     else if ( newt <=     0.4220 + 0.001 )
  8255         return     0.8156   * (1.0e4);
  8257     else if ( newt <=     0.4230 + 0.001 )
  8259         return     0.7968   * (1.0e4);
  8261     else if ( newt <=     0.4240 + 0.001 )
  8263         return     0.7811   * (1.0e4);
  8265     else if ( newt <=     0.4250 + 0.001 )
  8267         return     0.7648   * (1.0e4);
  8269     else if ( newt <=     0.4260 + 0.001 )
  8271         return     0.7525   * (1.0e4);
  8273     else if ( newt <=     0.4270 + 0.001 )
  8275         return     0.7431   * (1.0e4);
  8277     else if ( newt <=     0.4280 + 0.001 )
  8279         return     0.7359   * (1.0e4);
  8281     else if ( newt <=     0.4290 + 0.001 )
  8283         return     0.7315   * (1.0e4);
  8285     else if ( newt <=     0.4300 + 0.001 )
  8287         return     0.7271   * (1.0e4);
  8289     else if ( newt <=     0.4310 + 0.001 )
  8291         return     0.7222   * (1.0e4);
  8293     else if ( newt <=     0.4320 + 0.001 )
  8295         return     0.7166   * (1.0e4);
  8297     else if ( newt <=     0.4330 + 0.001 )
  8299         return     0.7158   * (1.0e4);
  8301     else if ( newt <=     0.4340 + 0.001 )
  8303         return     0.7130   * (1.0e4);
  8305     else if ( newt <=     0.4350 + 0.001 )
  8307         return     0.7092   * (1.0e4);
  8309     else if ( newt <=     0.4360 + 0.001 )
  8311         return     0.7047   * (1.0e4);
  8313     else if ( newt <=     0.4370 + 0.001 )
  8315         return     0.7121   * (1.0e4);
  8317     else if ( newt <=     0.4380 + 0.001 )
  8319         return     0.7147   * (1.0e4);
  8321     else if ( newt <=     0.4390 + 0.001 )
  8323         return     0.7143   * (1.0e4);
  8325     else if ( newt <=     0.4400 + 0.001 )
  8327         return     0.7113   * (1.0e4);
  8329     else if ( newt <=     0.4410 + 0.001 )
  8331         return     0.7017   * (1.0e4);
  8333     else if ( newt <=     0.4420 + 0.001 )
  8335         return     0.6919   * (1.0e4);
  8337     else if ( newt <=     0.4430 + 0.001 )
  8339         return     0.6821   * (1.0e4);
  8341     else if ( newt <=     0.4440 + 0.001 )
  8343         return     0.6731   * (1.0e4);
  8345     else if ( newt <=     0.4450 + 0.001 )
  8347         return     0.6662   * (1.0e4);
  8349     else if ( newt <=     0.4460 + 0.001 )
  8351         return     0.6611   * (1.0e4);
  8353     else if ( newt <=     0.4470 + 0.001 )
  8355         return     0.6585   * (1.0e4);
  8357     else if ( newt <=     0.4480 + 0.001 )
  8359         return     0.6582   * (1.0e4);
  8361     else if ( newt <=     0.4490 + 0.001 )
  8363         return     0.6559   * (1.0e4);
  8365     else if ( newt <=     0.4500 + 0.001 )
  8367         return     0.6572   * (1.0e4);
  8369     else if ( newt <=     0.4510 + 0.001 )
  8371         return     0.6613   * (1.0e4);
  8373     else if ( newt <=     0.4520 + 0.001 )
  8375         return     0.6679   * (1.0e4);
  8377     else if ( newt <=     0.4530 + 0.001 )
  8379         return     0.6770   * (1.0e4);
  8381     else if ( newt <=     0.4540 + 0.001 )
  8383         return     0.6883   * (1.0e4);
  8385     else if ( newt <=     0.4550 + 0.001 )
  8387         return     0.7018   * (1.0e4);
  8389     else if ( newt <=     0.4560 + 0.001 )
  8391         return     0.7170   * (1.0e4);
  8393     else if ( newt <=     0.4570 + 0.001 )
  8395         return     0.7279   * (1.0e4);
  8397     else if ( newt <=     0.4580 + 0.001 )
  8399         return     0.7396   * (1.0e4);
  8401     else if ( newt <=     0.4590 + 0.001 )
  8403         return     0.7502   * (1.0e4);
  8405     else if ( newt <=     0.4600 + 0.001 )
  8407         return     0.7594   * (1.0e4);
  8409     else if ( newt <=     0.4610 + 0.001 )
  8411         return     0.7671   * (1.0e4);
  8413     else if ( newt <=     0.4620 + 0.001 )
  8415         return     0.7735   * (1.0e4);
  8417     else if ( newt <=     0.4630 + 0.001 )
  8419         return     0.7790   * (1.0e4);
  8421     else if ( newt <=     0.4640 + 0.001 )
  8423         return     0.7842   * (1.0e4);
  8425     else if ( newt <=     0.4650 + 0.001 )
  8427         return     0.7864   * (1.0e4);
  8429     else if ( newt <=     0.4660 + 0.001 )
  8431         return     0.7893   * (1.0e4);
  8433     else if ( newt <=     0.4670 + 0.001 )
  8435         return     0.7925   * (1.0e4);
  8437     else if ( newt <=     0.4680 + 0.001 )
  8439         return     0.7959   * (1.0e4);
  8441     else if ( newt <=     0.4690 + 0.001 )
  8443         return     0.7963   * (1.0e4);
  8445     else if ( newt <=     0.4700 + 0.001 )
  8447         return     0.7973   * (1.0e4);
  8449     else if ( newt <=     0.4710 + 0.001 )
  8451         return     0.7983   * (1.0e4);
  8453     else if ( newt <=     0.4720 + 0.001 )
  8455         return     0.7989   * (1.0e4);
  8457     else if ( newt <=     0.4730 + 0.001 )
  8459         return     0.7977   * (1.0e4);
  8461     else if ( newt <=     0.4740 + 0.001 )
  8463         return     0.7960   * (1.0e4);
  8465     else if ( newt <=     0.4750 + 0.001 )
  8467         return     0.7937   * (1.0e4);
  8469     else if ( newt <=     0.4760 + 0.001 )
  8471         return     0.7905   * (1.0e4);
  8473     else if ( newt <=     0.4770 + 0.001 )
  8475         return     0.7949   * (1.0e4);
  8477     else if ( newt <=     0.4780 + 0.001 )
  8479         return     0.7963   * (1.0e4);
  8481     else if ( newt <=     0.4790 + 0.001 )
  8483         return     0.7959   * (1.0e4);
  8485     else if ( newt <=     0.4800 + 0.001 )
  8487         return     0.7940   * (1.0e4);
  8489     else if ( newt <=     0.4810 + 0.001 )
  8491         return     0.7962   * (1.0e4);
  8493     else if ( newt <=     0.4820 + 0.001 )
  8495         return     0.79562  * (1.0e4);
  8497     else if ( newt <=     0.4830 + 0.001 )
  8499         return     0.79342  * (1.0e4);
  8501     else if ( newt <=     0.4840 + 0.001 )
  8503         return     0.78971  * (1.0e4);
  8505     else if ( newt <=     0.4850 + 0.001 )
  8507         return     0.79061  * (1.0e4);
  8509     else if ( newt <=     0.4860 + 0.001 )
  8511         return     0.78931  * (1.0e4);
  8513     else if ( newt <=     0.4870 + 0.001 )
  8515         return     0.78711  * (1.0e4);
  8517     else if ( newt <=     0.4880 + 0.001 )
  8519         return     0.78441  * (1.0e4);
  8521     else if ( newt <=     0.4890 + 0.001 )
  8523         return     0.77531  * (1.0e4);
  8525     else if ( newt <=     0.4900 + 0.001 )
  8527         return     0.76761  * (1.0e4);
  8529     else if ( newt <=     0.4910 + 0.001 )
  8531         return     0.76051  * (1.0e4);
  8533     else if ( newt <=     0.4920 + 0.001 )
  8535         return     0.75401  * (1.0e4);
  8537     else if ( newt <=     0.4930 + 0.001 )
  8539         return     0.74501  * (1.0e4);
  8541     else if ( newt <=     0.4940 + 0.001 )
  8543         return     0.73741  * (1.0e4);
  8545     else if ( newt <=     0.4950 + 0.001 )
  8547         return     0.73031  * (1.0e4);
  8549     else if ( newt <=     0.4960 + 0.001 )
  8551         return     0.72381  * (1.0e4);
  8553     else if ( newt <=     0.4970 + 0.001 )
  8555         return     0.70921  * (1.0e4);
  8557     else if ( newt <=     0.4980 + 0.001 )
  8559         return     0.69701  * (1.0e4);
  8561     else if ( newt <=     0.4990 + 0.001 )
  8563         return     0.68591  * (1.0e4);
  8565     else if ( newt <=     0.5000 + 0.001 )
  8567         return     0.67541  * (1.0e4);
  8569     else if ( newt <=     0.5010 + 0.001 )
  8571         return     0.67231  * (1.0e4);
  8573     else if ( newt <=     0.5020 + 0.001 )
  8575         return     0.66751  * (1.0e4);
  8577     else if ( newt <=     0.5030 + 0.001 )
  8579         return     0.66181  * (1.0e4);
  8581     else if ( newt <=     0.5040 + 0.001 )
  8583         return     0.65491  * (1.0e4);
  8585     else if ( newt <=     0.5050 + 0.001 )
  8587         return     0.65251  * (1.0e4);
  8589     else if ( newt <=     0.5060 + 0.001 )
  8591         return     0.64732  * (1.0e4);
  8593     else if ( newt <=     0.5070 + 0.001 )
  8595         return     0.64022  * (1.0e4);
  8597     else if ( newt <=     0.5080 + 0.001 )
  8599         return     0.63142  * (1.0e4);
  8601     else if ( newt <=     0.5090 + 0.001 )
  8603         return     0.62132  * (1.0e4);
  8605     else if ( newt <=     0.5100 + 0.001 )
  8607         return     0.61002  * (1.0e4);
  8609     else if ( newt <=     0.5110 + 0.001 )
  8611         return     0.59812  * (1.0e4);
  8613     else if ( newt <=     0.5120 + 0.001 )
  8615         return     0.5858   * (1.0e4);
  8617     else if ( newt <=     0.5130 + 0.001 )
  8619         return     0.5707   * (1.0e4);
  8621     else if ( newt <=     0.5140 + 0.001 )
  8623         return     0.5569   * (1.0e4);
  8625     else if ( newt <=     0.5150 + 0.001 )
  8627         return     0.5439   * (1.0e4);
  8629     else if ( newt <=     0.5160 + 0.001 )
  8631         return     0.5318   * (1.0e4);
  8633     else if ( newt <=     0.5170 + 0.001 )
  8635         return     0.5291   * (1.0e4);
  8637     else if ( newt <=     0.5180 + 0.001 )
  8639         return     0.5248   * (1.0e4);
  8641     else if ( newt <=     0.5190 + 0.001 )
  8643         return     0.5199   * (1.0e4);
  8645     else if ( newt <=     0.5200 + 0.001 )
  8647         return     0.5141   * (1.0e4);
  8649     else if ( newt <=     0.5210 + 0.001 )
  8651         return     0.5099   * (1.0e4);
  8653     else if ( newt <=     0.5220 + 0.001 )
  8655         return     0.5037   * (1.0e4);
  8657     else if ( newt <=     0.5230 + 0.001 )
  8659         return     0.4958   * (1.0e4);
  8661     else if ( newt <=     0.5240 + 0.001 )
  8663         return     0.4864   * (1.0e4);
  8665     else if ( newt <=     0.5250 + 0.001 )
  8667         return     0.4759   * (1.0e4);
  8669     else if ( newt <=     0.5260 + 0.001 )
  8671         return     0.4649   * (1.0e4);
  8673     else if ( newt <=     0.5270 + 0.001 )
  8675         return     0.4539   * (1.0e4);
  8677     else if ( newt <=     0.5280 + 0.001 )
  8679         return     0.4434   * (1.0e4);
  8681     else if ( newt <=     0.5290 + 0.001 )
  8683         return     0.4309   * (1.0e4);
  8685     else if ( newt <=     0.5300 + 0.001 )
  8687         return     0.4204   * (1.0e4);
  8689     else if ( newt <=     0.5310 + 0.001 )
  8691         return     0.4115   * (1.0e4);
  8693     else if ( newt <=     0.5320 + 0.001 )
  8695         return     0.4040   * (1.0e4);
  8697     else if ( newt <=     0.5330 + 0.001 )
  8699         return     0.4048   * (1.0e4);
  8701     else if ( newt <=     0.5340 + 0.001 )
  8703         return     0.4048   * (1.0e4);
  8705     else if ( newt <=     0.5350 + 0.001 )
  8707         return     0.4047   * (1.0e4);
  8709     else if ( newt <=     0.5360 + 0.001 )
  8711         return     0.4043   * (1.0e4);
  8713     else if ( newt <=     0.5370 + 0.001 )
  8715         return     0.3908   * (1.0e4);
  8717     else if ( newt <=     0.5380 + 0.001 )
  8719         return     0.3801   * (1.0e4);
  8721     else if ( newt <=     0.5390 + 0.001 )
  8723         return     0.3705   * (1.0e4);
  8725     else if ( newt <=     0.5400 + 0.001 )
  8727         return     0.3620   * (1.0e4);
  8729     else if ( newt <=     0.5410 + 0.001 )
  8731         return     0.3447   * (1.0e4);
  8733     else if ( newt <=     0.5420 + 0.001 )
  8735         return     0.3312   * (1.0e4);
  8737     else if ( newt <=     0.5430 + 0.001 )
  8739         return     0.3200   * (1.0e4);
  8741     else if ( newt <=     0.5440 + 0.001 )
  8743         return     0.3107   * (1.0e4);
  8745     else if ( newt <=     0.5450 + 0.001 )
  8747         return     0.3246   * (1.0e4);
  8749     else if ( newt <=     0.5460 + 0.001 )
  8751         return     0.3339   * (1.0e4);
  8753     else if ( newt <=     0.5470 + 0.001 )
  8755         return     0.3412   * (1.0e4);
  8757     else if ( newt <=     0.5480 + 0.001 )
  8759         return     0.3461   * (1.0e4);
  8761     else if ( newt <=     0.5490 + 0.001 )
  8763         return     0.3512   * (1.0e4);
  8765     else if ( newt <=     0.5500 + 0.001 )
  8767         return     0.3529   * (1.0e4);
  8769     else if ( newt <=     0.5510 + 0.001 )
  8771         return     0.3520   * (1.0e4);
  8773     else if ( newt <=     0.5520 + 0.001 )
  8775         return     0.3489   * (1.0e4);
  8777     else if ( newt <=     0.5530 + 0.001 )
  8779         return     0.3303   * (1.0e4);
  8781     else if ( newt <=     0.5540 + 0.001 )
  8783         return     0.3151   * (1.0e4);
  8785     else if ( newt <=     0.5550 + 0.001 )
  8787         return     0.3019   * (1.0e4);
  8789     else if ( newt <=     0.5560 + 0.001 )
  8791         return     0.2913   * (1.0e4);
  8793     else if ( newt <=     0.5570 + 0.001 )
  8795         return     0.2778   * (1.0e4);
  8797     else if ( newt <=     0.5580 + 0.001 )
  8799         return     0.2686   * (1.0e4);
  8801     else if ( newt <=     0.5590 + 0.001 )
  8803         return     0.2624   * (1.0e4);
  8805     else if ( newt <=     0.5600 + 0.001 )
  8807         return     0.2583   * (1.0e4);
  8809     else if ( newt <=     0.5610 + 0.001 )
  8811         return     0.2423   * (1.0e4);
  8813     else if ( newt <=     0.5620 + 0.001 )
  8815         return     0.2298   * (1.0e4);
  8817     else if ( newt <=     0.5630 + 0.001 )
  8819         return     0.2179   * (1.0e4);
  8821     else if ( newt <=     0.5640 + 0.001 )
  8823         return     0.2057   * (1.0e4);
  8825     else if ( newt <=     0.5650 + 0.001 )
  8827         return     0.2154   * (1.0e4);
  8829     else if ( newt <=     0.5660 + 0.001 )
  8831         return     0.2183   * (1.0e4);
  8833     else if ( newt <=     0.5670 + 0.001 )
  8835         return     0.2177   * (1.0e4);
  8837     else if ( newt <=     0.5680 + 0.001 )
  8839         return     0.2140   * (1.0e4);
  8841     else if ( newt <=     0.5690 + 0.001 )
  8843         return     0.2090   * (1.0e4);
  8845     else if ( newt <=     0.5700 + 0.001 )
  8847         return     0.2016   * (1.0e4);
  8849     else if ( newt <=     0.5710 + 0.001 )
  8851         return     0.1924   * (1.0e4);
  8853     else if ( newt <=     0.5720 + 0.001 )
  8855         return     0.1819   * (1.0e4);
  8857     else if ( newt <=     0.5730 + 0.001 )
  8859         return     0.1777   * (1.0e4);
  8861     else if ( newt <=     0.5740 + 0.001 )
  8863         return     0.1711   * (1.0e4);
  8865     else if ( newt <=     0.5750 + 0.001 )
  8867         return     0.1635   * (1.0e4);
  8869     else if ( newt <=     0.5760 + 0.001 )
  8871         return     0.1554   * (1.0e4);
  8873     else if ( newt <=     0.5770 + 0.001 )
  8875         return     0.1343   * (1.0e4);
  8877     else if ( newt <=     0.5780 + 0.001 )
  8879         return     0.1166   * (1.0e4);
  8881     else if ( newt <=     0.5790 + 0.001 )
  8883         return     0.1007   * (1.0e4);
  8885     else if ( newt <=     0.5800 + 0.001 )
  8887         return     0.0864   * (1.0e4);
  8889     else if ( newt <=     0.5810 + 0.001 )
  8891         return     0.0822   * (1.0e4);
  8893     else if ( newt <=     0.5820 + 0.001 )
  8895         return     0.0772   * (1.0e4);
  8897     else if ( newt <=     0.5830 + 0.001 )
  8899         return     0.0723   * (1.0e4);
  8901     else if ( newt <=     0.5840 + 0.001 )
  8903         return     0.0672   * (1.0e4);
  8905     else if ( newt <=     0.5850 + 0.001 )
  8907         return     0.0492   * (1.0e4);
  8909     else if ( newt <=     0.5860 + 0.001 )
  8911         return     0.0340   * (1.0e4);
  8913     else if ( newt <=     0.5870 + 0.001 )
  8915         return     0.0197   * (1.0e4);
  8917     else if ( newt <=     0.5880 + 0.001 )
  8919         return     0.0061   * (1.0e4);
  8921     else if ( newt <=     0.5890 + 0.001 )
  8923         return    -0.0124   * (1.0e4);
  8925     else if ( newt <=     0.5900 + 0.001 )
  8927         return    -0.0287   * (1.0e4);
  8929     else if ( newt <=     0.5910 + 0.001 )
  8931         return    -0.0436   * (1.0e4);
  8933     else if ( newt <=     0.5920 + 0.001 )
  8935         return    -0.0571   * (1.0e4);
  8937     else if ( newt <=     0.5930 + 0.001 )
  8939         return    -0.0666   * (1.0e4);
  8941     else if ( newt <=     0.5940 + 0.001 )
  8943         return    -0.0756   * (1.0e4);
  8945     else if ( newt <=     0.5950 + 0.001 )
  8947         return    -0.0842   * (1.0e4);
  8949     else if ( newt <=     0.5960 + 0.001 )
  8951         return    -0.0926   * (1.0e4);
  8953     else if ( newt <=     0.5970 + 0.001 )
  8955         return    -0.0928   * (1.0e4);
  8957     else if ( newt <=     0.5980 + 0.001 )
  8959         return    -0.0958   * (1.0e4);
  8961     else if ( newt <=     0.5990 + 0.001 )
  8963         return    -0.1007   * (1.0e4);
  8965     else if ( newt <=     0.6000 + 0.001 )
  8967         return    -0.1074   * (1.0e4);
  8969     else if ( newt <=     0.6010 + 0.001 )
  8971         return    -0.1172   * (1.0e4);
  8973     else if ( newt <=     0.6020 + 0.001 )
  8975         return    -0.1281   * (1.0e4);
  8977     else if ( newt <=     0.6030 + 0.001 )
  8979         return    -0.1400   * (1.0e4);
  8981     else if ( newt <=     0.6040 + 0.001 )
  8983         return    -0.1526   * (1.0e4);
  8985     else if ( newt <=     0.6050 + 0.001 )
  8987         return    -0.1613   * (1.0e4);
  8989     else if ( newt <=     0.6060 + 0.001 )
  8991         return    -0.1714   * (1.0e4);
  8993     else if ( newt <=     0.6070 + 0.001 )
  8995         return    -0.1822   * (1.0e4);
  8997     else if ( newt <=     0.6080 + 0.001 )
  8999         return    -0.1936   * (1.0e4);
  9001     else if ( newt <=     0.6090 + 0.001 )
  9003         return    -0.2069   * (1.0e4);
  9005     else if ( newt <=     0.6100 + 0.001 )
  9007         return    -0.2203   * (1.0e4);
  9009     else if ( newt <=     0.6110 + 0.001 )
  9011         return    -0.2340   * (1.0e4);
  9013     else if ( newt <=     0.6120 + 0.001 )
  9015         return    -0.2481   * (1.0e4);
  9017     else if ( newt <=     0.6130 + 0.001 )
  9019         return    -0.2638   * (1.0e4);
  9021     else if ( newt <=     0.6140 + 0.001 )
  9023         return    -0.2792   * (1.0e4);
  9025     else if ( newt <=     0.6150 + 0.001 )
  9027         return    -0.2945   * (1.0e4);
  9029     else if ( newt <=     0.6160 + 0.001 )
  9031         return    -0.3096   * (1.0e4);
  9033     else if ( newt <=     0.6170 + 0.001 )
  9035         return    -0.3256   * (1.0e4);
  9037     else if ( newt <=     0.6180 + 0.001 )
  9039         return    -0.3408   * (1.0e4);
  9041     else if ( newt <=     0.6190 + 0.001 )
  9043         return    -0.3552   * (1.0e4);
  9045     else if ( newt <=     0.6200 + 0.001 )
  9047         return    -0.3690   * (1.0e4);
  9049     else if ( newt <=     0.6210 + 0.001 )
  9051         return    -0.3825   * (1.0e4);
  9053     else if ( newt <=     0.6220 + 0.001 )
  9055         return    -0.3961   * (1.0e4);
  9057     else if ( newt <=     0.6230 + 0.001 )
  9059         return    -0.4101   * (1.0e4);
  9061     else if ( newt <=     0.6240 + 0.001 )
  9063         return    -0.4247   * (1.0e4);
  9065     else if ( newt <=     0.6250 + 0.001 )
  9067         return    -0.4414   * (1.0e4);
  9069     else if ( newt <=     0.6260 + 0.001 )
  9071         return    -0.4583   * (1.0e4);
  9073     else if ( newt <=     0.6270 + 0.001 )
  9075         return    -0.4756   * (1.0e4);
  9077     else if ( newt <=     0.6280 + 0.001 )
  9079         return    -0.4932   * (1.0e4);
  9081     else if ( newt <=     0.6290 + 0.001 )
  9083         return    -0.5024   * (1.0e4);
  9085     else if ( newt <=     0.6300 + 0.001 )
  9087         return    -0.5139   * (1.0e4);
  9089     else if ( newt <=     0.6310 + 0.001 )
  9091         return    -0.5264   * (1.0e4);
  9093     else if ( newt <=     0.6320 + 0.001 )
  9095         return    -0.5396   * (1.0e4);
  9097     else if ( newt <=     0.6330 + 0.001 )
  9099         return    -0.5564   * (1.0e4);
  9101     else if ( newt <=     0.6340 + 0.001 )
  9103         return    -0.5729   * (1.0e4);
  9105     else if ( newt <=     0.6350 + 0.001 )
  9107         return    -0.5893   * (1.0e4);
  9109     else if ( newt <=     0.6360 + 0.001 )
  9111         return    -0.6055   * (1.0e4);
  9113     else if ( newt <=     0.6370 + 0.001 )
  9115         return    -0.6200   * (1.0e4);
  9117     else if ( newt <=     0.6380 + 0.001 )
  9119         return    -0.6345   * (1.0e4);
  9121     else if ( newt <=     0.6390 + 0.001 )
  9123         return    -0.6490   * (1.0e4);
  9125     else if ( newt <=     0.6400 + 0.001 )
  9127         return    -0.6633   * (1.0e4);
  9129     else if ( newt <=     0.6410 + 0.001 )
  9131         return    -0.6805   * (1.0e4);
  9133     else if ( newt <=     0.6420 + 0.001 )
  9135         return    -0.6971   * (1.0e4);
  9137     else if ( newt <=     0.6430 + 0.001 )
  9139         return    -0.7137   * (1.0e4);
  9141     else if ( newt <=     0.6440 + 0.001 )
  9143         return    -0.7303   * (1.0e4);
  9145     else if ( newt <=     0.6450 + 0.001 )
  9147         return    -0.7457   * (1.0e4);
  9149     else if ( newt <=     0.6460 + 0.001 )
  9151         return    -0.7617   * (1.0e4);
  9153     else if ( newt <=     0.6470 + 0.001 )
  9155         return    -0.7782   * (1.0e4);
  9157     else if ( newt <=     0.6480 + 0.001 )
  9159         return    -0.7951   * (1.0e4);
  9161     else if ( newt <=     0.6490 + 0.001 )
  9163         return    -0.8137   * (1.0e4);
  9165     else if ( newt <=     0.6500 + 0.001 )
  9167         return    -0.8322   * (1.0e4);
  9169     else if ( newt <=     0.6510 + 0.001 )
  9171         return    -0.8507   * (1.0e4);
  9173     else if ( newt <=     0.6520 + 0.001 )
  9175         return    -0.8692   * (1.0e4);
  9177     else if ( newt <=     0.6530 + 0.001 )
  9179         return    -0.8831   * (1.0e4);
  9181     else if ( newt <=     0.6540 + 0.001 )
  9183         return    -0.8979   * (1.0e4);
  9185     else if ( newt <=     0.6550 + 0.001 )
  9187         return    -0.9127   * (1.0e4);
  9189     else if ( newt <=     0.6560 + 0.001 )
  9191         return    -0.9275   * (1.0e4);
  9193     else if ( newt <=     0.6570 + 0.001 )
  9195         return    -0.9410   * (1.0e4);
  9197     else if ( newt <=     0.6580 + 0.001 )
  9199         return    -0.9551   * (1.0e4);
  9201     else if ( newt <=     0.6590 + 0.001 )
  9203         return    -0.9695   * (1.0e4);
  9205     else if ( newt <=     0.6600 + 0.001 )
  9207         return    -0.9843   * (1.0e4);
  9209     else if ( newt <=     0.6610 + 0.001 )
  9211         return    -0.9978   * (1.0e4);
  9213     else if ( newt <=     0.6620 + 0.001 )
  9215         return    -1.0116   * (1.0e4);
  9217     else if ( newt <=     0.6630 + 0.001 )
  9219         return    -1.0252   * (1.0e4);
  9221     else if ( newt <=     0.6640 + 0.001 )
  9223         return    -1.0384   * (1.0e4);
  9225     else if ( newt <=     0.6650 + 0.001 )
  9227         return    -1.0553   * (1.0e4);
  9229     else if ( newt <=     0.6660 + 0.001 )
  9231         return    -1.0704   * (1.0e4);
  9233     else if ( newt <=     0.6670 + 0.001 )
  9235         return    -1.0846   * (1.0e4);
  9237     else if ( newt <=     0.6680 + 0.001 )
  9239         return    -1.0977   * (1.0e4);
  9241     else if ( newt <=     0.6690 + 0.001 )
  9243         return    -1.1129   * (1.0e4);
  9245     else if ( newt <=     0.6700 + 0.001 )
  9247         return    -1.1267   * (1.0e4);
  9249     else if ( newt <=     0.6710 + 0.001 )
  9251         return    -1.1397   * (1.0e4);
  9253     else if ( newt <=     0.6720 + 0.001 )
  9255         return    -1.1522   * (1.0e4);
  9257     else if ( newt <=     0.6730 + 0.001 )
  9259         return    -1.1684   * (1.0e4);
  9261     else if ( newt <=     0.6740 + 0.001 )
  9263         return    -1.1833   * (1.0e4);
  9265     else if ( newt <=     0.6750 + 0.001 )
  9267         return    -1.1975   * (1.0e4);
  9269     else if ( newt <=     0.6760 + 0.001 )
  9271         return    -1.2113   * (1.0e4);
  9273     else if ( newt <=     0.6770 + 0.001 )
  9275         return    -1.2220   * (1.0e4);
  9277     else if ( newt <=     0.6780 + 0.001 )
  9279         return    -1.2334   * (1.0e4);
  9281     else if ( newt <=     0.6790 + 0.001 )
  9283         return    -1.2450   * (1.0e4);
  9285     else if ( newt <=     0.6800 + 0.001 )
  9287         return    -1.2571   * (1.0e4);
  9289     else if ( newt <=     0.6810 + 0.001 )
  9291         return    -1.2653   * (1.0e4);
  9293     else if ( newt <=     0.6820 + 0.001 )
  9295         return    -1.27522  * (1.0e4);
  9297     else if ( newt <=     0.6830 + 0.001 )
  9299         return    -1.28592  * (1.0e4);
  9301     else if ( newt <=     0.6840 + 0.001 )
  9303         return    -1.29731  * (1.0e4);
  9305     else if ( newt <=     0.6850 + 0.001 )
  9307         return    -1.31221  * (1.0e4);
  9309     else if ( newt <=     0.6860 + 0.001 )
  9311         return    -1.32661  * (1.0e4);
  9313     else if ( newt <=     0.6870 + 0.001 )
  9315         return    -1.34091  * (1.0e4);
  9317     else if ( newt <=     0.6880 + 0.001 )
  9319         return    -1.35481  * (1.0e4);
  9321     else if ( newt <=     0.6890 + 0.001 )
  9323         return    -1.35561  * (1.0e4);
  9325     else if ( newt <=     0.6900 + 0.001 )
  9327         return    -1.35921  * (1.0e4);
  9329     else if ( newt <=     0.6910 + 0.001 )
  9331         return    -1.36381  * (1.0e4);
  9333     else if ( newt <=     0.6920 + 0.001 )
  9335         return    -1.36911  * (1.0e4);
  9337     else if ( newt <=     0.6930 + 0.001 )
  9339         return    -1.37531  * (1.0e4);
  9341     else if ( newt <=     0.6940 + 0.001 )
  9343         return    -1.38211  * (1.0e4);
  9345     else if ( newt <=     0.6950 + 0.001 )
  9347         return    -1.38951  * (1.0e4);
  9349     else if ( newt <=     0.6960 + 0.001 )
  9351         return    -1.39721  * (1.0e4);
  9353     else if ( newt <=     0.6970 + 0.001 )
  9355         return    -1.34821  * (1.0e4);
  9357     else if ( newt <=     0.6980 + 0.001 )
  9359         return    -1.31321  * (1.0e4);
  9361     else if ( newt <=     0.6990 + 0.001 )
  9363         return    -1.28351  * (1.0e4);
  9365     else if ( newt <=     0.7000 + 0.001 )
  9367         return    -1.25771  * (1.0e4);
  9369     else if ( newt <=     0.7010 + 0.001 )
  9371         return    -1.28881  * (1.0e4);
  9373     else if ( newt <=     0.7020 + 0.001 )
  9375         return    -1.30841  * (1.0e4);
  9377     else if ( newt <=     0.7030 + 0.001 )
  9379         return    -1.32321  * (1.0e4);
  9381     else if ( newt <=     0.7040 + 0.001 )
  9383         return    -1.33291  * (1.0e4);
  9385     else if ( newt <=     0.7050 + 0.001 )
  9387         return    -1.37221  * (1.0e4);
  9389     else if ( newt <=     0.7060 + 0.001 )
  9391         return    -1.39612  * (1.0e4);
  9393     else if ( newt <=     0.7070 + 0.001 )
  9395         return    -1.40982  * (1.0e4);
  9397     else if ( newt <=     0.7080 + 0.001 )
  9399         return    -1.41422  * (1.0e4);
  9401     else if ( newt <=     0.7090 + 0.001 )
  9403         return    -1.40682  * (1.0e4);
  9405     else if ( newt <=     0.7100 + 0.001 )
  9407         return    -1.39472  * (1.0e4);
  9409     else if ( newt <=     0.7110 + 0.001 )
  9411         return    -1.37952  * (1.0e4);
  9413     else if ( newt <=     0.7120 + 0.001 )
  9415         return    -1.3640   * (1.0e4);
  9417     else if ( newt <=     0.7130 + 0.001 )
  9419         return    -1.3503   * (1.0e4);
  9421     else if ( newt <=     0.7140 + 0.001 )
  9423         return    -1.3404   * (1.0e4);
  9425     else if ( newt <=     0.7150 + 0.001 )
  9427         return    -1.3350   * (1.0e4);
  9429     else if ( newt <=     0.7160 + 0.001 )
  9431         return    -1.3346   * (1.0e4);
  9433     else if ( newt <=     0.7170 + 0.001 )
  9435         return    -1.3301   * (1.0e4);
  9437     else if ( newt <=     0.7180 + 0.001 )
  9439         return    -1.3313   * (1.0e4);
  9441     else if ( newt <=     0.7190 + 0.001 )
  9443         return    -1.3356   * (1.0e4);
  9445     else if ( newt <=     0.7200 + 0.001 )
  9447         return    -1.3417   * (1.0e4);
  9449     else if ( newt <=     0.7210 + 0.001 )
  9451         return    -1.3385   * (1.0e4);
  9453     else if ( newt <=     0.7220 + 0.001 )
  9455         return    -1.3378   * (1.0e4);
  9457     else if ( newt <=     0.7230 + 0.001 )
  9459         return    -1.3377   * (1.0e4);
  9461     else if ( newt <=     0.7240 + 0.001 )
  9463         return    -1.3380   * (1.0e4);
  9465     else if ( newt <=     0.7250 + 0.001 )
  9467         return    -1.3303   * (1.0e4);
  9469     else if ( newt <=     0.7260 + 0.001 )
  9471         return    -1.3255   * (1.0e4);
  9473     else if ( newt <=     0.7270 + 0.001 )
  9475         return    -1.3226   * (1.0e4);
  9477     else if ( newt <=     0.7280 + 0.001 )
  9479         return    -1.3214   * (1.0e4);
  9481     else if ( newt <=     0.7290 + 0.001 )
  9483         return    -1.2975   * (1.0e4);
  9485     else if ( newt <=     0.7300 + 0.001 )
  9487         return    -1.2807   * (1.0e4);
  9489     else if ( newt <=     0.7310 + 0.001 )
  9491         return    -1.2667   * (1.0e4);
  9493     else if ( newt <=     0.7320 + 0.001 )
  9495         return    -1.2546   * (1.0e4);
  9497     else if ( newt <=     0.7330 + 0.001 )
  9499         return    -1.2608   * (1.0e4);
  9501     else if ( newt <=     0.7340 + 0.001 )
  9503         return    -1.2631   * (1.0e4);
  9505     else if ( newt <=     0.7350 + 0.001 )
  9507         return    -1.2635   * (1.0e4);
  9509     else if ( newt <=     0.7360 + 0.001 )
  9511         return    -1.2615   * (1.0e4);
  9513     else if ( newt <=     0.7370 + 0.001 )
  9515         return    -1.2526   * (1.0e4);
  9517     else if ( newt <=     0.7380 + 0.001 )
  9519         return    -1.2416   * (1.0e4);
  9521     else if ( newt <=     0.7390 + 0.001 )
  9523         return    -1.2278   * (1.0e4);
  9525     else if ( newt <=     0.7400 + 0.001 )
  9527         return    -1.2112   * (1.0e4);
  9529     else if ( newt <=     0.7410 + 0.001 )
  9531         return    -1.1652   * (1.0e4);
  9533     else if ( newt <=     0.7420 + 0.001 )
  9535         return    -1.1240   * (1.0e4);
  9537     else if ( newt <=     0.7430 + 0.001 )
  9539         return    -1.0842   * (1.0e4);
  9541     else if ( newt <=     0.7440 + 0.001 )
  9543         return    -1.0457   * (1.0e4);
  9545     else if ( newt <=     0.7450 + 0.001 )
  9547         return    -1.0325   * (1.0e4);
  9549     else if ( newt <=     0.7460 + 0.001 )
  9551         return    -1.0144   * (1.0e4);
  9553     else if ( newt <=     0.7470 + 0.001 )
  9555         return    -0.9944   * (1.0e4);
  9557     else if ( newt <=     0.7480 + 0.001 )
  9559         return    -0.9724   * (1.0e4);
  9561     else if ( newt <=     0.7490 + 0.001 )
  9563         return    -0.9675   * (1.0e4);
  9565     else if ( newt <=     0.7500 + 0.001 )
  9567         return    -0.9542   * (1.0e4);
  9569     else if ( newt <=     0.7510 + 0.001 )
  9571         return    -0.9354   * (1.0e4);
  9573     else if ( newt <=     0.7520 + 0.001 )
  9575         return    -0.9113   * (1.0e4);
  9577     else if ( newt <=     0.7530 + 0.001 )
  9579         return    -0.8879   * (1.0e4);
  9581     else if ( newt <=     0.7540 + 0.001 )
  9583         return    -0.8588   * (1.0e4);
  9585     else if ( newt <=     0.7550 + 0.001 )
  9587         return    -0.8261   * (1.0e4);
  9589     else if ( newt <=     0.7560 + 0.001 )
  9591         return    -0.7912   * (1.0e4);
  9593     else if ( newt <=     0.7570 + 0.001 )
  9595         return    -0.7625   * (1.0e4);
  9597     else if ( newt <=     0.7580 + 0.001 )
  9599         return    -0.7323   * (1.0e4);
  9601     else if ( newt <=     0.7590 + 0.001 )
  9603         return    -0.7026   * (1.0e4);
  9605     else if ( newt <=     0.7600 + 0.001 )
  9607         return    -0.6739   * (1.0e4);
  9609     else if ( newt <=     0.7610 + 0.001 )
  9611         return    -0.6479   * (1.0e4);
  9613     else if ( newt <=     0.7620 + 0.001 )
  9615         return    -0.6224   * (1.0e4);
  9617     else if ( newt <=     0.7630 + 0.001 )
  9619         return    -0.5974   * (1.0e4);
  9621     else if ( newt <=     0.7640 + 0.001 )
  9623         return    -0.5724   * (1.0e4);
  9625     else if ( newt <=     0.7650 + 0.001 )
  9627         return    -0.5516   * (1.0e4);
  9629     else if ( newt <=     0.7660 + 0.001 )
  9631         return    -0.5295   * (1.0e4);
  9633     else if ( newt <=     0.7670 + 0.001 )
  9635         return    -0.5068   * (1.0e4);
  9637     else if ( newt <=     0.7680 + 0.001 )
  9639         return    -0.4838   * (1.0e4);
  9641     else if ( newt <=     0.7690 + 0.001 )
  9643         return    -0.4535   * (1.0e4);
  9645     else if ( newt <=     0.7700 + 0.001 )
  9647         return    -0.4255   * (1.0e4);
  9649     else if ( newt <=     0.7710 + 0.001 )
  9651         return    -0.3988   * (1.0e4);
  9653     else if ( newt <=     0.7720 + 0.001 )
  9655         return    -0.3736   * (1.0e4);
  9657     else if ( newt <=     0.7730 + 0.001 )
  9659         return    -0.3498   * (1.0e4);
  9661     else if ( newt <=     0.7740 + 0.001 )
  9663         return    -0.3275   * (1.0e4);
  9665     else if ( newt <=     0.7750 + 0.001 )
  9667         return    -0.3065   * (1.0e4);
  9669     else if ( newt <=     0.7760 + 0.001 )
  9671         return    -0.2865   * (1.0e4);
  9673     else if ( newt <=     0.7770 + 0.001 )
  9675         return    -0.2772   * (1.0e4);
  9677     else if ( newt <=     0.7780 + 0.001 )
  9679         return    -0.2659   * (1.0e4);
  9681     else if ( newt <=     0.7790 + 0.001 )
  9683         return    -0.2538   * (1.0e4);
  9685     else if ( newt <=     0.7800 + 0.001 )
  9687         return    -0.2409   * (1.0e4);
  9689     else if ( newt <=     0.7810 + 0.001 )
  9691         return    -0.2271   * (1.0e4);
  9693     else if ( newt <=     0.7820 + 0.001 )
  9695         return    -0.2127   * (1.0e4);
  9697     else if ( newt <=     0.7830 + 0.001 )
  9699         return    -0.1979   * (1.0e4);
  9701     else if ( newt <=     0.7840 + 0.001 )
  9703         return    -0.1829   * (1.0e4);
  9705     else if ( newt <=     0.7850 + 0.001 )
  9707         return    -0.1654   * (1.0e4);
  9709     else if ( newt <=     0.7860 + 0.001 )
  9711         return    -0.1494   * (1.0e4);
  9713     else if ( newt <=     0.7870 + 0.001 )
  9715         return    -0.1347   * (1.0e4);
  9717     else if ( newt <=     0.7880 + 0.001 )
  9719         return    -0.1215   * (1.0e4);
  9721     else if ( newt <=     0.7890 + 0.001 )
  9723         return    -0.1818   * (1.0e4);
  9725     else if ( newt <=     0.7900 + 0.001 )
  9727         return    -0.2258   * (1.0e4);
  9729     else if ( newt <=     0.7910 + 0.001 )
  9731         return    -0.2643   * (1.0e4);
  9733     else if ( newt <=     0.7920 + 0.001 )
  9735         return    -0.2985   * (1.0e4);
  9737     else if ( newt <=     0.7930 + 0.001 )
  9739         return    -0.2442   * (1.0e4);
  9741     else if ( newt <=     0.7940 + 0.001 )
  9743         return    -0.2085   * (1.0e4);
  9745     else if ( newt <=     0.7950 + 0.001 )
  9747         return    -0.1803   * (1.0e4);
  9749     else if ( newt <=     0.7960 + 0.001 )
  9751         return    -0.1595   * (1.0e4);
  9753     else if ( newt <=     0.7970 + 0.001 )
  9755         return    -0.1669   * (1.0e4);
  9757     else if ( newt <=     0.7980 + 0.001 )
  9759         return    -0.1773   * (1.0e4);
  9761     else if ( newt <=     0.7990 + 0.001 )
  9763         return    -0.1934   * (1.0e4);
 Real fZero(const Real &, const Real &, const Real &, const Real &, const ID &)
 
Real d0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real aortaFluxIn(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real aortaFlux6_(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real E(const Real &, const Real &, const Real &, const Real &, const ID &)
 
Real aortaFlux6(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real linearFlux8(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real p0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real linearFlux6_(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real linearPress2(Real t, const Real &, const Real &, const Real &, const ID &i)
 
void updateInverseJacobian(const UInt &iQuadPt)
 
Real aortaFlux3_(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real w0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real linearFlux3_(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real linearFlux7(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real f(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real aortaFlux3(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real linearFlux5(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real aortaPhisPress(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
double Real
Generic real data. 
 
Real u1(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real linearFlux9(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real aortaFlux8(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real aortaFlux4(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real u0(const Real &t, const Real &x, const Real &y, const Real &z, const ID &i)
 
Real linearFlux3(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real aortaFlux9(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real aortaFlux7(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real u2(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real linearFlux6(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real aortaFlux5(const Real &t, const Real &x=0, const Real &y=0, const Real &z=0, const ID &i=0)
 
Real linearFlux4(Real t, const Real &, const Real &, const Real &, const ID &i)
 
Real linearFluxIn(Real t, const Real &, const Real &, const Real &, const ID &i)