Difference between revisions of "File:Sqrt2tetatemap.jpg"

From TORI
Jump to: navigation, search
(Importing image file)
 
 
Line 1: Line 1:
  +
Range of validity of relation
−
Importing image file
  
  +
  +
$\mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z))=z$
  +
  +
for the [[tetration]]
  +
and [[arctetration]] to base $\sqrt{2}$
  +
in the plane $x=\Re(z)$, $y=\Im(z)$
  +
is shaded with lines
  +
  +
$u=\Re \Big( \mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z)) \Big) = \mathrm{const}$
  +
  +
and lines
  +
  +
$v=\Im \Big( \mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z)) \Big) = \mathrm{const}$
  +
  +
These lines provide the shading of the central part of the picture.
  +
  +
In addition, lines
  +
  +
$\Im\Big( \exp_{\sqrt{2}}^n (x\!+\!\mathrm i y)\Big)=\pm \frac{|P|}{2}=\pm \frac{2\, \pi}{\ln(2)}
  +
$
  +
  +
are drawn for $n=0,1,2,3,4$ as bounds if the shaded range.
  +
  +
Usage: this is figure 16.5 of the book [[Суперфункции]] (2014, In Russian) <ref>
  +
https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0 <br>
  +
http://www.ils.uec.ac.jp/~dima/BOOK/202.pdf <br>
  +
http://mizugadro.mydns.jp/BOOK/202.pdf
  +
Д.Кузнецов. Суперфункции. [[Lambert Academic Publishing]], 2014.
  +
</ref>; the English version is in preparation in 2015. (Numeration of figures in the English version may be different from that of the Russian version.)
  +
  +
This image is used also in figure 2
  +
The algorithm of the evaluation is also described in the article
  +
<ref>
  +
http://www.ams.org/journals/mcom/2010-79-271/S0025-5718-10-02342-2/home.html <br>
  +
http://mizugadro.mydns.jp/PAPERS/2010sqrt2.pdf offprint
  +
D.Kouznetsov, H.Trappmann. Portrait of the four regular super-exponentials to base sqrt(2). Mathematics of Computation, 2010, v.79, p.1727-1756.
  +
</ref>.
  +
(top right map)
  +
  +
==Refereces==
  +
<references/>
  +
  +
==[[C++]] generator of the shading==
  +
Files [[ado.cin]],
  +
[[conto.cin]],
  +
[[sqrt2f21e.cin]]
  +
[[sqrt2f21l.cin]]
  +
should be loaded in order to compile the code below.
  +
  +
<poem><nomathjax><nowiki>
  +
#include <math.h>
  +
#include <stdio.h>
  +
#include <stdlib.h>
  +
#define DB double
  +
#define DO(x,y) for(x=0;x<y;x++)
  +
#include <complex>
  +
typedef std::complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
#include "conto.cin"
  +
// #include "tq2e.cin"
  +
// #include "tq2L.cin"
  +
#include "sqrt2f21e.cin"
  +
#include "sqrt2f21l.cin"
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
  +
int M=401,M1=M+1;
  +
int N=501,N1=N+1;
  +
DB X[M1],Y[N1], g[M1*N1],f[M1*N1], w[M1*N1]; // w is working array.
  +
char v[M1*N1]; // v is working array
  +
FILE *o;o=fopen("sqrt2tetatema.eps","w"); ado(o,0,0,214,212);
  +
fprintf(o,"112 110 translate\n 10 10 scale\n");
  +
// DB sy=10.1/sinh(N/2./100.);
  +
DO(m,M1) X[m]=-10+.05*(m-.5);
  +
DO(n,N1) Y[n]=-10+.04*(n-.5);
  +
// DO(n,N1) Y[n]=sy*sinh((n-N/2.+.5)/100.);
  +
for(m=-10;m<11;m++) {M(m,-10)L(m,10)}
  +
for(n=-10;n<11;n++) {M( -10,n)L(10,n)} fprintf(o,".006 W 0 0 0 RGB S\n");
  +
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
  +
