gmane.comp.gcc.fortran http://blog.gmane.org/gmane.comp.gcc.fortran hourly 1 1901-01-01T00:00+00:00 Gmane http://gmane.org/img/gmane-25t.png http://gmane.org Re: Guru advice regarding Gfortran numeric libraries like AFNL or IMSL http://permalink.gmane.org/gmane.comp.gcc.fortran/38839 <pre> Here's where FreeBSD really shines with its ports system. You can build any or all of: gcc42 gcc44 gcc45 gcc46 gcc47 gcc48 llvm-gcc4 Or, if you are on i386 or amd64, you can install a pre-build package for all of the above. http://www.freebsd.org/cgi/ports.cgi?query=gcc&amp;stype=name&amp;sektion=lang </pre> Anton Shterenlikht 2012-05-24T08:02:49 Re: Guru advice regarding Gfortran numeric libraries like AFNL or IMSL http://permalink.gmane.org/gmane.comp.gcc.fortran/38838 <pre>Thank you, that is exactly what I wanted to know. I guess my best answer is probably to grab some fedora fortran packages and see if I can rebuild them on RHEL6.2. I suspect that road will lead me to dependency hell but it's worth a shot. (I had already enabled the testing repos to see if they had anything newer.) Failing that I can start again and give her a full Fedora 18 VM I guess. Thank you for such a detailed answer. You've given me the answer I need. 3 options. 1. Downgrade to 4.2 2. Upgrade to 4.6+ 3. Fix the library. I think at my current level of knowledge I should play to my strengths and go with 1 or 2. Thanks again for such detailed answers. I wasn't expecting it and it shows why GNU has gone so far. On 24/05/2012 2:38 PM, Tobias Burnus wrote: </pre> Franki Hauptle 2012-05-24T07:20:48 Re: Guru advice regarding Gfortran numeric libraries like AFNL or IMSL http://permalink.gmane.org/gmane.comp.gcc.fortran/38837 <pre>Hi Frank, Tobias Burnus wrote: Just for completeness: There is, of course, another solution: Changing the internal procedure to a noninternal one. I do not know the source code, but might be simple but it can also be difficult. It is relatively simple if the internal procedure doesn't use any local variable or argument from its host procedure. However, I assume that you don't want to change the library code - especially as you do not have Fortran experience. (It is unfortunate that some codes make heavily use of vendor extensions, which reduces the compatibility. It would really help if programmers used flags like -std=f2003 and multiple compilers; that not make code more portable but can also find true bugs.) Tobias </pre> Tobias Burnus 2012-05-24T06:48:48 Re: [OT] cure for nasty #ifdef's? http://permalink.gmane.org/gmane.comp.gcc.fortran/38836 <pre>Hi, Sorry for an off-topic answer... On May 24, 2012, at 1:09 AM, Tom Roche wrote: I'd say that the lack of readability results in: - bad code auditability, - bad code shareability, - bad code "documentability", three things which are arguably crucial for scientific codes nowadays (see e.g. http://www.nature.com/news/2010/101013/pdf/467775a.pdf) Best, Sylwester -- http://www.igf.fuw.edu.pl/~slayoo/ Please consider your environmental responsibility before printing this e-mail. </pre> Sylwester Arabas 2012-05-24T06:40:56 Re: Guru advice regarding Gfortran numeric libraries like AFNL or IMSL http://permalink.gmane.org/gmane.comp.gcc.fortran/38835 <pre>Hi Frank, Frank Hauptle wrote: Internal procedures are something like "bar" in the following example: subroutine foo (x) ... contains subroutine bar () end subroutine bar end subroutine foo "bar" is an internal procedure, contained in "foo". By host-association (here: host = "foo") it has access to "x", the dummy argument of "foo". If one passes "bar" as argument in a procedure call, it can be invoked in a different context - but the "x" has still to be available in "bar". Therefore, the Fortran standard didn't allow to pass such procedures until recently: Only Fortran 2008 allows it. (Intel's ifort compiler allowed it as vendor extension since years.) GCC internally supports this also since years and the early versions of gfortran didn't check this constraint. (The standard mandates that constraints are checked - at least with some -std= flag.) In gfortran, the constraint check was eventually implemented - and, for some reason, internal procedures were completely rejected instead of only for -std=f95 and -std=f2003. Since GCC/gfortran 4.6, internal procedures are (again) supported and are accepted with -std=f2008 (and the default -std=gnu). Thus, you have two options: downgrading the compiler - seemingly 4.2 works; upgrading to 4.6, 4.7 or the developer version 4.8; or patching your version of the compiler. I think there is no Red-Hat-packaged GCC 4.6/4.7 available for RHEL, yet; I believe there will be at some point a technical preview for a newer GCC which can be installed in parallel to 4.4 - but only RH knows when and in which 6.x version. (I am unfortunately not aware of a RH equivalent to SUSE's build server, which offers newer GCC version.) In principle, installing multiple GCC is not a problem and just a PATH and library path (LD_LIBRARY_PATH) problem. On the supercomputer, on which I have an account, they solved it by offering modules: "module load gcc/4.6" loads the wished-for version. However, I admit that installing multiple versions can confuse users (and admins) and modules are also not that easy to handle. gfortran's libgfortran is downwards compatible since GCC 4.3. Thus, GCC 4.6/4.7's libgfortran can also be used by GCC 4.4. - Thus, if GCC 4.6/4.7's libgfortran is in the LD_LIBRARY_PATH or /usr/local/lib*, /usr/bin/gfortran compiled programs will continue to work (using the new libgfortran library). See also: - http://gcc.gnu.org/wiki/GFortranBinaries , http://gcc.gnu.org/wiki/GFortranDistros - http://gcc.gnu.org/wiki/GFortran#news (release notes per version) Sorry for not being able to provide you with a simple solution. Tobias </pre> Tobias Burnus 2012-05-24T06:38:21 Guru advice regarding Gfortran numeric libraries like AFNL or IMSL http://permalink.gmane.org/gmane.comp.gcc.fortran/38834 <pre>Hi, I've been asked by an Academic to setup a gfortran install on Linux with the AFNL numeric libs installed. (I'm sadly a sysadmin, not a Fortran programmer) But no problem, set her up with a Scientific Linux VM, with spice and GCC/Gfortran 4.4.6 installed. (the default for RHEL/SL 6.2) [root&lt; at &gt;druggat afnl]# gfortran --version GNU Fortran (GCC) 4.4.6 20110731 (Red Hat 4.4.6-3) Copyright (C) 2010 Free Software Foundation, Inc. I swapped the makefile from ifort to gfortran and hit it up.. Got this: [root&lt; at &gt;druggat afnl]# make gfortran -c numtypes.o error.o constants.o ./src/min.f90 -o min.o ./src/min.f90:910.28: CALL LineSrch(X2, Fdir, Tol) 1 Error: Internal procedure 'fdir' is not allowed as an actual argument at (1) ./src/min.f90:829.24: CALL LineSrch(X2,Fdir, Tol) (and lots more of the same) Since the make file seems to suggest that it would compile successfully on gfortran 4.2 I'm a bit lost regarding what to do now. Are their any gfortran numerical libraries comparable to AFNL or IMSL kindest regards Franki IT Services University of Western Australia. </pre> Frank Hauptle 2012-05-24T04:26:34 [OT] cure for nasty #ifdef's? http://permalink.gmane.org/gmane.comp.gcc.fortran/38833 <pre> I'm pretty sure this is off-topic for this list (feel free to point me toward something more generally about fortran), but I figure folks on this list are likely to have answers (esp multi-platform coders): I work on an atmospheric model (now a family of models--more below), the code for which long predates me. (We've not yet gotten it to F90.) Like most Earth-system models it defines a 3-space over the surface of an Earth region, and various variables of interest within each "box" in that space, and then solves differential equations over discrete "timesteps." One area in which we advance is by creating "instrumented models" (IMs) derived from our "base model." E.g., we have a sensitivity model which essentially calculates derivatives of change in one quantity vs change(s) in one or more other quantities, in each box over each timestep. (E.g., one might be interested in the sensitivity of production of ozone to concentrations of other chemical species at various spacetimes.) Currently, this IM is developed as a separate body of code, in a separate repository: its developer takes the current version of the "base model", and adds necessary data structures and calculations. The other instrumented models (I believe we currently have 5) proceed similarly. Note each IM has well-separated concerns: * our IMs are purely additive: they only add code to the base model, never subtract. * our IMs are orthogonal: users can't (yet :-) simultaneously run a combination of 2 or more IMs. Separate development means each IM developer spends most of her/his time "catching up" to changes in the base model. So we'd like to, jointly, * have all the developers work on one set of codes but * (continue to) have users run only the model of interest * pay no runtime performance penalty How best to do this? My group seems to concur that 1 We should make the model-choice &lt; at &gt; compiletime. I.e., the user builds only the model(s) of interest. 2 The only way to do this (of which we're aware, in fortran--more below) in one set of code is to use the preprocessor. E.g., for IMs={foo, bar}, #ifdef foo SUBROUTINE VDIFF ( CGRID, JDATE, JTIME, TSTEP, &amp; SENGRID, NPMAX, IPT, IPARM, ISPLIT, &amp; ILAYER, STDATE, IDATE, ITIME, SEN_PAR ) #elif bar SUBROUTINE VDIFF ( CGRID, JDATE, JTIME, TSTEP, &amp; HGGRID, HGMAX, HGMIN, HGDATE, HGTIME, HGPAR ) #else C the base model SUBROUTINE VDIFF ( CGRID, JDATE, JTIME, TSTEP ) #endif 3 For all but the most trivial blocks (like the above), using the preprocessor would tend to degrade maintainability, e.g., the code would quickly get unreadable. I'd appreciate "second opinions" from outside my group about how to do this. (And feel free to say, "your analysis is flawless--you're totally screwed" :-) Regarding 1: can one make the choice of model/code &lt; at &gt; runtime, without paying a penalty at either buildtime (these models already require nontrivial time and expertise to build) or runtime (e.g., bloating the executable)? Regarding 2: is there an alternative to the preprocessor? E.g.: I spent a long time coding java; in a similar circumstance there, I'd try AspectJ or other vehicle for aspect-oriented programming. I believe there are attempts at AOP for C++, but know of nothing for fortran. Am I missing something? Regarding 3: is the readability concern warranted? or will we "just get over it"? Again, if there's a better vehicle for questions like this, please pass a pointer; otherwise, your assistance is appreciated, Tom Roche &lt;Tom_Roche&lt; at &gt;pobox.com&gt; </pre> Tom Roche 2012-05-23T23:09:09 Re: [Patch, fortran] PR 53456 Improve time resolution on targets without gettimeofday http://permalink.gmane.org/gmane.comp.gcc.fortran/38829 <pre> AFAICS it's not necessary. HAVE_CLOCK_GETTIME and HAVE_CLOCK_GETTIME_LIBRT are mutually exclusive, see the test for the latter in configure.ac. Also, POSIX states that CLOCK_REALTIME is mandatory, and in intrinsics/system_clock.c we assume its presense if clock_gettime is available. Thanks for the review. </pre> Janne Blomqvist 2012-05-23T17:24:42 Re: [Patch, fortran] PR 53456 Improve time resolution on targets without gettimeofday http://permalink.gmane.org/gmane.comp.gcc.fortran/38828 <pre>Hi Janne, On 05/23/2012 04:44 PM, Janne Blomqvist wrote: I concur and wrote essentially the same patch - except that I used a different condition (and forgot to change the comment): +#elif defined(HAVE_CLOCK_GETTIME) &amp;&amp; defined(CLOCK_REALTIME) &amp;&amp; !defined(HAVE_CLOCK_GETTIME_LIBRT) I think in practice it shouldn't matter as librt* mostly affects GLIBC which has gettimeofday. Plus I believe CLOCK_REALTIME should be always available and only CLOCK_MONOTONIC is unimplemented on some systems (e.g. on VxWorks 5.5 but not on 6.x according to the documentation). (* = In intrinsics/system_clock.c, weak refs are used for librt; for gf_gettime that's an overkill and a simple time() as fall back is sufficient.) Okay with considering my preprocessor condition. Tobias PS: In Thunderbird, the attachment is shown as inline. Thus, from my side, it is sufficient to just keep the attachment. </pre> Tobias Burnus 2012-05-23T15:36:34 [Patch, fortran] PR 53456 Improve time resolution on targets without gettimeofday http://permalink.gmane.org/gmane.comp.gcc.fortran/38827 <pre>Hi, some targets such as VXWorks don't provide gettimeofday but do provide clock_gettime. The attached patch allows such targets to get better resolution for the DATE_AND_TIME (up to the 1 millisecond limit of the API) intrinsic than the 1 second resolution provided by the current fallback of using the C standard time(). (SYSTEM_CLOCK already uses clock_gettime if available and thus nanosecond resultion has been available to vxworks also before this patch.) Ok for trunk? (Patch both inline below and as an attachment) 2012-05-23 Janne Blomqvist &lt;jb&lt; at &gt;gcc.gnu.org&gt; PR fortran/53456 * intrinsics/time_1.h (gf_gettime): Fallback for clock_gettime. diff --git a/libgfortran/intrinsics/time_1.h b/libgfortran/intrinsics/time_1.h index aaca56a..ca5b26b 100644 --- a/libgfortran/intrinsics/time_1.h +++ b/libgfortran/intrinsics/time_1.h &lt; at &gt;&lt; at &gt; -181,8 +181,8 &lt; at &gt;&lt; at &gt; gf_cputime (long *user_sec, long *user_usec, long *system_sec, long *system_usec #endif -/* Realtime clock with microsecond resolution, falling back to less - precise functions if the target does not support gettimeofday(). +/* Realtime clock with microsecond resolution, falling back to other + functions if the target does not support gettimeofday(). Arguments: secs - OUTPUT, seconds &lt; at &gt;&lt; at &gt; -204,6 +204,12 &lt; at &gt;&lt; at &gt; gf_gettime (time_t * secs, long * usecs) *secs = tv.tv_sec; *usecs = tv.tv_usec; return err; +#elif defined(HAVE_CLOCK_GETTIME) + struct timespec ts; + int err = clock_gettime (CLOCK_REALTIME, &amp;ts); + *secs = ts.tv_sec; + *usecs = ts.tv_nsec / 1000; + return err; #else time_t t = time (NULL); *secs = t; </pre> Janne Blomqvist 2012-05-23T14:44:54 Re: [Patch, libgfortran] PR53445/53444 - Fix compilation on VxWorks http://permalink.gmane.org/gmane.comp.gcc.fortran/38826 <pre>On Tue, May 22, 2012 at 10:33 PM, Janne Blomqvist &lt;blomqvist.janne&lt; at &gt;gmail.com&gt; wrote: Googling some more, it seems the vxworks STATUS is just a typedef for int, so I guess the original patch works. Also the error checking is not useful here, so we could do just #elif defined(HAVE_STRERROR_R_2ARGS) strerror_r (errnum, buf); return buf; </pre> Janne Blomqvist 2012-05-22T19:42:36 Re: [Patch, libgfortran] PR53445/53444 - Fix compilation on VxWorks http://permalink.gmane.org/gmane.comp.gcc.fortran/38825 <pre> For the a) patch (strerror_r): The configure.ac diff occurs twice in the patch, and the patch file has DOS line endings. Also, based on some googling the vxworks 2-arg strerror_r returns OK or ERROR (an enum, I presume). So the trick with builtin_choose_expr is both wrong and unnecessary. Thus I'd replace +#elif defined(HAVE_STRERROR_R_2ARGS) + return + __builtin_choose_expr (__builtin_classify_type (strerror_r (0, buf)) + == 5, + /* char*-returning strerror_r() */ + strerror_r (errnum, buf), + /* int-returning strerror_r () */ + (strerror_r (errnum, buf), buf)); with #elif defined(HAVE_STRERROR_R_2ARGS) if (strerror_r (errnum, buf) == OK) return buf; return NULL; Otherwise Ok. For the b) patch (sticky bit): Also DOS line endings. If the patch is verbatim by the OP it would perhaps be nice to include his name in the ChangeLog as well (it's small enough that AFAICS it goes below the trivial limit, thus no requirement for copyright assignment etc.). Otherwise Ok. Thanks for the patches! </pre> Janne Blomqvist 2012-05-22T19:33:27 Re: [Patch, libgfortran] PR53445/53444 - Fix compilation on VxWorks http://permalink.gmane.org/gmane.comp.gcc.fortran/38824 <pre> I forgot to attach that patch. Now with patch and automake 1.11.1 for the generated files. Tobias </pre> Tobias Burnus 2012-05-22T14:57:20 [Patch, Fortran] PR53389 realloc-on-assignment: Memory leak and wrong double evaluation (4.6-4.8 regression) http://permalink.gmane.org/gmane.comp.gcc.fortran/38821 <pre>First: I have another rather simple patch, which still needs to be reviewed: http://gcc.gnu.org/ml/fortran/2012-05/msg00086.html * * * When a realloc on assignment is done, gfortran first obtains the descriptor by calling realloc_lhs_loop_for_fcn_call and then does the actual function call using gfc_conv_function_expr. However, realloc_lhs_loop_for_fcn_call already causes that the RHS is evaluated (via gfc_add_loop_ss_code). The double evaluation of an inner function call, causes memory leaks - and double evaluation (correctness/performance issue). The following code solves the issue. I am not 100% sure that it correctly does so, but it solved the problem and the testsuite passed. (It might obsolete the check in "se-&gt;direct_byref" of gfc_conv_procedure_call.) The test case is a bit pointless as it will only fail if one uses valgrind on it; thus, one could consider committing the patch even without. However, it might not harm to have it in the test suite. Build and regtested on x86-64-linux. OK for the trunk and 4.6/4.7? Tobias </pre> Tobias Burnus 2012-05-21T17:51:23 Re: [PATCH] Fix bessel_7.f90 failures on sparc-*-linux-gnu http://permalink.gmane.org/gmane.comp.gcc.fortran/38818 <pre> OK. Ciao! Steven </pre> Steven Bosscher 2012-05-17T18:33:10 [PATCH] Fix bessel_7.f90 failures on sparc-*-linux-gnu http://permalink.gmane.org/gmane.comp.gcc.fortran/38817 <pre> I've noticed this failure for some time but never got around to inspecting things. The X=34.53 case seem to need a precision allowance bump. Ok for mainline and 4.7 branch? BTW, for the curious, I've attached a log file of this testcase run with all the print statements uncommented. * gfortran.dg/bessel_7.f90: Bump allowed precision to avoid failure on sparc-*-linux-gnu. diff --git a/gcc/testsuite/gfortran.dg/bessel_7.f90 b/gcc/testsuite/gfortran.dg/bessel_7.f90 index 78f1ff2..7e63ed1 100644 --- a/gcc/testsuite/gfortran.dg/bessel_7.f90 +++ b/gcc/testsuite/gfortran.dg/bessel_7.f90 &lt; at &gt;&lt; at &gt; -16,7 +16,7 &lt; at &gt;&lt; at &gt; implicit none real,parameter :: values(*) = [0.0, 0.5, 1.0, 0.9, 1.8,2.0,3.0,4.0,4.25,8.0,34.53, 475.78] real,parameter :: myeps(size(values)) = epsilon(0.0) &amp; - * [2, 3, 4, 5, 8, 2, 12, 6, 7, 6, 31, 168 ] + * [2, 3, 4, 5, 8, 2, 12, 6, 7, 6, 36, 168 ] ! The following is sufficient for me - the values above are a bit ! more tolerant ! * [0, 0, 0, 3, 3, 0, 9, 0, 2, 1, 22, 130 ] YN for X = 0.00000000 -- Epsilon = 1.19209290E-07 0 -Infinity -Infinity NaN NaN T T 1 -Infinity -Infinity NaN NaN T T 2 -Infinity -Infinity NaN NaN T T 3 -Infinity -Infinity NaN NaN T T 4 -Infinity -Infinity NaN NaN T T 5 -Infinity -Infinity NaN NaN T T 6 -Infinity -Infinity NaN NaN T T 7 -Infinity -Infinity NaN NaN T T 8 -Infinity -Infinity NaN NaN T T 9 -Infinity -Infinity NaN NaN T T 10 -Infinity -Infinity NaN NaN T T 11 -Infinity -Infinity NaN NaN T T 12 -Infinity -Infinity NaN NaN T T 13 -Infinity -Infinity NaN NaN T T 14 -Infinity -Infinity NaN NaN T T 15 -Infinity -Infinity NaN NaN T T 16 -Infinity -Infinity NaN NaN T T 17 -Infinity -Infinity NaN NaN T T 18 -Infinity -Infinity NaN NaN T T 19 -Infinity -Infinity NaN NaN T T 20 -Infinity -Infinity NaN NaN T T 21 -Infinity -Infinity NaN NaN T T 22 -Infinity -Infinity NaN NaN T T 23 -Infinity -Infinity NaN NaN T T 24 -Infinity -Infinity NaN NaN T T 25 -Infinity -Infinity NaN NaN T T 26 -Infinity -Infinity NaN NaN T T 27 -Infinity -Infinity NaN NaN T T 28 -Infinity -Infinity NaN NaN T T 29 -Infinity -Infinity NaN NaN T T 30 -Infinity -Infinity NaN NaN T T 31 -Infinity -Infinity NaN NaN T T 32 -Infinity -Infinity NaN NaN T T 33 -Infinity -Infinity NaN NaN T T 34 -Infinity -Infinity NaN NaN T T 35 -Infinity -Infinity NaN NaN T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 0.500000000 -- Epsilon = 1.19209290E-07 0 -0.444518715E+00 -0.444518715E+00 0.00E+00 -0.0000000000 T T 1 -0.147147238E+01 -0.147147238E+01 0.00E+00 -0.0000000000 T T 2 -0.544137096E+01 -0.544137096E+01 0.00E+00 -0.0000000000 T T 3 -0.420594940E+02 -0.420594940E+02 0.00E+00 -0.0000000000 T T 4 -0.499272552E+03 -0.499272552E+03 0.00E+00 -0.0000000000 T T 5 -0.794630127E+04 -0.794630127E+04 0.00E+00 -0.0000000000 T T 6 -0.158426750E+06 -0.158426766E+06 0.16E-01 -0.8273350596 T T 7 -0.379429575E+07 -0.379429625E+07 0.50E+00 -1.1054235697 T T 8 -0.106081856E+09 -0.106081872E+09 0.16E+02 -1.2652279139 T T 9 -0.339082522E+10 -0.339082573E+10 0.51E+03 -1.2666436434 T T 10 -0.121963626E+12 -0.121963651E+12 0.25E+05 -1.6903271675 T T 11 -0.487515436E+13 -0.487515541E+13 0.10E+07 -1.8042696714 T T 12 -0.214384836E+15 -0.214384870E+15 0.34E+08 -1.3129426241 T T 13 -0.102855973E+17 -0.102855984E+17 0.11E+10 -0.8757098913 T T 14 -0.534636670E+18 -0.534636739E+18 0.69E+11 -1.0782289505 T T 15 -0.299293684E+20 -0.299293728E+20 0.44E+13 -1.2326852083 T T 16 -0.179522742E+22 -0.179522770E+22 0.28E+15 -1.3152557611 T T 17 -0.114864624E+24 -0.114864642E+24 0.18E+17 -1.3155984879 T T 18 -0.780899906E+25 -0.780900079E+25 0.17E+19 -1.8577425480 T T 19 -0.562133056E+27 -0.562133203E+27 0.15E+21 -2.2022190094 T T 20 -0.427143054E+29 -0.427143149E+29 0.94E+22 -1.8548390865 T T 21 -0.341658243E+31 -0.341658304E+31 0.60E+24 -1.4841153622 T T 22 -0.286950207E+33 -0.286950265E+33 0.58E+26 -1.6963843107 T T 23 -0.252482006E+35 -0.252482055E+35 0.50E+28 -1.6452014446 T T 24 -0.232254744E+37 -0.232254791E+37 0.48E+30 -1.7169440985 T T 25 -0.222939297E+39 -0.222939358E+39 0.61E+32 -2.2895181179 T T 26 -Infinity -Infinity NaN NaN T T 27 -Infinity -Infinity NaN NaN T T 28 -Infinity -Infinity NaN NaN T T 29 -Infinity -Infinity NaN NaN T T 30 -Infinity -Infinity NaN NaN T T 31 -Infinity -Infinity NaN NaN T T 32 -Infinity -Infinity NaN NaN T T 33 -Infinity -Infinity NaN NaN T T 34 -Infinity -Infinity NaN NaN T T 35 -Infinity -Infinity NaN NaN T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 1.00000000 -- Epsilon = 1.19209290E-07 0 0.882569551E-01 0.882569551E-01 0.00E+00 0.0000000000 T T 1 -0.781212807E+00 -0.781212807E+00 0.00E+00 -0.0000000000 T T 2 -0.165068257E+01 -0.165068257E+01 0.00E+00 -0.0000000000 T T 3 -0.582151747E+01 -0.582151747E+01 0.00E+00 -0.0000000000 T T 4 -0.332784233E+02 -0.332784195E+02 -0.38E-05 0.9615840316 T T 5 -0.260405884E+03 -0.260405853E+03 -0.31E-04 0.9830806851 T T 6 -0.257078052E+04 -0.257078027E+04 -0.24E-03 0.7966452241 T T 7 -0.305889609E+05 -0.305889570E+05 -0.39E-02 1.0712361336 T T 8 -0.425674688E+06 -0.425674625E+06 -0.63E-01 1.2316635847 T T 9 -0.678020600E+07 -0.678020500E+07 -0.10E+01 1.2372201681 T T 10 -0.121618032E+09 -0.121618016E+09 -0.16E+02 1.1036005020 T T 11 -0.242558054E+10 -0.242558029E+10 -0.26E+03 0.8853483200 T T 12 -0.532411556E+11 -0.532411474E+11 -0.82E+04 1.2907209396 T T 13 -0.127536215E+13 -0.127536189E+13 -0.26E+06 1.7242343426 T T 14 -0.331061741E+14 -0.331061679E+14 -0.63E+07 1.5941605568 T T 15 -0.925697523E+15 -0.925697388E+15 -0.13E+09 1.2162718773 T T 16 -0.277378201E+17 -0.277378158E+17 -0.43E+10 1.2989051342 T T 17 -0.886684522E+18 -0.886684385E+18 -0.14E+12 1.3002612591 T T 18 -0.301195355E+20 -0.301195311E+20 -0.44E+13 1.2249022722 T T 19 -0.108341661E+22 -0.108341640E+22 -0.21E+15 1.6345397234 T T 20 -0.411397116E+23 -0.411397026E+23 -0.90E+16 1.8366162777 T T 21 -0.164450504E+25 -0.164450461E+25 -0.43E+18 2.2053914070 T T 22 -0.690280714E+26 -0.690280530E+26 -0.18E+20 2.2417330742 T T 23 -0.303559060E+28 -0.303558971E+28 -0.89E+21 2.4468517303 T T 24 -0.139568144E+30 -0.139568096E+30 -0.47E+23 2.8383326530 T T 25 -0.669623528E+31 -0.669623286E+31 -0.24E+25 3.0289273262 T T 26 -0.334672177E+33 -0.334672060E+33 -0.12E+27 2.9089829922 T T 27 -0.173962564E+35 -0.173962514E+35 -0.50E+28 2.3877766132 T T 28 -0.939063205E+36 -0.939062888E+36 -0.32E+30 2.8309659958 T T 29 -0.525701448E+38 -0.525701245E+38 -0.20E+32 3.2364602089 T T 30 -Infinity -Infinity NaN NaN T T 31 -Infinity -Infinity NaN NaN T T 32 -Infinity -Infinity NaN NaN T T 33 -Infinity -Infinity NaN NaN T T 34 -Infinity -Infinity NaN NaN T T 35 -Infinity -Infinity NaN NaN T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 0.899999976 -- Epsilon = 1.19209290E-07 0 0.562828034E-02 0.562828034E-02 0.00E+00 0.0000000000 T T 1 -0.873126566E+00 -0.873126566E+00 0.00E+00 -0.0000000000 T T 2 -0.194590962E+01 -0.194590962E+01 0.00E+00 -0.0000000000 T T 3 -0.777536106E+01 -0.777536154E+01 0.48E-06 -0.5144455433 T T 4 -0.498898315E+02 -0.498898354E+02 0.38E-05 -0.6414132714 T T 5 -0.435689819E+03 -0.435689850E+03 0.31E-04 -0.5875740051 T T 6 -0.479110840E+04 -0.479110889E+04 0.49E-03 -0.8549169898 T T 7 -0.634457578E+05 -0.634457617E+05 0.39E-02 -0.5164726973 T T 8 -0.982142938E+06 -0.982143000E+06 0.63E-01 -0.5338204503 T T 9 -0.173968740E+08 -0.173968740E+08 0.00E+00 -0.0000000000 T T 10 -0.346955328E+09 -0.346955328E+09 0.00E+00 -0.0000000000 T T 11 -0.769272166E+10 -0.769272218E+10 0.51E+03 -0.5583156943 T T 12 -0.187697365E+12 -0.187697381E+12 0.16E+05 -0.7322369814 T T 13 -0.499757036E+13 -0.499757089E+13 0.52E+06 -0.8800369501 T T 14 -0.144186565E+15 -0.144186565E+15 0.00E+00 -0.0000000000 T T 15 -0.448080677E+16 -0.448080704E+16 0.27E+09 -0.5025433898 T T 16 -0.149216046E+18 -0.149216063E+18 0.17E+11 -0.9658156037 T T 17 -0.530097865E+19 -0.530097920E+19 0.55E+12 -0.8699687719 T T 18 -0.200109990E+21 -0.200110008E+21 0.18E+14 -0.7374641895 T T 19 -0.799909888E+22 -0.799909944E+22 0.56E+15 -0.5903623104 T T 20 -0.337539640E+24 -0.337539676E+24 0.36E+17 -0.8953954577 T T 21 -0.149937638E+26 -0.149937649E+26 0.12E+19 -0.6450285912 T T 22 -0.699371482E+27 -0.699371482E+27 0.00E+00 -0.0000000000 T T 23 -0.341765029E+29 -0.341765006E+29 -0.24E+22 0.5795513988 T T 24 -0.174609977E+31 -0.174609946E+31 -0.30E+24 1.4519795179 T T 25 -0.930911470E+32 -0.930911277E+32 -0.19E+26 1.7430151701 T T 26 -0.516998454E+34 -0.516998330E+34 -0.12E+28 2.0086314678 T T 27 -0.298617124E+36 -0.298617065E+36 -0.59E+29 1.6692295074 T T 28 -0.179118586E+38 -0.179118548E+38 -0.38E+31 1.7810251713 T T 29 -Infinity -Infinity NaN NaN T T 30 -Infinity -Infinity NaN NaN T T 31 -Infinity -Infinity NaN NaN T T 32 -Infinity -Infinity NaN NaN T T 33 -Infinity -Infinity NaN NaN T T 34 -Infinity -Infinity NaN NaN T T 35 -Infinity -Infinity NaN NaN T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 1.79999995 -- Epsilon = 1.19209290E-07 0 0.477431715E+00 0.477431715E+00 0.00E+00 0.0000000000 T T 1 -0.223664910E+00 -0.223664910E+00 0.00E+00 -0.0000000000 T T 2 -0.725948274E+00 -0.725948274E+00 0.00E+00 -0.0000000000 T T 3 -0.138955355E+01 -0.138955355E+01 0.00E+00 -0.0000000000 T T 4 -0.390589714E+01 -0.390589714E+01 0.00E+00 -0.0000000000 T T 5 -0.159699898E+02 -0.159699898E+02 0.00E+00 -0.0000000000 T T 6 -0.848162689E+02 -0.848162689E+02 0.00E+00 -0.0000000000 T T 7 -0.549471802E+03 -0.549471863E+03 0.61E-04 -0.9318039417 T T 8 -0.418885352E+04 -0.418885400E+04 0.49E-03 -0.9778332114 T T 9 -0.366847852E+05 -0.366847852E+05 0.00E+00 -0.0000000000 T T 10 -0.362659000E+06 -0.362659000E+06 0.00E+00 -0.0000000000 T T 11 -0.399285975E+07 -0.399286000E+07 0.25E+00 -0.5252255797 T T 12 -0.484389640E+08 -0.484389640E+08 0.00E+00 -0.0000000000 T T 13 -0.641860032E+09 -0.641860032E+09 0.00E+00 -0.0000000000 T T 14 -0.922287309E+10 -0.922287206E+10 -0.10E+04 0.9313729405 T T 15 -0.142825062E+12 -0.142825046E+12 -0.16E+05 0.9622887969 T T 16 -0.237119498E+13 -0.237119498E+13 0.00E+00 -0.0000000000 T T 17 -0.420117553E+14 -0.420117553E+14 0.00E+00 -0.0000000000 T T 18 -0.791184180E+15 -0.791184180E+15 0.00E+00 -0.0000000000 T T 19 -0.157816716E+17 -0.157816726E+17 0.11E+10 -0.5707379580 T T 20 -0.332377452E+18 -0.332377486E+18 0.34E+11 -0.8671779037 T T 21 -0.737038433E+19 -0.737038488E+19 0.55E+12 -0.6257049441 T T 22 -0.171643265E+21 -0.171643282E+21 0.18E+14 -0.8597713113 T T 23 -0.418835413E+22 -0.418835441E+22 0.28E+15 -0.5637496710 T T 24 -0.106864078E+24 -0.106864078E+24 0.00E+00 -0.0000000000 T T 25 -0.284552067E+25 -0.284552067E+25 0.00E+00 -0.0000000000 T T 26 -0.789353796E+26 -0.789353796E+26 0.00E+00 -0.0000000000 T T 27 -0.227750999E+28 -0.227750999E+28 0.00E+00 -0.0000000000 T T 28 -0.682463701E+29 -0.682463654E+29 -0.47E+22 0.5804569721 T T 29 -0.212094293E+31 -0.212094293E+31 0.00E+00 -0.0000000000 T T 30 -0.682732539E+32 -0.682732443E+32 -0.97E+25 1.1883077621 T T 31 -0.227365437E+34 -0.227365406E+34 -0.31E+27 1.1418395042 T T 32 -0.782464909E+35 -0.782464810E+35 -0.99E+28 1.0617313385 T T 33 -0.277982403E+37 -0.277982371E+37 -0.32E+30 0.9563396573 T T 34 -0.101848637E+39 -0.101848627E+39 -0.10E+32 0.8352649212 T T 35 -Infinity -Infinity NaN NaN T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 2.00000000 -- Epsilon = 1.19209290E-07 0 0.510375679E+00 0.510375679E+00 0.00E+00 0.0000000000 T T 1 -0.107032441E+00 -0.107032441E+00 0.00E+00 -0.0000000000 T T 2 -0.617408097E+00 -0.617408097E+00 0.00E+00 -0.0000000000 T T 3 -0.112778378E+01 -0.112778378E+01 0.00E+00 -0.0000000000 T T 4 -0.276594329E+01 -0.276594329E+01 0.00E+00 -0.0000000000 T T 5 -0.993598938E+01 -0.993598938E+01 0.00E+00 -0.0000000000 T T 6 -0.469140053E+02 -0.469140053E+02 0.00E+00 -0.0000000000 T T 7 -0.271548035E+03 -0.271548035E+03 0.00E+00 -0.0000000000 T T 8 -0.185392224E+04 -0.185392212E+04 -0.12E-03 0.5523424745 T T 9 -0.145598301E+05 -0.145598291E+05 -0.98E-03 0.5626439452 T T 10 -0.129184547E+06 -0.129184539E+06 -0.78E-02 0.5073052645 T T 11 -0.127728563E+07 -0.127728550E+07 -0.13E+00 0.8209409118 T T 12 -0.139209570E+08 -0.139209550E+08 -0.20E+01 1.2051769495 T T 13 -0.165774192E+09 -0.165774176E+09 -0.16E+02 0.8096418977 T T 14 -0.214114355E+10 -0.214114342E+10 -0.13E+03 0.5014805198 T T 15 -0.298102354E+11 -0.298102354E+11 0.00E+00 -0.0000000000 T T 16 -0.445012378E+12 -0.445012378E+12 0.00E+00 -0.0000000000 T T 17 -0.709038755E+13 -0.709038755E+13 0.00E+00 -0.0000000000 T T 18 -0.120091572E+15 -0.120091572E+15 0.00E+00 -0.0000000000 T T 19 -0.215455786E+16 -0.215455786E+16 0.00E+00 -0.0000000000 T T 20 -0.408165087E+17 -0.408165087E+17 0.00E+00 -0.0000000000 T T 21 -0.814175647E+18 -0.814175647E+18 0.00E+00 -0.0000000000 T T 22 -0.170568723E+20 -0.170568723E+20 0.00E+00 -0.0000000000 T T 23 -0.374437014E+21 -0.374437014E+21 0.00E+00 -0.0000000000 T T 24 -0.859499435E+22 -0.859499435E+22 0.00E+00 -0.0000000000 T T 25 -0.205905422E+24 -0.205905440E+24 0.18E+17 -0.7339084148 T T 26 -0.513904039E+25 -0.513904096E+25 0.58E+18 -0.9409739971 T T 27 -0.133409141E+27 -0.133409159E+27 0.18E+20 -1.1599093676 T T 28 -0.359690790E+28 -0.359690820E+28 0.30E+21 -0.6883357167 T T 29 -0.100580012E+30 -0.100580021E+30 0.94E+22 -0.7877128124 T T 30 -0.291322350E+31 -0.291322380E+31 0.30E+24 -0.8702734709 T T 31 -0.872961272E+32 -0.872961272E+32 0.00E+00 -0.0000000000 T T 32 -0.270326676E+34 -0.270326676E+34 0.00E+00 -0.0000000000 T T 33 -0.864172368E+35 -0.864172368E+35 0.00E+00 -0.0000000000 T T 34 -0.284906564E+37 -0.284906564E+37 0.00E+00 -0.0000000000 T T 35 -0.967818173E+38 -0.967818173E+38 0.00E+00 -0.0000000000 T T 36 -Infinity -Infinity NaN NaN T T 37 -Infinity -Infinity NaN NaN T T 38 -Infinity -Infinity NaN NaN T T 39 -Infinity -Infinity NaN NaN T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 3.00000000 -- Epsilon = 1.19209290E-07 0 0.376850069E+00 0.376850069E+00 0.00E+00 0.0000000000 T T 1 0.324674398E+00 0.324674398E+00 0.00E+00 0.0000000000 T T 2 -0.160400465E+00 -0.160400465E+00 0.00E+00 -0.0000000000 T T 3 -0.538541675E+00 -0.538541675E+00 0.00E+00 -0.0000000000 T T 4 -0.916682899E+00 -0.916682899E+00 0.00E+00 -0.0000000000 T T 5 -0.190594614E+01 -0.190594614E+01 0.00E+00 -0.0000000000 T T 6 -0.543647099E+01 -0.543647099E+01 0.00E+00 -0.0000000000 T T 7 -0.198399372E+02 -0.198399372E+02 0.00E+00 -0.0000000000 T T 8 -0.871499100E+02 -0.871499023E+02 -0.76E-05 0.7343667746 T T 9 -0.444959595E+03 -0.444959564E+03 -0.31E-04 0.5753331184 T T 10 -0.258260767E+04 -0.258260742E+04 -0.24E-03 0.7929969430 T T 11 -0.167724258E+05 -0.167724238E+05 -0.20E-02 0.9768413901 T T 12 -0.120415188E+06 -0.120415172E+06 -0.16E-01 1.0885006189 T T 13 -0.946549063E+06 -0.946548938E+06 -0.13E+00 1.1077883244 T T 14 -0.808301050E+07 -0.808300950E+07 -0.10E+01 1.0378073454 T T 15 -0.744948880E+08 -0.744948640E+08 -0.24E+02 2.7025558949 T T 16 -0.736865856E+09 -0.736865600E+09 -0.26E+03 2.9143481255 T T 17 -0.778540800E+10 -0.778540493E+10 -0.31E+04 3.3100132942 T T 18 -0.874977608E+11 -0.874977280E+11 -0.33E+05 3.1415421963 T T 19 -0.104218775E+13 -0.104218729E+13 -0.46E+06 3.6925120354 T T 20 -0.131135470E+14 -0.131135418E+14 -0.52E+07 3.3538191319 T T 21 -0.173805113E+15 -0.173805029E+15 -0.84E+08 4.0487155914 T T 22 -0.242015797E+16 -0.242015689E+16 -0.11E+10 3.7217402458 T T 23 -0.353218454E+17 -0.353218304E+17 -0.15E+11 3.5700509548 T T 24 -0.539181501E+18 -0.539181226E+18 -0.27E+12 4.2765617371 T T 25 -0.859158221E+19 -0.859157781E+19 -0.44E+13 4.2941436768 T T 26 -0.142653867E+21 -0.142653779E+21 -0.88E+14 5.1724481583 T T 27 -0.246407557E+22 -0.246407417E+22 -0.14E+16 4.7912149429 T T 28 -0.442107071E+23 -0.442106801E+23 -0.27E+17 5.1271200180 T T 29 -0.822802536E+24 -0.822801887E+24 -0.65E+18 6.6117520332 T T 30 -0.158633052E+26 -0.158632926E+26 -0.13E+20 6.7063875198 T T 31 -0.316443306E+27 -0.316443048E+27 -0.26E+21 6.8460764885 T T 32 -0.652396523E+28 -0.652395992E+28 -0.53E+22 6.8310971260 T T 33 -0.138861489E+30 -0.138861366E+30 -0.12E+24 7.4172191620 T T 34 -0.304842886E+31 -0.304842614E+31 -0.27E+25 7.4850726128 T T 35 -0.689588648E+32 -0.689587971E+32 -0.68E+26 8.2354526520 T T 36 -0.160599184E+34 -0.160599029E+34 -0.15E+28 8.0826950073 T T 37 -0.384748447E+35 -0.384748100E+35 -0.35E+29 7.5573697090 T T 38 -0.947440236E+36 -0.947439365E+36 -0.87E+30 7.7163219452 T T 39 -0.239633450E+38 -0.239633221E+38 -0.23E+32 7.9875674248 T T 40 -Infinity -Infinity NaN NaN T T 41 -Infinity -Infinity NaN NaN T T 42 -Infinity -Infinity NaN NaN T T 43 -Infinity -Infinity NaN NaN T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 4.00000000 -- Epsilon = 1.19209290E-07 0 -0.169407371E-01 -0.169407371E-01 0.00E+00 -0.0000000000 T T 1 0.397925735E+00 0.397925735E+00 0.00E+00 0.0000000000 T T 2 0.215903610E+00 0.215903610E+00 0.00E+00 0.0000000000 T T 3 -0.182022125E+00 -0.182022125E+00 0.00E+00 -0.0000000000 T T 4 -0.488936782E+00 -0.488936812E+00 0.30E-07 -0.5113135576 T T 5 -0.795851469E+00 -0.795851469E+00 0.00E+00 -0.0000000000 T T 6 -0.150069189E+01 -0.150069189E+01 0.00E+00 -0.0000000000 T T 7 -0.370622420E+01 -0.370622420E+01 0.00E+00 -0.0000000000 T T 8 -0.114710932E+02 -0.114710922E+02 -0.95E-06 0.6974052191 T T 9 -0.421781502E+02 -0.421781464E+02 -0.38E-05 0.7586866617 T T 10 -0.178330582E+03 -0.178330566E+03 -0.15E-04 0.7177680731 T T 11 -0.849474731E+03 -0.849474670E+03 -0.61E-04 0.6027253866 T T 12 -0.449378027E+04 -0.449378027E+04 0.00E+00 -0.0000000000 T T 13 -0.261132070E+05 -0.261132070E+05 0.00E+00 -0.0000000000 T T 14 -0.165242063E+06 -0.165242063E+06 0.00E+00 -0.0000000000 T T 15 -0.113058125E+07 -0.113058125E+07 0.00E+00 -0.0000000000 T T 16 -0.831411750E+07 -0.831411700E+07 -0.50E+00 0.5044797659 T T 17 -0.653823600E+08 -0.653823560E+08 -0.40E+01 0.5132031441 T T 18 -0.547435968E+09 -0.547435904E+09 -0.64E+02 0.9807008505 T T 19 -0.486154138E+10 -0.486154086E+10 -0.51E+03 0.8834579587 T T 20 -0.456372060E+11 -0.456372019E+11 -0.41E+04 0.7528887391 T T 21 -0.451510534E+12 -0.451510501E+12 -0.33E+05 0.6087962389 T T 22 -0.469522344E+13 -0.469522291E+13 -0.52E+06 0.9367065430 T T 23 -0.511959473E+14 -0.511959431E+14 -0.42E+07 0.6872491837 T T 24 -0.584058174E+15 -0.584058107E+15 -0.67E+08 0.9638593793 T T 25 -0.695750235E+16 -0.695750127E+16 -0.11E+10 1.2946023941 T T 26 -0.863847202E+17 -0.863847116E+17 -0.86E+10 0.8341474533 T T 27 -0.111604388E+19 -0.111604375E+19 -0.14E+12 1.0330431461 T T 28 -0.149802082E+20 -0.149802060E+20 -0.22E+13 1.2314077616 T T 29 -0.208606875E+21 -0.208606840E+21 -0.35E+14 1.4148522615 T T 30 -0.300981953E+22 -0.300981896E+22 -0.56E+15 1.5689865351 T T 31 -0.449386870E+23 -0.449386780E+23 -0.90E+16 1.6813545227 T T 32 -0.693539859E+24 -0.693539643E+24 -0.22E+18 2.6146857738 T T 33 -0.110516992E+26 -0.110516957E+26 -0.35E+19 2.6253175735 T T 34 -0.181659496E+27 -0.181659440E+27 -0.55E+20 2.5554816723 T T 35 -0.307715982E+28 -0.307715864E+28 -0.12E+22 3.2183964252 T T 36 -0.536686373E+29 -0.536686184E+29 -0.19E+23 2.9524939060 T T 37 -0.962958301E+30 -0.962957923E+30 -0.38E+24 3.2910318375 T T 38 -0.177610594E+32 -0.177610522E+32 -0.73E+25 3.4258782864 T T 39 -0.336497171E+33 -0.336497016E+33 -0.15E+27 3.8576080799 T T 40 -0.654393364E+34 -0.654393055E+34 -0.31E+28 3.9672596455 T T 41 -0.130542174E+36 -0.130542104E+36 -0.69E+29 4.4547843933 T T 42 -0.266957075E+37 -0.266956917E+37 -0.16E+31 4.9791827202 T T 43 -0.559304431E+38 -0.559304076E+38 -0.35E+32 5.3235244751 T T 44 -Infinity -Infinity NaN NaN T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 4.25000000 -- Epsilon = 1.19209290E-07 0 -0.111918867E+00 -0.111918867E+00 0.00E+00 -0.0000000000 T T 1 0.358568966E+00 0.358568966E+00 0.00E+00 0.0000000000 T T 2 0.280657202E+00 0.280657202E+00 0.00E+00 0.0000000000 T T 3 -0.944210142E-01 -0.944209993E-01 -0.15E-07 1.3238577843 T T 4 -0.413957447E+00 -0.413957447E+00 0.00E+00 -0.0000000000 T T 5 -0.684792995E+00 -0.684792995E+00 0.00E+00 -0.0000000000 T T 6 -0.119732022E+01 -0.119732022E+01 0.00E+00 -0.0000000000 T T 7 -0.269587588E+01 -0.269587588E+01 0.00E+00 -0.0000000000 T T 8 -0.768321228E+01 -0.768321228E+01 0.00E+00 -0.0000000000 T T 9 -0.262291584E+02 -0.262291603E+02 0.19E-05 -0.6100081205 T T 10 -0.103404991E+03 -0.103404999E+03 0.76E-05 -0.6189256310 T T 11 -0.460382599E+03 -0.460382629E+03 0.31E-04 -0.5560592413 T T 12 -0.227975195E+04 -0.227975220E+04 0.24E-03 -0.8983433247 T T 13 -0.124135107E+05 -0.124135127E+05 0.20E-02 -1.3198522329 T T 14 -0.736617266E+05 -0.736617422E+05 0.16E-01 -1.7793772221 T T 15 -0.472887281E+06 -0.472887406E+06 0.13E+00 -2.2173910141 T T 16 -0.326436625E+07 -0.326436700E+07 0.75E+00 -1.9273132086 T T 17 -0.241058700E+08 -0.241058760E+08 0.60E+01 -2.0879416466 T T 18 -0.189582592E+09 -0.189582640E+09 0.48E+02 -2.1238932610 T T 19 -0.158177024E+10 -0.158177075E+10 0.51E+03 -2.7152915001 T T 20 -0.139533046E+11 -0.139533087E+11 0.41E+04 -2.4624803066 T T 21 -0.129743454E+12 -0.129743495E+12 0.41E+05 -2.6482830048 T T 22 -0.126821728E+13 -0.126821768E+13 0.39E+06 -2.6009225845 T T 23 -0.130000355E+14 -0.130000397E+14 0.42E+07 -2.7064826488 T T 24 -0.139438051E+15 -0.139438093E+15 0.42E+08 -2.5232977867 T T 25 -0.156182969E+16 -0.156183023E+16 0.54E+09 -2.8835408688 T T 26 -0.182350288E+17 -0.182350353E+17 0.64E+10 -2.9637022018 T T 27 -0.221549119E+18 -0.221549205E+18 0.86E+11 -3.2524433136 T T 28 -0.279674197E+19 -0.279674307E+19 0.11E+13 -3.2978987694 T T 29 -0.366296405E+20 -0.366296537E+20 0.13E+14 -3.0216093063 T T 30 -0.497090122E+21 -0.497090298E+21 0.18E+15 -2.9687564373 T T 31 -0.698111310E+22 -0.698111591E+22 0.28E+16 -3.3822445869 T T 32 -0.101345025E+24 -0.101345070E+24 0.45E+17 -3.7277538776 T T 33 -0.151915567E+25 -0.151915639E+25 0.72E+18 -3.9789400101 T T 34 -0.234902499E+26 -0.234902591E+26 0.92E+19 -3.2937602997 T T 35 -0.374324836E+27 -0.374324983E+27 0.15E+21 -3.3071277142 T T 36 -0.614186026E+28 -0.614186262E+28 0.24E+22 -3.2249252796 T T 37 -0.103676024E+30 -0.103676062E+30 0.38E+23 -3.0567593575 T T 38 -0.179904059E+31 -0.179904135E+31 0.76E+24 -3.5231294632 T T 39 -0.320674029E+32 -0.320674174E+32 0.15E+26 -3.7949583530 T T 40 -0.586732127E+33 -0.586732398E+33 0.27E+27 -3.8716642857 T T 41 -0.110123024E+35 -0.110123073E+35 0.50E+28 -3.7719972134 T T 42 -0.211885916E+36 -0.211886015E+36 0.99E+29 -3.9208245277 T T 43 -0.417685042E+37 -0.417685232E+37 0.19E+31 -3.8188428879 T T 44 -0.843079123E+38 -0.843079528E+38 0.41E+32 -4.0361852646 T T 45 -Infinity -Infinity NaN NaN T T 46 -Infinity -Infinity NaN NaN T T 47 -Infinity -Infinity NaN NaN T T 48 -Infinity -Infinity NaN NaN T T 49 -Infinity -Infinity NaN NaN T T 50 -Infinity -Infinity NaN NaN T T 51 -Infinity -Infinity NaN NaN T T 52 -Infinity -Infinity NaN NaN T T 53 -Infinity -Infinity NaN NaN T T 54 -Infinity -Infinity NaN NaN T T 55 -Infinity -Infinity NaN NaN T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 8.00000000 -- Epsilon = 1.19209290E-07 0 0.223521501E+00 0.223521501E+00 0.00E+00 0.0000000000 T T 1 -0.158060491E+00 -0.158060491E+00 0.00E+00 -0.0000000000 T T 2 -0.263036609E+00 -0.263036609E+00 0.00E+00 -0.0000000000 T T 3 0.265421867E-01 0.265421867E-01 0.00E+00 0.0000000000 T T 4 0.282943249E+00 0.282943249E+00 0.00E+00 0.0000000000 T T 5 0.256401062E+00 0.256401062E+00 0.00E+00 0.0000000000 T T 6 0.375580788E-01 0.375580788E-01 0.00E+00 0.0000000000 T T 7 -0.200063944E+00 -0.200063944E+00 0.00E+00 -0.0000000000 T T 8 -0.387669981E+00 -0.387669981E+00 0.00E+00 -0.0000000000 T T 9 -0.575276017E+00 -0.575276017E+00 0.00E+00 -0.0000000000 T T 10 -0.906701088E+00 -0.906701028E+00 -0.60E-07 0.5514496565 T T 11 -0.169147670E+01 -0.169147646E+01 -0.24E-06 1.1823987961 T T 12 -0.374485993E+01 -0.374485922E+01 -0.72E-06 1.6021960974 T T 13 -0.954310322E+01 -0.954310036E+01 -0.29E-05 2.5149052143 T T 14 -0.272702255E+02 -0.272702160E+02 -0.95E-05 2.9336023331 T T 15 -0.859026871E+02 -0.859026566E+02 -0.31E-04 2.9801163673 T T 16 -0.294864838E+03 -0.294864746E+03 -0.92E-04 2.6045832634 T T 17 -0.109355664E+04 -0.109355627E+04 -0.37E-03 2.8091824055 T T 18 -0.435275098E+04 -0.435274951E+04 -0.15E-02 2.8230421543 T T 19 -0.184938223E+05 -0.184938164E+05 -0.59E-02 2.6577522755 T T 20 -0.834929063E+05 -0.834928750E+05 -0.31E-01 3.1397159100 T T 21 -0.398970719E+06 -0.398970563E+06 -0.16E+00 3.2852535248 T T 22 -0.201110338E+07 -0.201110263E+07 -0.75E+00 3.1283602715 T T 23 -0.106620980E+08 -0.106620930E+08 -0.50E+01 3.9338448048 T T 24 -0.592959600E+08 -0.592959320E+08 -0.28E+02 3.9611639977 T T 25 -0.345113664E+09 -0.345113504E+09 -0.16E+03 3.8890876770 T T 26 -0.209766438E+10 -0.209766362E+10 -0.77E+03 3.0712497234 T T 27 -0.132897044E+11 -0.132897004E+11 -0.41E+04 2.5854403973 T T 28 -0.876078367E+11 -0.876078203E+11 -0.16E+05 1.5687974691 T T 29 -0.599965172E+12 -0.599965041E+12 -0.13E+06 1.8326257467 T T 30 -0.426213953E+13 -0.426213874E+13 -0.79E+06 1.5478305817 T T 31 -0.313660812E+14 -0.313660749E+14 -0.63E+07 1.6825996637 T T 32 -0.238824995E+15 -0.238824945E+15 -0.50E+08 1.7678738832 T T 33 -0.187923382E+16 -0.187923355E+16 -0.27E+09 1.1982542276 T T 34 -0.152648539E+17 -0.152648517E+17 -0.21E+10 1.1801226139 T T 35 -0.127872025E+18 -0.127872000E+18 -0.26E+11 1.6905400753 T T 36 -0.110361535E+19 -0.110361514E+19 -0.21E+12 1.5670154095 T T 37 -0.980466624E+19 -0.980466404E+19 -0.22E+13 1.8814250231 T T 38 -0.895895444E+20 -0.895895268E+20 -0.18E+14 1.6472229958 T T 39 -0.841295984E+21 -0.841295843E+21 -0.14E+15 1.4033011198 T T 40 -0.811304614E+22 -0.811304501E+22 -0.11E+16 1.1641414165 T T 41 -0.802891654E+23 -0.802891473E+23 -0.18E+17 1.8821434975 T T 42 -0.814850908E+24 -0.814850692E+24 -0.22E+18 2.2254240513 T T 43 -0.847564537E+25 -0.847564364E+25 -0.17E+19 1.7116230726 T T 44 -0.902983341E+26 -0.902983249E+26 -0.92E+19 0.8568403125 T T 45 -0.984805996E+27 -0.984805996E+27 0.00E+00 -0.0000000000 T T 46 -0.109887685E+29 -0.109887685E+29 0.00E+00 -0.0000000000 T T 47 -0.125386036E+30 -0.125386036E+30 0.00E+00 -0.0000000000 T T 48 -0.146229718E+31 -0.146229703E+31 -0.15E+24 0.8668898344 T T 49 -0.174221806E+32 -0.174221782E+32 -0.24E+25 1.1641716957 T T 50 -0.211959417E+33 -0.211959378E+33 -0.39E+26 1.5310410261 T T 51 -0.263207066E+34 -0.263207004E+34 -0.62E+27 1.9727041721 T T 52 -0.333469405E+35 -0.333469331E+35 -0.74E+28 1.8684641123 T T 53 -0.430878163E+36 -0.430878044E+36 -0.12E+30 2.3136959076 T T 54 -0.567578894E+37 -0.567578704E+37 -0.19E+31 2.8103117943 T T 55 -0.761922712E+38 -0.761922458E+38 -0.25E+32 2.7913129330 T T 56 -Infinity -Infinity NaN NaN T T 57 -Infinity -Infinity NaN NaN T T 58 -Infinity -Infinity NaN NaN T T 59 -Infinity -Infinity NaN NaN T T 60 -Infinity -Infinity NaN NaN T T 61 -Infinity -Infinity NaN NaN T T 62 -Infinity -Infinity NaN NaN T T 63 -Infinity -Infinity NaN NaN T T 64 -Infinity -Infinity NaN NaN T T 65 -Infinity -Infinity NaN NaN T T 66 -Infinity -Infinity NaN NaN T T 67 -Infinity -Infinity NaN NaN T T 68 -Infinity -Infinity NaN NaN T T 69 -Infinity -Infinity NaN NaN T T 70 -Infinity -Infinity NaN NaN T T 71 -Infinity -Infinity NaN NaN T T 72 -Infinity -Infinity NaN NaN T T 73 -Infinity -Infinity NaN NaN T T 74 -Infinity -Infinity NaN NaN T T 75 -Infinity -Infinity NaN NaN T T 76 -Infinity -Infinity NaN NaN T T 77 -Infinity -Infinity NaN NaN T T 78 -Infinity -Infinity NaN NaN T T 79 -Infinity -Infinity NaN NaN T T 80 -Infinity -Infinity NaN NaN T T 81 -Infinity -Infinity NaN NaN T T 82 -Infinity -Infinity NaN NaN T T 83 -Infinity -Infinity NaN NaN T T 84 -Infinity -Infinity NaN NaN T T 85 -Infinity -Infinity NaN NaN T T 86 -Infinity -Infinity NaN NaN T T 87 -Infinity -Infinity NaN NaN T T 88 -Infinity -Infinity NaN NaN T T 89 -Infinity -Infinity NaN NaN T T 90 -Infinity -Infinity NaN NaN T T 91 -Infinity -Infinity NaN NaN T T 92 -Infinity -Infinity NaN NaN T T 93 -Infinity -Infinity NaN NaN T T 94 -Infinity -Infinity NaN NaN T T 95 -Infinity -Infinity NaN NaN T T 96 -Infinity -Infinity NaN NaN T T 97 -Infinity -Infinity NaN NaN T T 98 -Infinity -Infinity NaN NaN T T 99 -Infinity -Infinity NaN NaN T T ** -Infinity -Infinity NaN NaN T T YN for X = 34.5299988 -- Epsilon = 1.19209290E-07 0 0.989498049E-01 0.989498049E-01 0.00E+00 0.0000000000 T T 1 0.944138989E-01 0.944138989E-01 0.00E+00 0.0000000000 T T 2 -0.934812874E-01 -0.934812874E-01 0.00E+00 -0.0000000000 T T 3 -0.105242893E+00 -0.105242893E+00 0.00E+00 -0.0000000000 T T 4 0.751940757E-01 0.751940757E-01 0.00E+00 0.0000000000 T T 5 0.122664049E+00 0.122664049E+00 0.00E+00 0.0000000000 T T 6 -0.396701694E-01 -0.396701656E-01 -0.37E-08 0.7877455950 T T 7 -0.136450380E+00 -0.136450380E+00 0.00E+00 -0.0000000000 T T 8 -0.156528950E-01 -0.156528987E-01 0.37E-08 -1.9964357615 T T 9 0.129197374E+00 0.129197374E+00 0.00E+00 0.0000000000 T T 10 0.830016583E-01 0.830016583E-01 0.00E+00 0.0000000000 T T 11 -0.811222792E-01 -0.811222717E-01 -0.75E-08 0.7704418898 T T 12 -0.134686872E+00 -0.134686857E+00 -0.15E-07 0.9280785918 T T 13 -0.124915345E-01 -0.124915317E-01 -0.28E-08 1.8762706518 T T 14 0.125281140E+00 0.125281125E+00 0.15E-07 0.9977559447 T T 15 0.114080645E+00 0.114080630E+00 0.15E-07 1.0957161188 T T 16 -0.261667632E-01 -0.261667594E-01 -0.37E-08 1.1942631006 T T 17 -0.138330176E+00 -0.138330162E+00 -0.15E-07 0.9036350846 T T 18 -0.110040188E+00 -0.110040180E+00 -0.75E-08 0.5679743290 T T 19 0.236053951E-01 0.236053839E-01 0.11E-07 3.9715497494 T T 20 0.136017740E+00 0.136017725E+00 0.15E-07 0.9189977646 T T 21 0.133959308E+00 0.133959323E+00 -0.15E-07 -0.9331191778 T T 22 0.269214716E-01 0.269215107E-01 -0.39E-07-12.1882266998 T T 23 -0.996545032E-01 -0.996544734E-01 -0.30E-07 2.5086674690 T T 24 -0.159678698E+00 -0.159678712E+00 0.15E-07 -0.7828220129 T T 25 -0.122314148E+00 -0.122314207E+00 0.60E-07 -4.0878343582 T T 26 -0.174341798E-01 -0.174342543E-01 0.75E-07-35.8491210938 T T 27 0.960593745E-01 0.960593224E-01 0.52E-07 4.5544748306 T T 28 0.167657360E+00 0.167657346E+00 0.15E-07 0.7455682158 T T 29 0.175843686E+00 0.175843716E+00 -0.30E-07 -1.4217172861 T T 30 0.127707079E+00 0.127707124E+00 -0.45E-07 -2.9364073277 T T 31 0.460626259E-01 0.460626781E-01 -0.52E-07 -9.4979391098 T T 32 -0.449997857E-01 -0.449997336E-01 -0.52E-07 9.7222681046 T T 33 -0.129467964E+00 -0.129467919E+00 -0.45E-07 2.8964693546 T T 34 -0.202462882E+00 -0.202462852E+00 -0.30E-07 1.2347942591 T T 35 -0.269242615E+00 -0.269242644E+00 0.30E-07 -0.9285305738 T T 36 -0.343351841E+00 -0.343352020E+00 0.18E-06 -4.3686966896 T T 37 -0.446695179E+00 -0.446695507E+00 0.33E-06 -6.1563234329 T T 38 -0.613944530E+00 -0.613945067E+00 0.54E-06 -7.3296523094 T T 39 -0.904587328E+00 -0.904588163E+00 0.83E-06 -7.7383351326 T T 40 -0.142943263E+01 -0.142943406E+01 0.14E-05 -8.3949394226 T T 41 -0.240715909E+01 -0.240716195E+01 0.29E-05 -9.9702587128 T T 42 -0.428696012E+01 -0.428696632E+01 0.62E-05-12.1298074722 T T 43 -0.802158833E+01 -0.802160072E+01 0.12E-04-12.9650135040 T T 44 -0.156915112E+02 -0.156915359E+02 0.25E-04-13.2555751801 T T 45 -0.319683628E+02 -0.319684124E+02 0.50E-04-13.0128650665 T T 46 -0.676317673E+02 -0.676318741E+02 0.11E-03-13.2482118607 T T 47 -0.148226318E+03 -0.148226578E+03 0.26E-03-14.6802539825 T T 48 -0.335880371E+03 -0.335880981E+03 0.61E-03-15.2435226440 T T 49 -0.785585327E+03 -0.785586792E+03 0.15E-02-15.6418399811 T T 50 -0.189369873E+04 -0.189370227E+04 0.35E-02-15.6814804077 T T 51 -0.469862793E+04 -0.469863672E+04 0.88E-02-15.6913890839 T T 52 -0.119858271E+05 -0.119858516E+05 0.24E-01-17.0868473053 T T 53 -0.314011699E+05 -0.314012344E+05 0.64E-01-17.2182121277 T T 54 -0.844093125E+05 -0.844094844E+05 0.17E+00-17.0809593201 T T 55 -0.232607109E+06 -0.232607578E+06 0.47E+00-16.9047279358 T T 56 -0.656592250E+06 -0.656593500E+06 0.13E+01-15.9699726105 T T 57 -0.189708700E+07 -0.189709038E+07 0.34E+01-14.9236974716 T T 58 -0.560659700E+07 -0.560660700E+07 0.10E+02-14.9620313644 T T 59 -0.169377020E+08 -0.169377320E+08 0.30E+02-14.8578739166 T T 60 -0.522749240E+08 -0.522750160E+08 0.92E+02-14.7633295059 T T 61 -0.164730160E+09 -0.164730448E+09 0.29E+03-14.6659183502 T T 62 -0.529743040E+09 -0.529744000E+09 0.96E+03-15.2018299103 T T 63 -0.173761958E+10 -0.173762278E+10 0.32E+04-15.4484596252 T T 64 -0.581083238E+10 -0.581084314E+10 0.11E+05-15.5217542648 T T 65 -0.198026813E+11 -0.198027162E+11 0.35E+05-14.7483959198 T T 66 -0.687431434E+11 -0.687432663E+11 0.12E+06-14.9948358536 T T 67 -0.242986058E+12 -0.242986500E+12 0.44E+06-15.2718706131 T T 68 -0.874208821E+12 -0.874210460E+12 0.16E+07-15.7215242386 T T 69 -0.320017636E+13 -0.320018265E+13 0.63E+07-16.4917659760 T T 70 -0.119153749E+14 -0.119153990E+14 0.24E+08-16.9789142609 T T 71 -0.451100583E+14 -0.451101548E+14 0.96E+08-17.9392490387 T T 72 -0.173593703E+15 -0.173594072E+15 0.37E+09-17.8360424042 T T 73 -0.678825420E+15 -0.678826896E+15 0.15E+10-18.2446002960 T T 74 -0.269662099E+16 -0.269662716E+16 0.62E+10-19.2060337067 T T 75 -0.108792370E+17 -0.108792617E+17 0.25E+11-19.0422897339 T T 76 -0.445633136E+17 -0.445634210E+17 0.11E+12-20.2121391296 T T 77 -0.185287105E+18 -0.185287569E+18 0.46E+12-21.0004367828 T T 78 -0.781796748E+18 -0.781798741E+18 0.20E+13-21.3832588196 T T 79 -0.334672263E+19 -0.334673143E+19 0.88E+13-22.0475330353 T T 80 -0.145319065E+20 -0.145319450E+20 0.38E+14-22.2144298553 T T 81 -0.639890699E+20 -0.639892458E+20 0.18E+15-23.0623683929 T T 82 -0.285677468E+21 -0.285678260E+21 0.79E+15-23.2458934784 T T 83 -0.129283413E+22 -0.129283779E+22 0.37E+16-23.7427062988 T T 84 -0.592951133E+22 -0.592952822E+22 0.17E+17-23.8925247192 T T 85 -0.275562223E+23 -0.275563011E+23 0.79E+17-23.9920883179 T T 86 -0.129736807E+24 -0.129737176E+24 0.37E+18-23.8781299591 T T 87 -0.618685638E+24 -0.618687367E+24 0.17E+19-23.4482727051 T T 88 -0.298787967E+25 -0.298788831E+25 0.86E+19-24.2765960693 T T 89 -0.146105892E+26 -0.146106307E+26 0.42E+20-23.8300209045 T T 90 -0.723287981E+26 -0.723290010E+26 0.20E+21-23.5337467194 T T 91 -0.362429231E+27 -0.362430227E+27 0.10E+22-23.0557975769 T T 92 -0.183795566E+28 -0.183796068E+28 0.50E+22-22.9004230499 T T 93 -0.943148437E+28 -0.943151093E+28 0.27E+23-23.6261024475 T T 94 -0.489658733E+29 -0.489660103E+29 0.14E+24-23.4614086151 T T 95 -0.257165213E+30 -0.257165930E+30 0.72E+24-23.4142875671 T T 96 -0.136607590E+31 -0.136607968E+31 0.38E+25-23.1987590790 T T 97 -0.733873885E+31 -0.733875940E+31 0.21E+26-23.4918403625 T T 98 -0.398651824E+32 -0.398652985E+32 0.12E+27-24.4212017059 T T 99 -0.218944957E+33 -0.218945596E+33 0.64E+27-24.4561748505 T T ** -0.121559677E+34 -0.121560033E+34 0.36E+28-24.5605354309 T T YN for X = 475.779999 -- Epsilon = 1.19209290E-07 0 -0.210703555E-01 -0.210703555E-01 0.00E+00 -0.0000000000 T T 1 0.298792869E-01 0.298792869E-01 0.00E+00 0.0000000000 T T 2 0.211959574E-01 0.211959574E-01 0.00E+00 0.0000000000 T T 3 -0.297010876E-01 -0.297010876E-01 0.00E+00 -0.0000000000 T T 4 -0.215705149E-01 -0.215705149E-01 0.00E+00 -0.0000000000 T T 5 0.293383896E-01 0.293383896E-01 0.00E+00 0.0000000000 T T 6 0.221871529E-01 0.221871529E-01 0.00E+00 0.0000000000 T T 7 -0.287787914E-01 -0.287787914E-01 0.00E+00 -0.0000000000 T T 8 -0.230339784E-01 -0.230339784E-01 0.00E+00 -0.0000000000 T T 9 0.280041825E-01 0.280041825E-01 0.00E+00 0.0000000000 T T 10 0.240934491E-01 0.240934491E-01 0.00E+00 0.0000000000 T T 11 -0.269913841E-01 -0.269913841E-01 0.00E+00 -0.0000000000 T T 12 -0.253415275E-01 -0.253415275E-01 0.00E+00 -0.0000000000 T T 13 0.257130694E-01 0.257130694E-01 0.00E+00 0.0000000000 T T 14 0.267466716E-01 0.267466716E-01 0.00E+00 0.0000000000 T T 15 -0.241390076E-01 -0.241390076E-01 0.00E+00 -0.0000000000 T T 16 -0.282687414E-01 -0.282687414E-01 0.00E+00 -0.0000000000 T T 17 0.222377088E-01 0.222377088E-01 0.00E+00 0.0000000000 T T 18 0.298578832E-01 0.298578832E-01 0.00E+00 0.0000000000 T T 19 -0.199785046E-01 -0.199785046E-01 0.00E+00 -0.0000000000 T T 20 -0.314535424E-01 -0.314535424E-01 0.00E+00 -0.0000000000 T T 21 0.173341278E-01 0.173341278E-01 0.00E+00 0.0000000000 T T 22 0.329837315E-01 0.329837315E-01 0.00E+00 0.0000000000 T T 23 -0.142838014E-01 -0.142838014E-01 0.00E+00 -0.0000000000 T T 24 -0.343647376E-01 -0.343647376E-01 0.00E+00 -0.0000000000 T T 25 0.108168470E-01 0.108168479E-01 -0.93E-09 -0.7222529650 T T 26 0.355014876E-01 0.355014876E-01 0.00E+00 0.0000000000 T T 27 -0.693673966E-02 -0.693674106E-02 0.14E-08 -1.6893743277 T T 28 -0.362887941E-01 -0.362887941E-01 0.00E+00 -0.0000000000 T T 29 0.266549550E-02 0.266549736E-02 -0.19E-08 -5.8619494438 T T 30 0.366137326E-01 0.366137326E-01 0.00E+00 0.0000000000 T T 31 0.195181498E-02 0.195181323E-02 0.17E-08 7.5050339699 T T 32 -0.363593884E-01 -0.363593884E-01 0.00E+00 -0.0000000000 T T 33 -0.684273290E-02 -0.684273103E-02 -0.19E-08 2.2834444046 T T 34 0.354101658E-01 0.354101658E-01 0.00E+00 0.0000000000 T T 35 0.119036669E-01 0.119036660E-01 0.93E-09 0.6563103795 T T 36 -0.336588174E-01 -0.336588174E-01 0.00E+00 -0.0000000000 T T 37 -0.169972703E-01 -0.169972703E-01 0.00E+00 -0.0000000000 T T 38 0.310151633E-01 0.310151633E-01 0.00E+00 0.0000000000 T T 39 0.219515618E-01 0.219515599E-01 0.19E-08 0.7117944360 T T 40 -0.274163950E-01 -0.274163950E-01 0.00E+00 -0.0000000000 T T 41 -0.265614893E-01 -0.265614875E-01 -0.19E-08 0.5882576704 T T 42 0.228385609E-01 0.228385609E-01 0.00E+00 0.0000000000 T T 43 0.305936877E-01 0.305936858E-01 0.19E-08 0.5107262731 T T 44 -0.173085742E-01 -0.173085742E-01 0.00E+00 -0.0000000000 T T 45 -0.337950736E-01 -0.337950699E-01 -0.37E-08 0.9246909618 T T 46 0.109157944E-01 0.109157953E-01 -0.93E-09 -0.7157060504 T T 47 0.359058231E-01 0.359058194E-01 0.37E-08 0.8703323603 T T 48 -0.382186985E-02 -0.382187171E-02 0.19E-08 -4.0883131027 T T 49 -0.366769768E-01 -0.366769731E-01 -0.37E-08 0.8520331383 T T 50 -0.373276439E-02 -0.373276183E-02 -0.26E-08 5.7556204796 T T 51 0.358924195E-01 0.358924158E-01 0.37E-08 0.8706573844 T T 52 0.114275543E-01 0.114275496E-01 0.47E-08 3.4182729721 T T 53 -0.333944894E-01 -0.333944857E-01 -0.37E-08 0.9357831478 T T 54 -0.188675821E-01 -0.188675746E-01 -0.75E-08 3.3125600815 T T 55 0.291116294E-01 0.291116275E-01 0.19E-08 0.5367270708 T T 56 0.255981702E-01 0.255981628E-01 0.75E-08 2.4415807724 T T 57 -0.230857451E-01 -0.230857451E-01 0.00E+00 -0.0000000000 T T 58 -0.311296657E-01 -0.311296582E-01 -0.75E-08 2.0077311993 T T 59 0.154960155E-01 0.154960174E-01 -0.19E-08 -1.0083236694 T T 60 0.349728912E-01 0.349728838E-01 0.75E-08 1.7870984077 T T 61 -0.667524338E-02 -0.667524710E-02 0.37E-08 -4.6814770699 T T 62 -0.366845652E-01 -0.366845578E-01 -0.75E-08 1.7037137747 T T 63 -0.288565969E-02 -0.288565271E-02 -0.70E-08 20.3051490784 T T 64 0.359203592E-01 0.359203555E-01 0.37E-08 0.8699801564 T T 65 0.125493826E-01 0.125493743E-01 0.84E-08 5.6028652191 T T 66 -0.324914232E-01 -0.324914195E-01 -0.37E-08 0.9617922902 T T 67 -0.215637758E-01 -0.215637647E-01 -0.11E-07 4.3475689888 T T 68 0.264181420E-01 0.264181420E-01 0.00E+00 0.0000000000 T T 69 0.291153062E-01 0.291152969E-01 0.93E-08 2.6832964420 T T 70 -0.179732461E-01 -0.179732479E-01 0.19E-08 -0.8693476915 T T 71 -0.344039984E-01 -0.344039910E-01 -0.75E-08 1.8166493177 T T 72 0.770512270E-02 0.770512689E-02 -0.42E-08 -4.5627112389 T T 73 0.367360376E-01 0.367360301E-01 0.75E-08 1.7013266087 T T 74 0.356786419E-02 0.356785767E-02 0.65E-08 15.3277978897 T T 75 -0.356261879E-01 -0.356261842E-01 -0.37E-08 0.8771637082 T T 76 -0.147997951E-01 -0.147997877E-01 -0.75E-08 4.2230315208 T T 77 0.308980178E-01 0.308980159E-01 0.19E-08 0.5056958795 T T 78 0.248008352E-01 0.248008259E-01 0.93E-08 3.1500954628 T T 79 -0.227662548E-01 -0.227662548E-01 0.00E+00 -0.0000000000 T T 80 -0.323611945E-01 -0.323611870E-01 -0.75E-08 1.9313255548 T T 81 0.118835121E-01 0.118835149E-01 -0.28E-08 -1.9722704887 T T 82 0.364074521E-01 0.364074484E-01 0.37E-08 0.8583407402 T T 83 0.666032720E-03 0.666027889E-03 0.48E-08 60.8488769531 T T 84 -0.361750722E-01 -0.361750722E-01 0.00E+00 -0.0000000000 T T 85 -0.134396097E-01 -0.134396041E-01 -0.56E-08 3.4878244400 T T 86 0.313729905E-01 0.313729942E-01 -0.37E-08 -0.9960797429 T T 87 0.247813109E-01 0.247813053E-01 0.56E-08 1.8915463686 T T 88 -0.223100875E-01 -0.223100930E-01 0.56E-08 -2.1010675430 T T 89 -0.330342352E-01 -0.330342278E-01 -0.75E-08 1.8919765949 T T 90 0.995123666E-02 0.995124504E-02 -0.84E-08 -7.0657048225 T T 91 0.367990471E-01 0.367990434E-01 0.37E-08 0.8492067456 T T 92 0.412549311E-02 0.412548333E-02 0.98E-08 19.8839874268 T T 93 -0.352035798E-01 -0.352035798E-01 0.00E+00 -0.0000000000 T T 94 -0.178878754E-01 -0.178878643E-01 -0.11E-07 5.2409801483 T T 95 0.281353537E-01 0.281353574E-01 -0.37E-08 -1.1107022762 T T 96 0.291235670E-01 0.291235559E-01 0.11E-07 3.2190425396 T T 97 -0.163825992E-01 -0.163826086E-01 0.93E-08 -4.7687792778 T T 98 -0.358035974E-01 -0.358035900E-01 -0.75E-08 1.7456346750 T T 99 0.163312408E-02 0.163313653E-02 -0.12E-07-63.9830055237 T T ** 0.364832357E-01 0.364832357E-01 0.00E+00 0.0000000000 T T </pre> David Miller 2012-05-16T23:56:04 Re: 32b intel fortran vs. 64b linux gfortran http://permalink.gmane.org/gmane.comp.gcc.fortran/38816 <pre>All CPUs on which you are likely to run Fortran nowadays have IEEE754 compliant firmware for sqrt. tan(), however, is difficult enough that numerical standardization proposals haven't been implemented. You will find numerical differences even with the same compiler between scalar and vector libraries, for example, or possibly differing options for various instruction architecture targets. In your example, yes, the correct interpretation is that your decimal value is rounded to the closest representable value in single precision binary before passing to the print and tan functions. </pre> Tim Prince 2012-05-16T14:12:18 Re: 32b intel fortran vs. 64b linux gfortran http://permalink.gmane.org/gmane.comp.gcc.fortran/38815 <pre> On ma., 2012-05-14 at 14:31 +0100, N.M. Maclaren wrote: On ma., 2012-05-14 at 14:38 +0200, Tobias Burnus wrote: What Nick tell us here is probably correct with respect to what he refers, but the "error" I experience with my software seem quite odd (at least for me). I get different result in the 7th digit for the tan() function, above a certain level. When working with sqrt(), I experience no difference at all between the two compilers and platforms, pushing it really hard by taking the sqrt 15 times a = sqrt(r) b = sqrt(a) c = sqrt(b) . . . and going back again, b = c*c a = b*b r = a*a ending up with exactly the same error. What seem odd to me is that the varable "chi1" holds the same value in advance of the call to "tan(chi1)" for both ifort and gfortran compiled source, but produce a different result. Win7 ifort compiled source: chi1 = 0.52359876 print '(f15.10)', chi1 print '(f15.10)', tan(chi1) 0.5235987306 0.5773501992 Increasing one number in the last digit 6 to 7: chi1 = 0.52359877 print '(f15.10)', chi1 print '(f15.10)', tan(chi1) 0.5235987902 0.5773502588 ############################## While for Linux 64b gfortran compiled source, the last calculation produce: chi1 = 0.52359877 print '(f15.10)', chi1 print '(f15.10)', tan(chi1) 0.5235987902 0.5773503184 So, something different happens when working around 30 deg (pi/6). I'm far from an expert, but from earlier experience I remember problems with angles like pi/6 and pi/3 and so on. Cheers, K.E </pre> Edvardsen Kåre 2012-05-16T13:05:17 Re: [Patch, libgfortran] Pass mode in "open" for O_CREAT and on VxWorks http://permalink.gmane.org/gmane.comp.gcc.fortran/38814 <pre> Ah, bummer. We have something roughly similiar for snprintf (see libgfortran.h), but it seems that it works slightly differently due to using a variadic macro etc. So it seems this idea will not work, sorry. Interestingly, according to that document open() also has the POSIXly correct varargs prototype. So apparently the old 3-argument prototype was used only in some older vxworks version? Since the document above is from 2007, one wonders how long ago did the prototype change, and how much new development is done on these older vxworks versions? IMHO lets do only the second part now, and wait for someone with vxworks experience to pop in and explain what changes are necessary, and which vxworks versions are worth considering to support. </pre> Janne Blomqvist 2012-05-16T11:26:13 Re: [Patch, libgfortran] Pass mode in "open" for O_CREAT and on VxWorks http://permalink.gmane.org/gmane.comp.gcc.fortran/38812 <pre> IMHO it would be cleaner if you instead somewhere in the beginning of unix.c did #ifdef __VXWORKS__ /* open is not a variadic function on vxworks (or something...) */ #define open(path, flags) open(path, flags, 0666) #endif Untested, but AFAICS it should work (unless I'm missing something wrt macro expansion, which of course is quite possible). Good catch. This part is ok as is. Ok with the above change. I suspect there isn't that much demand for gfortran on vxworks, but the patch isn't intrusive either so if you want to backport it go ahead. </pre> Janne Blomqvist 2012-05-16T06:45:08 [Fortran, patch] PR 48831 - Constant expression (PARAMETER array element) rejected as nonconstant http://permalink.gmane.org/gmane.comp.gcc.fortran/38811 <pre>Hi All, in attachment there's a patch for PR 48831, it also includes a new test case suggested by Tobias Burnus. The patch is bootstrapped and tested on x86_64-unknown-linux-gnu. Regards. Alessandro Fanfarillo </pre> Alessandro Fanfarillo 2012-05-15T20:39:35 Search Engine Search the mailing list at Gmane query http://search.gmane.org/?group=$group=gmane.comp.gcc.fortran