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)