Difference between revisions of "File:IterEq2plotU.png"

From TORI
Jump to: navigation, search
(Importing image file)
 
Line 1: Line 1:
  +
[[Explicit plot]] of $c$th [[iteration]] of [[exponential]] to [[base sqrt(2)]] for various values of the number $c$ of iterations.
−
Importing image file
  
  +
  +
For evaluation of the non-integer iteration, the plotter uses the implementation through the [[superfunction]] $F$ of the exponential to base $\sqrt{2}$, constructed at the fixed point $L\!=\!4$, and the corresponding [[Abel function]] $G$:
  +
  +
:$ \exp_b^{c}(x)=F\big(c+G(x)\big)$
  +
  +
In the publication
  +
<ref>
  +
http://tori.ils.uec.ac.jp/PAPERS/2010sqrt2.pdf D.Kouznetsov, H.Trappmnn. Portrait of the four regular super-exponentials to base sqrt(2). [[Mathematics of Computation]], v.271, July 2010, p.1727-1756.
  +
</ref>,
  +
these F and G are referred as $F_{4,5}$ and $F_{4,5}^{~-1}$, respectively.
  +
  +
==References==
  +
<references/>
  +
  +
==[[C++]] generator of curves==
  +
  +
// Files [[F45E.cin]], [[F45L.cin]] and [[ado.cin]] should be loaded in the working directory in order to compile the [[C++]] code below:
  +
  +
#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 "ado.cin"
  +
  +
#include "F45E.cin"
  +
#include "F45L.cin"
  +
  +
DB B=sqrt(2.);
  +
DB F(DB z) { return exp( exp( log(B)*z));}
  +
DB G(DB z) { return log( log(z) )/log(B);}
  +
  +