f[m*N1+n]=9999;}
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
// if( abs(z-2.)>.2 || abs(z-4.)>.2)
  +
{ c=F21L(z);
  +
c=F21E(c);
  +
if(abs(c-z)<.1) {
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}
  +
}}}
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
//p=2; q=1.1;
  +
p=.1;q=.1;
  +
for(m=-10;m<10;m++) for(n=2;n<10;n+=2)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q, q);
  +
fprintf(o,".01 W 0 .6 0 RGB S\n");
  +
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q, q);
  +
fprintf(o,".01 W .9 0 0 RGB S\n");
  +
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q, q);
  +
fprintf(o,".01 W 0 0 .9 RGB S\n");
  +
  +
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".06 W .9 0 0 RGB S\n");
  +
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".06 W 0 0 .9 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
  +
for(m=-10;m<11;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");
  +
  +
  +
//#include "plofu.cin"
  +
  +
// M(-13.2,0)L(2,0) fprintf(o,".05 W 0 .8 0 RGB S\n");
  +
  +
fprintf(o,"0 setlinejoin 0 setlinecap\n");
  +
  +
M(2,0)L(10.1,0)fprintf(o,".05 W 1 1 1 RGB S\n");
  +
DO(n,27){M(2+.3*n,0)L(2+.3*(n+.5) ,0)} fprintf(o,".1 W 0 0 0 RGB S\n");
  +
  +
//M(2,0)L(10.1,0)fprintf(o,".1 W 0 0 0 RGB [.19 .19] 0 setdash S\n");
  +
//M(-10,0)L(-2,0)fprintf(o,".04 W 1 0 1 RGB S\n"); // fails at some printers
  +
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
  +
system("epstopdf sqrt2tetatema.eps");
  +
system( "open sqrt2tetatema.pdf"); //for linux
  +
// getchar(); system("killall Preview"); // For macintosh
  +
}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[C++]] generator boundaries of the shading==
  +
Files [[ado.cin]],
  +
[[conto.cin]],
  +
[[sqrt2f21e.cin]]
  +
should be loaded in order to compile the code below.
  +
  +
<poem><nomathjax><nowiki>
  +
  +
#include <math.h>
  +
#include <stdio.h>
  +
#include <stdlib.h>
  +
#define DB double
  +
#define DO(x,y) for(x=0;x<y;x++)
  +
#include <complex>
  +
typedef std::complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
#include "conto.cin"
  +
// #include "tq2e.cin"
  +
// #include "tq2L.cin"
  +
#include "sqrt2f21e.cin"
  +
#include "sqrt2f21l.cin"
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
  +
int M=201,M1=M+1;
  +
int N=401,N1=N+1;
  +
DB X[M1],Y[N1], g[M1*N1],f[M1*N1], w[M1*N1]; // w is working array.
  +
char v[M1*N1]; // v is working array
  +
FILE *o;o=fopen("sqrt2tetatemb.eps","w"); ado(o,0,0,214,212);
  +
fprintf(o,"112 110 translate\n 10 10 scale\n");
  +
// DB sy=10.1/sinh(N/2./100.);
  +
  +
DO(m,M1) X[m]=0+.05*(m-.5);
  +
DO(n,N1) Y[n]=-10+.05*(n-.5);
  +
for(m=-10;m<11;m++) {M(m,-10)L(m,10)}
  +
for(n=-10;n<11;n++) {M( -10,n)L(10,n)} fprintf(o,".006 W 0 0 0 RGB S\n");
  +
  +
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
  +
f[m*N1+n]=9999;}
  +
  +
M(-10, 2./log(2.)*M_PI) L( 10, 2./log(2.)*M_PI)
  +
M(-10,-2./log(2.)*M_PI) L( 10,-2./log(2.)*M_PI)
  +
  +
p=96.;
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=2;q=1;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
  +
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
c=exp(.5*log(2.)*c);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=96.;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 0 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 0 RGB S\n");
  +
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=96.;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 .6 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 .6 RGB S\n");
  +
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=96.;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 .5 0 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 .5 0 RGB S\n");
  +
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=96.;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 0 0 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");
  +
  +
