Difference between revisions of "File:QexpMapT400.jpg"

From TORI
Jump to navigation Jump to search
 
(5 intermediate revisions by the same user not shown)
Line 10: Line 10:
 
==C++ generator of curves==
  
==C++ generator of curves==
   
  +
// FIles [[ado.cin]] and [[conto.cin]] and [[superexp.cin]] and [[superlo.cin]] should be loaded in the working directory in order to compile the code below. Actually, the last two evaluate [[tetration]] tet and [[arctetration]] ate; routines [[fsexp.cin]] and [[fslog.cin]] can be used instead.
−
==Latex generator of labels==
  
   
  +
#include <math.h>
  +
#include <stdio.h>
  +
#include <stdlib.h>
  +
#define DB double
  +
#define DO(x,y) for(x=0;x<y;x++)
  +
using namespace std;
  +
#include <complex>
  +
typedef complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
//#include "fac.cin"
  +
//#include "sinc.cin"
  +
//#include "facp.cin"
  +
//#include "afacc.cin"
  +
//#include "superfac.cin"
  +
#include "superex.cin"
  +
#include "superlo.cin"
  +
DB xL=0.31813150520476413;
  +
DB yL=1.3372357014306895;
  +
  +
#include "conto.cin"
  +
  +
int main(){ int j,k,m,n,n1; DB x,y, p,q, t; z_type z,c,d;
  +
int M=160,M1=M+1;
  +
int N=165,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("fig4b.eps","w");ado(o,402,402);
  +
FILE *o;o=fopen("qexpMap.eps","w");ado(o,402,402);
  +
fprintf(o,"201 201 translate\n 20 20 scale\n");
  +
DO(m,M1) X[m]=-8.+.1*m;
  +
DO(n,N1){ y=-8.+.1*n; if(y< -yL) Y[n]=y; else break;}
  +
Y[n]=-yL-.001; n++;
  +
Y[n]=-yL+.001; n++;
  +
for(j=n;j<N1;j++){y=-8.+.1*(j-2); if(y<-.02) Y[j]=y; else break;}
  +
Y[j]= -.02; j++;
  +
Y[j]= +.02; j++;
  +
for(k=j;k<N1;k++){y=-8.+.1*(k-3); if(y<yL) Y[k]=y; else break;}
  +
Y[k]= yL-.001; k++;
  +
Y[k]= yL+.001; k++;
  +
for(n=k;n<N1;n++){y=-8.+.1*(n-5); Y[n]=y;}
  +
for(m=-8;m<9;m++){if(m==0){M(m,-8.5)L(m,8.5)} else{M(m,-8)L(m,8)}}
  +
for(n=-8;n<9;n++){ M( -8,n)L(8,n)}
  +
fprintf(o,".008 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]; //printf("%5.2f\n",x);
  +
DO(n,N1){y=Y[n]; z=z_type(x,y);
  +
// c=afacc(z);
  +
// c=fac(c);
  +
// c=arcsuperfac(z);
  +
// c=superfac(.5+c);
  +
c=FSLOG(z);
  +
c=FSEXP(.5+c);
  +
// d=z;
  +
// p=abs(c-d)/(abs(c)+abs(d)); p=-log(p)/log(10.)-1.;
  +
p=Re(c);q=Im(c);
  +
if(p>-999 && p<999) g[m*N1+n]=p;
  +
if(q>-999 && q<999 && fabs(q)> 1.e-14) f[m*N1+n]=q;
  +
}}
  +
fprintf(o,"1 setlinejoin 1 setlinecap\n"); p=6.;q=6.;
  +
//#include"plofu.cin"
  +
//p=2;q=1;
  +
for(m=-4;m<4;m++) for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M-20,N,(m+.1*n),-q, q); fprintf(o,".02 W 0 .6 0 RGB S\n");
  +
for(m= 0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M-20,N,-(m+.1*n),-q, q); fprintf(o,".02 W .9 0 0 RGB S\n");
  +
for(m= 0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M-20,N, (m+.1*n),-q, q); fprintf(o,".02 W 0 0 .9 RGB S\n");
  +
for(m=1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p); fprintf(o,".05 W .9 0 0 RGB S\n");
  +
