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, search

// 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);
                        }

//