DO(m,M1){x=X[m];
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
c=exp(.5*log(2.)*z);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
c=exp(.5*log(2.)*c);
  +
p=Re(c); q=Im(c);
  +
if(p>-99 && p<99) g[m*N1+n]=p;
  +
if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
  +
}}
  +
p=96.;
  +
conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 0 0 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");
  +
  +
/*
  +
p=4;q=1;
  +
for(m=-10;m<10;m++) for(n=2;n<10;n+=2)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q, q);
  +
fprintf(o,".01 W 0 .6 0 RGB S\n");
  +
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q, q);
  +
fprintf(o,".01 W .9 0 0 RGB S\n");
  +
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q, q);
  +
fprintf(o,".01 W 0 0 .9 RGB S\n");
  +
  +
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".06 W .9 0 0 RGB S\n");
  +
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".06 W 0 0 .9 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
  +
for(m=-10;m<11;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");
  +
*/
  +
  +
// #include "plofu.cin"
  +
  +
// M(-13.2,0)L(2,0) fprintf(o,".05 W 0 .8 0 RGB S\n");
  +
  +
//M(2,0)L(10.1,0)fprintf(o,".05 W 1 1 1 RGB S\n");
  +
//DO(n,27){M(2+.3*n,0)L(2+.3*(n+.5) ,0)} fprintf(o,".1 W 0 0 0 RGB S\n");
  +
  +
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
  +
system("epstopdf sqrt2tetatemb.eps");
  +
system( "open sqrt2tetatemb.pdf"); //for linux
  +
// getchar(); system("killall Preview"); // For macintosh
  +
}
  +
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] combiner==
  +
Files generated with codes above
  +
should be loaded in order to compile the code below.
  +
  +
<poem><nomathjax><nowiki>
  +
\documentclass[12pt]{article}
  +
\paperwidth 422px
  +
\paperheight 418px
  +
\textwidth 1394px
  +
\textheight 1300px
  +
\topmargin -94px
  +
\oddsidemargin -76px
  +
\usepackage{graphics}
  +
\usepackage{rotating}
  +
\newcommand \sx {\scalebox}
  +
\newcommand \rot {\begin{rotate}}
  +
\newcommand \ero {\end{rotate}}
  +
\newcommand \ing {\includegraphics}
  +
\newcommand \rmi {\mathrm{i}}
  +
\parindent 0pt
  +
\pagestyle{empty}
  +
\begin{document}\parindent 0pt
  +
  +
\sx{2}{\begin{picture}(204,204)
  +
\put(0,0){\ing{sqrt2tetatema}}
  +
\put(0,0){\ing{sqrt2tetatemb}}
  +
\put(6,206){\sx{.8}{$y$}}
  +
\put(6,188){\sx{.8}{$8$}}
  +
\put(6,168){\sx{.8}{$6$}}
  +
\put(6,148){\sx{.8}{$4$}}
  +
\put(6,128){\sx{.8}{$2$}}
  +
\put(6,108){\sx{.8}{$0$}}
  +
\put(-1, 88){\sx{.8}{$-2$}}
  +
\put(-1, 68){\sx{.8}{$-4$}}
  +
\put(-1, 48){\sx{.8}{$-6$}}
  +
\put(-1, 28){\sx{.8}{$-8$}}
  +
\put(24,2){\sx{.8}{$-8$}}
  +
\put(44,2){\sx{.8}{$-6$}}
  +
\put(64,2){\sx{.8}{$-4$}}
  +
\put(84,2){\sx{.8}{$-2$}}
  +
\put(110.5,2){\sx{.8}{$0$}}
  +
\put(130.5,2){\sx{.8}{$2$}}
  +
\put(150.5,2){\sx{.8}{$4$}}
  +
\put(170.5,2){\sx{.8}{$6$}}
  +
\put(190.5,2){\sx{.8}{$8$}}
  +
\put(208.6,2){\sx{.8}{$x$}}
  +
  +
\put(26,107.4){$v\!=\!0$}
  +
%\put(176,108){\bf cut}
  +
\put(115.4,57.4){\rot{90}$u\!=\!0$\ero}
  +
\end{picture}}
  +