for(m=1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".05 W 0 0 .9 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".05 W .6 0 .6 RGB S\n");
  +
for(m=-8;m<9;m++) conto(o,g,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".05 W 0 0 0 RGB S\n");
  +
fprintf(o,"0 setlinejoin 0 setlinecap\n");
  +
M(xL, yL)L(-8.1, yL)
  +
M(xL,-yL)L(-8.1,-yL) fprintf(o," .06 W 1 1 1 RGB S\n");
  +
DO(m,21){M(xL-.4*m, yL)L(xL-.4*(m+.4), yL)}
  +
DO(m,21){M(xL-.4*m,-yL)L(xL-.4*(m+.4),-yL)} fprintf(o,".1 W 0 0 0 RGB S\n");
  +
//M(xL, yL)L(-8, yL)
  +
//M(xL,-yL)L(-8,-yL) fprintf(o,"[.2 .2]0 setdash .12 W 0 0 0 RGB S\n");
  +
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
  +
// system( "open fig4b.eps"); //for macintosh
  +
system("epstopdf qexpMap.eps");
  +
system( "open qexpMap.pdf"); //for unix
  +
// system( "xpdf fig4b.pdf &"); //for unix
  +
// getchar(); system("killall Preview"); //for macintosh
  +
}
  +
  +
==Latex generator of labels==
  +
<poem>
  +
\documentclass[12pt]{article}
  +
\usepackage{graphicx}
  +
\usepackage{geometry}
  +
\paperwidth 340pt
  +
\paperheight 340pt
  +
\usepackage{rotating}
  +
\textwidth 800mm
  +
\textheight 400mm
  +
\topmargin -61mm
  +
\oddsidemargin -77pt
  +
\parindent 0pt
  +
\newcommand \sx {\scalebox}
  +
\newcommand \rmi {{\rm i}} %imaginary unity is always roman font
  +
\newcommand \rot {\begin{rotate}}
  +
\newcommand \ero {\end{rotate}}
  +
\newcommand \ax {
  +
\put( 10,346){\sx{1.4}{$y$}}
  +
\put( 10,307){\sx{1.3}{$6$}}
  +
\put( 10,267){\sx{1.3}{$4$}}
  +
\put( 10,227){\sx{1.3}{$2$}}
  +
\put( 10,187){\sx{1.3}{$0$}}
  +
\put( 0,147){\sx{1.3}{$-2$}}
  +
\put( 0,107){\sx{1.3}{$-4$}}
  +
\put( 0, 67){\sx{1.3}{$-6$}}
  +
\put( 0, 27){\sx{1.3}{$-8$}}
  +
\put( 50, 17){\sx{1.3}{$-6$}}
  +
\put( 90, 17){\sx{1.3}{$-4$}}
  +
\put(130, 17){\sx{1.3}{$-2$}}
  +
\put(178, 17){\sx{1.3}{$0$}}
  +
\put(218, 17){\sx{1.3}{$2$}}
  +
\put(258, 17){\sx{1.3}{$4$}}
  +
\put(298, 17){\sx{1.3}{$6$}}
  +
\put(334, 17){\sx{1.4}{$x$}}
  +
}
  +
\begin {document}
  +
{\begin{picture}(420,420)
  +
%\put(-20,-10){\includegraphics{fig4b}}\ax
  +
\put(-20,-10){\includegraphics{qexpMap}}\ax
  +
\put( 52,188){\sx{1.1}{$v\!=\!0$}}
  +
\put( 81,176){\sx{1.1}{$u\!=\!-1.6$}}
  +
{\put(40,213){\sx{1.8}{\bf cut}}}
  +
{\put(40,160){\sx{1.8}{\bf cut}}}
  +
\put(300,318){\rot{19}\sx{1.3}{$u\!=\!-8$}\ero}
  +
\put(310,266){\rot{28}\sx{1.3}{$u\!= 8$}\ero}
  +
\put(300,218){\rot{-8}\sx{1.3}{$v\!= 8$}\ero}
  +
\put(300,156){\rot{ 8}\sx{1.3}{$v\!=\!-8$}\ero}
  +
\put(308,108){\rot{-30}\sx{1.3}{$u\!= 8$}\ero}
  +
\put(300, 56){\rot{-22}\sx{1.3}{$u\!=\!-8$}\ero}
  +
\end{picture}}
  +
\end{document}
  +
$</poem>
   
 
==References==
  
==References==

Latest revision as of 18:26, 11 July 2013

Complex map of function $\sqrt(\exp)= \exp^{1/2}$, Halfiteration of exponential to base $\mathrm e$.

