Sqrt2f43l.cin

From TORI
Revision as of 18:46, 30 July 2019 by T (talk | contribs) (Text replacement - "\$([^\$]+)\$" to "\\(\1\\)")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
The printable version is no longer supported and may have rendering errors. Please update your browser bookmarks and please use the default browser print function instead.

// Sqrt2f43e.cin is routine for evaluarion of abelexponent to base \(\sqrt{2}\) that approach 4 at \(-\infty\) and has value zero at 3.



z_type f43L(z_type z){ int n; z_type e,s;
DB Uco[21]={1,
 -.44858743119526122890, .21208912005491969757,
 -.10218436750697167872, 0.49698683037371830337e-1,
 -0.2430759032611955221e-1, 0.11933088396510860210e-1,
 -0.587369764200886206e-2, 0.289686728710575713e-2,
 -0.1430908106079253664e-2, 0.7076637148565759223e-3,
 -0.3503296158729878e-3, 0.17357560046634138e-3,
 -0.86061011929162626e-4, 0.426959089012891e-4,
 -0.2119302906819844809e-4, 0.1052442259960e-4,
 -0.52285174354086e-5, 0.259844999161e-5,
 -0.129178211214818578e-5, 0.4e-6 };
 z-=4.;
 s=Uco[19]; for(n=18; n>=0; n--){ s*=z; s+=Uco[n]; }
 return log(-s*z)/.32663425997828098238-1.90057764535874;
  }
z_type F43L(z_type z){ DB b=sqrt(2.);
                        if(abs(z-4.)>.4) return F43L(log(z)/log(b))+1. ;
                         return f43L(z);
                        }

//