\end{document}
  +
</nowiki></nomathjax></poem>
  +
  +
  +
[[Category:Complex map]]
  +
[[Category:Agreement]]
  +
[[Category:Base sqrt2]]
  +
[[Category:Book]]
  +
[[Category:BookMap]]
  +
[[Category:Tetration]]
  +
[[Category:Arctetration]]
  +
[[Category:Inverse function]]
  +
[[Category:C++]]
  +
[[Category:Latex]]

Latest revision as of 08:52, 1 December 2018

Range of validity of relation

$\mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z))=z$

for the tetration and arctetration to base $\sqrt{2}$ in the plane $x=\Re(z)$, $y=\Im(z)$ is shaded with lines

$u=\Re \Big( \mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z)) \Big) = \mathrm{const}$

and lines

$v=\Im \Big( \mathrm{tet}_{\sqrt{2}}(\mathrm{ate}_{\sqrt{2}}(z)) \Big) = \mathrm{const}$

These lines provide the shading of the central part of the picture.

In addition, lines

$\Im\Big( \exp_{\sqrt{2}}^n (x\!+\!\mathrm i y)\Big)=\pm \frac{|P|}{2}=\pm \frac{2\, \pi}{\ln(2)} $

are drawn for $n=0,1,2,3,4$ as bounds if the shaded range.

Usage: this is figure 16.5 of the book Суперфункции (2014, In Russian) [1]; the English version is in preparation in 2015. (Numeration of figures in the English version may be different from that of the Russian version.)

This image is used also in figure 2 The algorithm of the evaluation is also described in the article [2]. (top right map)

Refereces

  1. ↑ https://www.morebooks.de/store/ru/book/Суперфункции/isbn/978-3-659-56202-0
    http://www.ils.uec.ac.jp/~dima/BOOK/202.pdf
    http://mizugadro.mydns.jp/BOOK/202.pdf Д.Кузнецов. Суперфункции. Lambert Academic Publishing, 2014.
  2. ↑ http://www.ams.org/journals/mcom/2010-79-271/S0025-5718-10-02342-2/home.html
    http://mizugadro.mydns.jp/PAPERS/2010sqrt2.pdf offprint D.Kouznetsov, H.Trappmann. Portrait of the four regular super-exponentials to base sqrt(2). Mathematics of Computation, 2010, v.79, p.1727-1756.

C++ generator of the shading

Files ado.cin, conto.cin, sqrt2f21e.cin sqrt2f21l.cin should be loaded in order to compile the code below.


 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
 #define DB double
 #define DO(x,y) for(x=0;x<y;x++)
 #include <complex>
 typedef std::complex<double> z_type;
 #define Re(x) x.real()
 #define Im(x) x.imag()
 #define I z_type(0.,1.)
 #include "conto.cin"
// #include "tq2e.cin"
// #include "tq2L.cin"
#include "sqrt2f21e.cin"
#include "sqrt2f21l.cin"
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
 int M=401,M1=M+1;
 int N=501,N1=N+1;
 DB X[M1],Y[N1], g[M1*N1],f[M1*N1], w[M1*N1]; // w is working array.
 char v[M1*N1]; // v is working array
 FILE *o;o=fopen("sqrt2tetatema.eps","w"); ado(o,0,0,214,212);
 fprintf(o,"112 110 translate\n 10 10 scale\n");
// DB sy=10.1/sinh(N/2./100.);
 DO(m,M1) X[m]=-10+.05*(m-.5);
 DO(n,N1) Y[n]=-10+.04*(n-.5);
// DO(n,N1) Y[n]=sy*sinh((n-N/2.+.5)/100.);
 for(m=-10;m<11;m++) {M(m,-10)L(m,10)}
 for(n=-10;n<11;n++) {M( -10,n)L(10,n)} fprintf(o,".006 W 0 0 0 RGB S\n");
 DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
                        f[m*N1+n]=9999;}
 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
 // if( abs(z-2.)>.2 || abs(z-4.)>.2)
        { c=F21L(z);
                c=F21E(c);
                if(abs(c-z)<.1) {
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
                }
        }}}