$f=\exp^{1/2}(x+\rm i y$

is shown with lines $u=\Re(f)=\mathrm const$ and lines $v=\Im(f)=\mathrm const$.

C++ generator of curves

// FIles ado.cin and conto.cin and superexp.cin and superlo.cin should be loaded in the working directory in order to compile the code below. Actually, the last two evaluate tetration tet and arctetration ate; routines fsexp.cin and fslog.cin can be used instead. 

#include <math.h> 
#include <stdio.h>
#include <stdlib.h>
#define DB double
#define DO(x,y) for(x=0;x<y;x++)
using namespace std;
#include <complex>
typedef complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
//#include "fac.cin"
//#include "sinc.cin"
//#include "facp.cin"
//#include "afacc.cin"
//#include "superfac.cin"
#include "superex.cin"
#include "superlo.cin"
DB xL=0.31813150520476413;
DB yL=1.3372357014306895;

#include "conto.cin"

int main(){ int j,k,m,n,n1; DB x,y, p,q, t; z_type z,c,d;
 int M=160,M1=M+1;
 int N=165,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("fig4b.eps","w");ado(o,402,402);
FILE *o;o=fopen("qexpMap.eps","w");ado(o,402,402);
fprintf(o,"201 201 translate\n 20 20 scale\n");
DO(m,M1) X[m]=-8.+.1*m;
DO(n,N1){   y=-8.+.1*n; if(y< -yL) Y[n]=y; else break;}
        Y[n]=-yL-.001; n++;
        Y[n]=-yL+.001; n++;
for(j=n;j<N1;j++){y=-8.+.1*(j-2); if(y<-.02) Y[j]=y; else break;}
        Y[j]= -.02; j++;
        Y[j]= +.02; j++;
for(k=j;k<N1;k++){y=-8.+.1*(k-3); if(y<yL) Y[k]=y; else break;}
        Y[k]= yL-.001; k++;
        Y[k]= yL+.001; k++;
for(n=k;n<N1;n++){y=-8.+.1*(n-5); Y[n]=y;}
for(m=-8;m<9;m++){if(m==0){M(m,-8.5)L(m,8.5)} else{M(m,-8)L(m,8)}}
for(n=-8;n<9;n++){      M(  -8,n)L(8,n)}
fprintf(o,".008 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]; //printf("%5.2f\n",x);
DO(n,N1){y=Y[n]; z=z_type(x,y); 
//      c=afacc(z);
//      c=fac(c);
//      c=arcsuperfac(z);
//      c=superfac(.5+c);
        c=FSLOG(z);
        c=FSEXP(.5+c);
//      d=z;
//      p=abs(c-d)/(abs(c)+abs(d));  p=-log(p)/log(10.)-1.;
        p=Re(c);q=Im(c);        
        if(p>-999 && p<999)     g[m*N1+n]=p;
        if(q>-999 && q<999 && fabs(q)> 1.e-14) f[m*N1+n]=q;
                        }}
fprintf(o,"1 setlinejoin 1 setlinecap\n"); p=6.;q=6.;
//#include"plofu.cin"
//p=2;q=1;
for(m=-4;m<4;m++) for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M-20,N,(m+.1*n),-q, q); fprintf(o,".02 W 0 .6 0 RGB S\n");
for(m= 0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M-20,N,-(m+.1*n),-q, q); fprintf(o,".02 W .9 0 0 RGB S\n");
for(m= 0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M-20,N, (m+.1*n),-q, q); fprintf(o,".02 W 0 0 .9 RGB S\n");
for(m=1;m<9;m++)  conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p); fprintf(o,".05 W .9 0 0 RGB S\n");
for(m=1;m<9;m++)  conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".05 W 0 0 .9 RGB S\n");
                  conto(o,f,w,v,X,Y,M,N, (0.  ),-p,p); fprintf(o,".05 W .6 0 .6 RGB S\n");
for(m=-8;m<9;m++) conto(o,g,w,v,X,Y,M,N, (0.+m),-p,p); fprintf(o,".05 W 0 0 0 RGB S\n");
fprintf(o,"0 setlinejoin 0 setlinecap\n");
M(xL, yL)L(-8.1, yL) 
M(xL,-yL)L(-8.1,-yL) fprintf(o,"                  .06 W 1 1 1 RGB S\n");
DO(m,21){M(xL-.4*m, yL)L(xL-.4*(m+.4), yL)}
DO(m,21){M(xL-.4*m,-yL)L(xL-.4*(m+.4),-yL)} fprintf(o,".1 W 0 0 0 RGB S\n");
//M(xL, yL)L(-8, yL)
//M(xL,-yL)L(-8,-yL) fprintf(o,"[.2 .2]0 setdash .12 W 0 0 0 RGB S\n");
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
//      system(    "open fig4b.eps");           //for macintosh
       system("epstopdf qexpMap.eps");
       system(    "open qexpMap.pdf");         //for  unix