main(){ int m,n; double x,y,t; FILE *o;
  +
o=fopen("IterEq2plot.eps","w"); ado(o,1420,1420);
  +
fprintf(o,"1 1 translate 100 100 scale\n");
  +
#define M(x,y) fprintf(o,"%6.3f %6.3f M\n",0.+x,0.+y);
  +
#define L(x,y) fprintf(o,"%6.3f %6.3f L\n",0.+x,0.+y);
  +
for(m=0;m<15;m++) {M(m,0)L(m,14)}
  +
for(m=0;m<15;m++) {M(0,m)L(14,m)}
  +
fprintf(o,"2 setlinecap .01 W S\n");
  +
DO(m,82){x=0.001+.1*m; y=exp(log(B)*x); y=exp(log(B)*y); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 0 1 0 RGB S\n");
  +
DO(m,82){x=0.001+.1*m; y=exp(log(B)*x); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 0 1 0 RGB S\n");
  +
DO(m,141){x=0.001+.1*m; y=log(x)/log(B); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 1 0 1 RGB S\n");
  +
DO(m,131){x=1.421+.1*m;y=log(x)/log(B);y=log(y)/log(B); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 1 0 1 RGB S\n");
  +
for(n=-10;n<11;n++){t=.1*n; M(2,2); DO(m,122){x=2.05+.1*m; y=Re(F45E(t+F45L(x+1.e-14*I))); L(x,y); if(y>14.1)break;} }
  +
fprintf(o,"1 setlinecap 1 setlinejoin .02 W 0 0 0 RGB S\n");
  +
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
  +
system("epstopdf IterEq2plot.eps");
  +
system( "open IterEq2plot.pdf");
  +
getchar(); system("killall Preview");
  +
}
  +
  +
==[[Latex]] generator of labels==
  +
  +
%<nowiki> %<br>
  +
% file IterEq2plot.pdf should be generated with the code above in order to compile the Latex document below. %<br>
  +
% Copyleft 2012 by Dmitrii Kouznetsov <br> %
  +
\documentclass[12pt]{article} % <br>
  +
\usepackage{geometry} % <br>
  +
\usepackage{graphicx} % <br>
  +
\usepackage{rotating} % <br>
  +
\paperwidth 1418pt % <br>
  +
\paperheight 1418pt % <br>
  +
\topmargin -94pt % <br>
  +
\oddsidemargin -81pt % <br>
  +
\textwidth 1600pt % <br>
  +
\textheight 1600pt % <br>
  +
\pagestyle {empty} % <br>
  +
\newcommand \sx {\scalebox} % <br>
  +
\newcommand \rot {\begin{rotate}} % <br>
  +
\newcommand \ero {\end{rotate}} % <br>
  +
\newcommand \ing {\includegraphics} % <br>
  +
\parindent 0pt% <br>
  +
\pagestyle{empty} % <br>
  +
\begin{document} % <br>
  +
\begin{picture}(1412,1412) % <br>
  +
%\put(10,10){\ing{IterPowPlot}} % <br>
  +
\put(10,10){\ing{IterEq2plot}} % <br>
  +
\put(11,1374){\sx{4.7}{$y\!=\!\exp_{b}^{~ c}(x)$}} % <br>
  +
\put(11,1298){\sx{4}{$13$}} % <br>
  +
\put(11,1198){\sx{4}{$12$}} % <br>
  +
\put(11,1098){\sx{4}{$11$}} % <br>
  +
\put(11,998){\sx{4}{$10$}} % <br>
  +
\put(11,898){\sx{4}{$9$}} % <br>
  +
\put(11,798){\sx{4}{$8$}} % <br>
  +
\put(11,698){\sx{4}{$7$}} % <br>
  +
\put(11,598){\sx{4}{$6$}} % <br>
  +
\put(11,498){\sx{4}{$5$}} % <br>
  +
\put(11,398){\sx{4}{$4$}} % <br>
  +
\put(11,298){\sx{4}{$3$}} % <br>
  +
\put(11,198){\sx{4}{$2$}} % <br>
  +
\put(11,098){\sx{4}{$1$}} % <br>
  +
% <br>
  +
\put(100,16){\sx{4}{$1$}} % <br>
  +
\put(200,16){\sx{4}{$2$}} % <br>
  +
\put(301,16){\sx{4}{$3$}} % <br>
  +
\put(401,16){\sx{4}{$4$}} % <br>
  +
\put(502,16){\sx{4}{$5$}} % <br>
  +
\put(602,16){\sx{4}{$6$}} % <br>
  +
\put(703,16){\sx{4}{$7$}} % <br>
  +
\put(803,16){\sx{4}{$8$}} % <br>
  +
\put(903,16){\sx{4}{$9$}} % <br>
  +
\put(993,16){\sx{4}{$10$}} % <br>
  +
\put(1093,16){\sx{4}{$11$}} % <br>
  +
\put(1193,16){\sx{4}{$12$}} % <br>
  +
\put(1293,16){\sx{4}{$13$}} % <br>
  +
\put(1390,16){\sx{4}{$x$}} % <br>
  +
% <br>
  +
\put(595,1250){\sx{5}{\rot{84}$c\!=\!2$\ero}} % <br>
  +
\put(754,1250){\sx{5}{\rot{78}$c\!=\!1$\ero}} % <br>
  +
\put(809,1250){\sx{5}{\rot{74}$c\!=\!0.8$\ero}} % <br>
  +
\put(882,1250){\sx{5}{\rot{72}$c\!=\!0.6$\ero}} % <br>
  +
\put(971,1250){\sx{5}{\rot{64}$c\!=\!0.4$\ero}} % <br>
  +
\put(1094,1250){\sx{5}{\rot{54}$c\!=\!0.2$\ero}} % <br>
  +
% <br>
  +
\put(1282,1262){\sx{5}{\rot{44}$c\!=\!0$\ero}} % <br>
  +
% <br>
  +
\put(1256,1060){\sx{5}{\rot{34}$c\!=\!-0.2$\ero}} % <br>
  +
\put(1235, 926){\sx{5}{\rot{26}$c\!=\!-0.4$\ero}} % <br>
  +
\put(1230, 840){\sx{5}{\rot{19}$c\!=\!-0.6$\ero}} % <br>
  +
\put(1220, 770){\sx{5}{\rot{15}$c\!=\!-0.8$\ero}} % <br>
  +
\put(1230, 702){\sx{5}{\rot{11}$c\!=\!-1$\ero}} % <br>
  +
\put(1234, 542){\sx{5}{\rot{5}$c\!=\!-2$\ero}} % <br>
  +
%<br>
  +
\put(560, 1032){\sx{5.4}{\rot{83}$y\!=\!b^{b^x}$\ero}} % <br>
  +
\put(674, 1062){\sx{5}{\rot{72}$y\!=\!b^x$\ero}} % <br>
  +
\put(890,610){\sx{5}{\rot{12}$y\!=\!\log_b(x)$\ero}} % <br>
  +
\put(825,484){\sx{5}{\rot{6}$y\!=\!\log_b^{~2}(x)$\ero}} % <br>
  +
\put(600,200){\sx{11}{$b\!=\!\sqrt{2}$}}
  +
\end{picture} % <br>
  +
\end{document} % <br>
  +
%</nowiki>
  +
  +
[[Category:Iteration]]
  +
[[Category:SuperFunctions]]
  +
[[Category:AbelFunctions]]
  +
[[Category:Explicit plot]]
  +
[[Category:C++]]
  +
[[Category:Latex]]

Revision as of 09:43, 21 June 2013

Explicit plot of $c$th iteration of exponential to base sqrt(2) for various values of the number $c$ of iterations.

For evaluation of the non-integer iteration, the plotter uses the implementation through the superfunction $F$ of the exponential to base $\sqrt{2}$, constructed at the fixed point $L\!=\!4$, and the corresponding Abel function $G$:

$ \exp_b^{c}(x)=F\big(c+G(x)\big)$

In the publication [1], these F and G are referred as $F_{4,5}$ and $F_{4,5}^{~-1}$, respectively.

References

  1. ↑ http://tori.ils.uec.ac.jp/PAPERS/2010sqrt2.pdf D.Kouznetsov, H.Trappmnn. Portrait of the four regular super-exponentials to base sqrt(2). Mathematics of Computation, v.271, July 2010, p.1727-1756.

C++ generator of curves

// Files F45E.cin, F45L.cin and ado.cin should be loaded in the working directory in order to compile the C++ code below: 

#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 "ado.cin"

#include "F45E.cin"
#include "F45L.cin"

DB B=sqrt(2.);
DB F(DB z) { return exp( exp( log(B)*z));}
DB G(DB z) { return log( log(z) )/log(B);}

main(){ int m,n; double x,y,t; FILE *o;
o=fopen("IterEq2plot.eps","w"); ado(o,1420,1420);
fprintf(o,"1 1 translate 100 100 scale\n");
#define M(x,y) fprintf(o,"%6.3f %6.3f M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%6.3f %6.3f L\n",0.+x,0.+y);
for(m=0;m<15;m++) {M(m,0)L(m,14)}
for(m=0;m<15;m++) {M(0,m)L(14,m)}
fprintf(o,"2 setlinecap .01 W S\n");
DO(m,82){x=0.001+.1*m; y=exp(log(B)*x); y=exp(log(B)*y); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1  setlinejoin .04 W 0 1 0 RGB S\n");
DO(m,82){x=0.001+.1*m; y=exp(log(B)*x); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1  setlinejoin .04 W 0 1 0 RGB S\n");
DO(m,141){x=0.001+.1*m; y=log(x)/log(B); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 1 0 1 RGB S\n");
DO(m,131){x=1.421+.1*m;y=log(x)/log(B);y=log(y)/log(B); if(m==0)M(x,y) else L(x,y);} fprintf(o,"1 setlinecap 1 setlinejoin .04 W 1 0 1 RGB S\n");
for(n=-10;n<11;n++){t=.1*n; M(2,2); DO(m,122){x=2.05+.1*m; y=Re(F45E(t+F45L(x+1.e-14*I))); L(x,y); if(y>14.1)break;} }
fprintf(o,"1 setlinecap 1 setlinejoin .02 W 0 0 0 RGB S\n");
fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
      system("epstopdf IterEq2plot.eps");    
      system(    "open IterEq2plot.pdf");
      getchar(); system("killall Preview");
}

Latex generator of labels

% %<br> % file IterEq2plot.pdf should be generated with the code above in order to compile the Latex document below. %<br> % Copyleft 2012 by Dmitrii Kouznetsov <br> % \documentclass[12pt]{article} % <br> \usepackage{geometry} % <br> \usepackage{graphicx} % <br> \usepackage{rotating} % <br> \paperwidth 1418pt % <br> \paperheight 1418pt % <br> \topmargin -94pt % <br> \oddsidemargin -81pt % <br> \textwidth 1600pt % <br> \textheight 1600pt % <br> \pagestyle {empty} % <br> \newcommand \sx {\scalebox} % <br> \newcommand \rot {\begin{rotate}} % <br> \newcommand \ero {\end{rotate}} % <br> \newcommand \ing {\includegraphics} % <br> \parindent 0pt% <br> \pagestyle{empty} % <br> \begin{document} % <br> \begin{picture}(1412,1412) % <br> %\put(10,10){\ing{IterPowPlot}} % <br> \put(10,10){\ing{IterEq2plot}} % <br> \put(11,1374){\sx{4.7}{$y\!=\!\exp_{b}^{~ c}(x)$}} % <br> \put(11,1298){\sx{4}{$13$}} % <br> \put(11,1198){\sx{4}{$12$}} % <br> \put(11,1098){\sx{4}{$11$}} % <br> \put(11,998){\sx{4}{$10$}} % <br> \put(11,898){\sx{4}{$9$}} % <br> \put(11,798){\sx{4}{$8$}} % <br> \put(11,698){\sx{4}{$7$}} % <br> \put(11,598){\sx{4}{$6$}} % <br> \put(11,498){\sx{4}{$5$}} % <br> \put(11,398){\sx{4}{$4$}} % <br> \put(11,298){\sx{4}{$3$}} % <br> \put(11,198){\sx{4}{$2$}} % <br> \put(11,098){\sx{4}{$1$}} % <br> % <br> \put(100,16){\sx{4}{$1$}} % <br> \put(200,16){\sx{4}{$2$}} % <br> \put(301,16){\sx{4}{$3$}} % <br> \put(401,16){\sx{4}{$4$}} % <br> \put(502,16){\sx{4}{$5$}} % <br> \put(602,16){\sx{4}{$6$}} % <br> \put(703,16){\sx{4}{$7$}} % <br> \put(803,16){\sx{4}{$8$}} % <br> \put(903,16){\sx{4}{$9$}} % <br> \put(993,16){\sx{4}{$10$}} % <br> \put(1093,16){\sx{4}{$11$}} % <br> \put(1193,16){\sx{4}{$12$}} % <br> \put(1293,16){\sx{4}{$13$}} % <br> \put(1390,16){\sx{4}{$x$}} % <br> % <br> \put(595,1250){\sx{5}{\rot{84}$c\!=\!2$\ero}} % <br> \put(754,1250){\sx{5}{\rot{78}$c\!=\!1$\ero}} % <br> \put(809,1250){\sx{5}{\rot{74}$c\!=\!0.8$\ero}} % <br> \put(882,1250){\sx{5}{\rot{72}$c\!=\!0.6$\ero}} % <br> \put(971,1250){\sx{5}{\rot{64}$c\!=\!0.4$\ero}} % <br> \put(1094,1250){\sx{5}{\rot{54}$c\!=\!0.2$\ero}} % <br> % <br> \put(1282,1262){\sx{5}{\rot{44}$c\!=\!0$\ero}} % <br> % <br> \put(1256,1060){\sx{5}{\rot{34}$c\!=\!-0.2$\ero}} % <br> \put(1235, 926){\sx{5}{\rot{26}$c\!=\!-0.4$\ero}} % <br> \put(1230, 840){\sx{5}{\rot{19}$c\!=\!-0.6$\ero}} % <br> \put(1220, 770){\sx{5}{\rot{15}$c\!=\!-0.8$\ero}} % <br> \put(1230, 702){\sx{5}{\rot{11}$c\!=\!-1$\ero}} % <br> \put(1234, 542){\sx{5}{\rot{5}$c\!=\!-2$\ero}} % <br> %<br> \put(560, 1032){\sx{5.4}{\rot{83}$y\!=\!b^{b^x}$\ero}} % <br> \put(674, 1062){\sx{5}{\rot{72}$y\!=\!b^x$\ero}} % <br> \put(890,610){\sx{5}{\rot{12}$y\!=\!\log_b(x)$\ero}} % <br> \put(825,484){\sx{5}{\rot{6}$y\!=\!\log_b^{~2}(x)$\ero}} % <br> \put(600,200){\sx{11}{$b\!=\!\sqrt{2}$}} \end{picture} % <br> \end{document} % <br> %

File history

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

Date/TimeThumbnailDimensionsUserComment
current17:50, 20 June 2013Thumbnail for version as of 17:50, 20 June 20132,944 × 2,944 (986 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:IterEq2plotU.png&oldid=5233"