fprintf(o,"1 setlinejoin 2 setlinecap\n");
//p=2; q=1.1;
p=.1;q=.1;
for(m=-10;m<10;m++) for(n=2;n<10;n+=2)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q, q);
                                                        fprintf(o,".01 W 0 .6 0 RGB S\n");
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q, q);
                                                        fprintf(o,".01 W .9 0 0 RGB S\n");
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q, q);
                                                        fprintf(o,".01 W 0 0 .9 RGB S\n");

for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".06 W .9 0 0 RGB S\n");
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".06 W 0 0 .9 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (0. ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
for(m=-10;m<11;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");


//#include "plofu.cin"

// M(-13.2,0)L(2,0) fprintf(o,".05 W 0 .8 0 RGB S\n");

fprintf(o,"0 setlinejoin 0 setlinecap\n");

M(2,0)L(10.1,0)fprintf(o,".05 W 1 1 1 RGB S\n");
DO(n,27){M(2+.3*n,0)L(2+.3*(n+.5) ,0)} fprintf(o,".1 W 0 0 0 RGB S\n");

//M(2,0)L(10.1,0)fprintf(o,".1 W 0 0 0 RGB [.19 .19] 0 setdash S\n");
//M(-10,0)L(-2,0)fprintf(o,".04 W 1 0 1 RGB S\n"); // fails at some printers
 fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
        system("epstopdf sqrt2tetatema.eps");
        system( "open sqrt2tetatema.pdf"); //for linux
// getchar(); system("killall Preview"); // For macintosh
}

C++ generator boundaries of the shading

Files ado.cin, conto.cin, sqrt2f21e.cin should be loaded in order to compile the code below.



 #include <math.h>
 #include <stdio.h>
 #include <stdlib.h>
 #define DB double
 #define DO(x,y) for(x=0;x<y;x++)
 #include <complex>
 typedef std::complex<double> z_type;
 #define Re(x) x.real()
 #define Im(x) x.imag()
 #define I z_type(0.,1.)
 #include "conto.cin"
// #include "tq2e.cin"
// #include "tq2L.cin"
#include "sqrt2f21e.cin"
#include "sqrt2f21l.cin"
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
 int M=201,M1=M+1;
 int N=401,N1=N+1;
 DB X[M1],Y[N1], g[M1*N1],f[M1*N1], w[M1*N1]; // w is working array.
 char v[M1*N1]; // v is working array
 FILE *o;o=fopen("sqrt2tetatemb.eps","w"); ado(o,0,0,214,212);
 fprintf(o,"112 110 translate\n 10 10 scale\n");
// DB sy=10.1/sinh(N/2./100.);

 DO(m,M1) X[m]=0+.05*(m-.5);
 DO(n,N1) Y[n]=-10+.05*(n-.5);
 for(m=-10;m<11;m++) {M(m,-10)L(m,10)}
 for(n=-10;n<11;n++) {M( -10,n)L(10,n)} fprintf(o,".006 W 0 0 0 RGB S\n");

 DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
                        f[m*N1+n]=9999;}

M(-10, 2./log(2.)*M_PI) L( 10, 2./log(2.)*M_PI)
M(-10,-2./log(2.)*M_PI) L( 10,-2./log(2.)*M_PI)

p=96.;
 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=2;q=1;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");

 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                c=exp(.5*log(2.)*c);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=96.;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 0 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .6 0 0 RGB S\n");

 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=96.;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 .6 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 .6 RGB S\n");

 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=96.;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 .5 0 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 .5 0 RGB S\n");

 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=96.;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 0 0 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");

 DO(m,M1){x=X[m];
 DO(n,N1){y=Y[n]; z=z_type(x,y);
                c=exp(.5*log(2.)*z);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                c=exp(.5*log(2.)*c);
                p=Re(c); q=Im(c);
                if(p>-99 && p<99) g[m*N1+n]=p;
                if(q>-99 && q<99 && fabs(q)>1.e-18) f[m*N1+n]=q;
        }}
