Sqrt2f21l.cin

// Sqrt2f21l.cin is code for evaluation of arctetration to base $\sqrt{2}$ of complex argument.

//

//z_type tq2L(z_type z){ int n; z_type e,s,k;
z_type f21L(z_type z){ int n; z_type e,s,k;
DB TcL[23]={1., //coeff. of expansion of exp(-q(z+1.2 ...) by powers of (2-F).
 -.56472283831773236365, 0.29964618138408807683,
 -.15593239048925425850, 0.8035187974815443609e-1,
 -0.411584960662439279e-1, 0.2099852095441203541e-1,
 -0.1068258032026355653e-1, 0.542288102231591005e-2,
 -0.2748252661868267e-2, 0.13909151872677962e-2,
 -0.703181586212482131e-3, 0.35517006776480e-3,
 -0.1792537427481520668e-3, 0.9040887657183e-4,
 -0.45572543028501136e-4, 0.2296022632181e-4,
 -0.1156277075032e-4, 0.5820169657e-5,
  -.2928968839277e-5, 0.144e-5, -.71e-6
};
 z=2.-z; s=TcL[22];
         for(n=21; n>=0; n--){ s*=z; s+=TcL[n]; }
        return -log(s*z)/0.36651292058166432701 -1.251551478822190;};
                        //.32663425997828098238;
z_type F21L(z_type z){ DB b=sqrt(2.);
                        if(abs(z-2.)>9999.) return 9999.;
                        if(abs(z-2.)>.4) return F21L(exp(z*log(b)))-1. ;
                         return f21L(z);
                        }
//

//