//      system(    "xpdf fig4b.pdf &");         //for  unix
//      getchar(); system("killall Preview");   //for macintosh
}

Latex generator of labels

\documentclass[12pt]{article}
\usepackage{graphicx}
\usepackage{geometry}
\paperwidth 340pt
\paperheight 340pt
\usepackage{rotating}
\textwidth 800mm
\textheight 400mm
\topmargin -61mm
\oddsidemargin -77pt
\parindent 0pt
\newcommand \sx {\scalebox}
\newcommand \rmi Template:\rm i %imaginary unity is always roman font
\newcommand \rot {\begin{rotate}}'"`UNIQ--item-1--QINU`"' \newcommand \ero {\end{rotate}}
\newcommand \ax {
\put( 10,346){\sx{1.4}{$y$}}
\put( 10,307){\sx{1.3}{$6$}}
\put( 10,267){\sx{1.3}{$4$}}
\put( 10,227){\sx{1.3}{$2$}}
\put( 10,187){\sx{1.3}{$0$}}
\put( 0,147){\sx{1.3}{$-2$}}
\put( 0,107){\sx{1.3}{$-4$}}
\put( 0, 67){\sx{1.3}{$-6$}}
\put( 0, 27){\sx{1.3}{$-8$}}
\put( 50, 17){\sx{1.3}{$-6$}}
\put( 90, 17){\sx{1.3}{$-4$}}
\put(130, 17){\sx{1.3}{$-2$}}
\put(178, 17){\sx{1.3}{$0$}}
\put(218, 17){\sx{1.3}{$2$}}
\put(258, 17){\sx{1.3}{$4$}}
\put(298, 17){\sx{1.3}{$6$}}
\put(334, 17){\sx{1.4}{$x$}}
}
\begin {document}
{\begin{picture}(420,420)'"`UNIQ--item-1--QINU`"' %\put(-20,-10){\includegraphics{fig4b}}\ax'"`UNIQ--item-1--QINU`"' \put(-20,-10){\includegraphics{qexpMap}}\ax'"`UNIQ--item-1--QINU`"' \put( 52,188){\sx{1.1}{$v\!=\!0$}}'"`UNIQ--item-1--QINU`"' \put( 81,176){\sx{1.1}{$u\!=\!-1.6$}}'"`UNIQ--item-1--QINU`"' {\put(40,213){\sx{1.8}{\bf cut}}}'"`UNIQ--item-1--QINU`"' {\put(40,160){\sx{1.8}{\bf cut}}}'"`UNIQ--item-1--QINU`"' \put(300,318){\rot{19}\sx{1.3}{$u\!=\!-8$}\ero}'"`UNIQ--item-1--QINU`"' \put(310,266){\rot{28}\sx{1.3}{$u\!= 8$}\ero}'"`UNIQ--item-1--QINU`"' \put(300,218){\rot{-8}\sx{1.3}{$v\!= 8$}\ero}'"`UNIQ--item-1--QINU`"' \put(300,156){\rot{ 8}\sx{1.3}{$v\!=\!-8$}\ero}'"`UNIQ--item-1--QINU`"' \put(308,108){\rot{-30}\sx{1.3}{$u\!= 8$}\ero}'"`UNIQ--item-1--QINU`"' \put(300, 56){\rot{-22}\sx{1.3}{$u\!=\!-8$}\ero}'"`UNIQ--item-1--QINU`"' \end{picture}}
\end{document}
$

References

Categrory:Tetration

File history

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

Date/TimeThumbnailDimensionsUserComment
current18:03, 11 July 2013Thumbnail for version as of 18:03, 11 July 20131,881 × 1,881 (1.83 MB)T (talk | contribs)Complex map of function $\sqrt(\exp)= \exp^{1/2}$, Halfiteration of exponential to base $\mathrm e$. $f=\exp^{1/2}(x+\rm i y$ is shown with lines $u=\Re(f)=\mathrm const$ and lines $v=\Im(f)=\mathrm const$. ==References== [[Categrory:T...

You cannot overwrite this file.

The following page uses 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:QexpMapT400.jpg&oldid=6574"