p=96.;
                  conto(o,f,w,v,X,Y,M,N, (2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W .0 0 0 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (-2./log(2.)*M_PI ),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");

/*
p=4;q=1;
for(m=-10;m<10;m++) for(n=2;n<10;n+=2)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q, q);
                                                        fprintf(o,".01 W 0 .6 0 RGB S\n");
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q, q);
                                                        fprintf(o,".01 W .9 0 0 RGB S\n");
for(m=0;m<10;m++) for(n=2;n<10;n+=2)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q, q);
                                                        fprintf(o,".01 W 0 0 .9 RGB S\n");

for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".06 W .9 0 0 RGB S\n");
for(m=1;m<11;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".06 W 0 0 .9 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (0. ),-p,p);fprintf(o,".06 W .6 0 .6 RGB S\n");
for(m=-10;m<11;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".06 W 0 0 0 RGB S\n");
*/

// #include "plofu.cin"

// M(-13.2,0)L(2,0) fprintf(o,".05 W 0 .8 0 RGB S\n");

//M(2,0)L(10.1,0)fprintf(o,".05 W 1 1 1 RGB S\n");
//DO(n,27){M(2+.3*n,0)L(2+.3*(n+.5) ,0)} fprintf(o,".1 W 0 0 0 RGB S\n");

 fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
        system("epstopdf sqrt2tetatemb.eps");
        system( "open sqrt2tetatemb.pdf"); //for linux
// getchar(); system("killall Preview"); // For macintosh
}

Latex combiner

Files generated with codes above should be loaded in order to compile the code below.


\documentclass[12pt]{article}
\paperwidth 422px
\paperheight 418px
\textwidth 1394px
\textheight 1300px
\topmargin -94px
\oddsidemargin -76px
\usepackage{graphics}
\usepackage{rotating}
\newcommand \sx {\scalebox}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\newcommand \ing {\includegraphics}
\newcommand \rmi {\mathrm{i}}
\parindent 0pt
\pagestyle{empty}
\begin{document}\parindent 0pt

\sx{2}{\begin{picture}(204,204)
\put(0,0){\ing{sqrt2tetatema}}
\put(0,0){\ing{sqrt2tetatemb}}
\put(6,206){\sx{.8}{$y$}}
\put(6,188){\sx{.8}{$8$}}
\put(6,168){\sx{.8}{$6$}}
\put(6,148){\sx{.8}{$4$}}
\put(6,128){\sx{.8}{$2$}}
\put(6,108){\sx{.8}{$0$}}
\put(-1, 88){\sx{.8}{$-2$}}
\put(-1, 68){\sx{.8}{$-4$}}
\put(-1, 48){\sx{.8}{$-6$}}
\put(-1, 28){\sx{.8}{$-8$}}
\put(24,2){\sx{.8}{$-8$}}
\put(44,2){\sx{.8}{$-6$}}
\put(64,2){\sx{.8}{$-4$}}
\put(84,2){\sx{.8}{$-2$}}
\put(110.5,2){\sx{.8}{$0$}}
\put(130.5,2){\sx{.8}{$2$}}
\put(150.5,2){\sx{.8}{$4$}}
\put(170.5,2){\sx{.8}{$6$}}
\put(190.5,2){\sx{.8}{$8$}}
\put(208.6,2){\sx{.8}{$x$}}

\put(26,107.4){$v\!=\!0$}
%\put(176,108){\bf cut}
\put(115.4,57.4){\rot{90}$u\!=\!0$\ero}
\end{picture}}
\end{document}

File history

Click on a date/time to view the file as it appeared at that time.

Date/TimeThumbnailDimensionsUserComment
current06:14, 1 December 2018Thumbnail for version as of 06:14, 1 December 20181,758 × 1,741 (1,008 KB)Maintenance script (talk | contribs)Importing image file
  • You cannot overwrite this file.

There are no pages that link to this file.

Metadata

This file contains additional information, probably added from the digital camera or scanner used to create or digitize it. If the file has been modified from its original state, some details may not fully reflect the modified file.
Retrieved from "https://mizugadro.mydns.jp/t/index.php?title=File:Sqrt2tetatemap.jpg&oldid=21607"