Difference between revisions of "File:Straro20at.jpg"

From TORI
Jump to navigation Jump to search
(Importing image file)
 
 
Line 1: Line 1:
  +
[[Complex map]] of function [[StraRo]]:
−
Importing image file
  
  +
  +
$u\!+\!\mathrm i v=\mathrm{Straro}_q(x\!+\!\mathrm i y)$
  +
  +
For positive $q$ and complez $z$, function $\mathrm{Straro}_q(z)$ is such that
  +
  +
$\mathrm{StraRo}_q(z)^2=27(z\!-\!q)^2+4(1\!+\!4qz)^3$
  +
  +
where $q\!>\!0$
  +
and the cut lines are chosen in the special way:
  +
  +
From the negative branch point, there is cut line along the negative part of the real axis to $-\infty$.
  +
  +
Also, from each complex branch point, there is straight line to zero.
  +
  +
In the rest of the complex plane, for positive $q$, function $\mathrm{StraRo}$ is holomorphic.
  +
  +
The figure shows the map for $q\!=\!2$.
  +
  +
The special choice of the cuts above is used to define function [[ArcNem]] (with the same structure of cut lines), and it is used for evaluation of function [[AuNem]], which is [[Abel function]] for the [[Nemtsov function]]. This function is proposed al illustration of the exotic iterates of the [[Nemtsov function]] constructed as zero.
  +
  +
The upper complex branch point is expressed with function [[NemBran]],
  +
  +
$x_0\!+\!\mathrm i y_0=\mathrm{NemBran}(q)$
  +
  +
The complex brach points are shown with additional grid lines
  +
$x=x_0$ and
  +
$y=\pm y_0$
  +
  +
  +
For $q=2$, the evaluation suggests
  +
  +
$x_0=0.089910637427530 ~$, $~ y_0= 0.296978217844034$
  +
  +
<!--
  +
$0.07226562500000~$, $~ y_0=0.340398926912761$
  +
!-->
  +
  +
==References==
  +
<references/>
  +
  +
[[Category:Book]]
  +
[[Category:BookMap]]
  +
[[Category:C++]]
  +
[[Category:Branch point]]
  +
[[Category:Grid line]]
  +
[[Category:Latex]]
  +
[[Category:Nemtsov function]]
  +
[[Category:Superfunctions]]
  +
  +
==[[C++]] generator of map==
  +
Files
  +
[[ado.cin]],
  +
[[conto.cin]],
  +
[[nembran.cin]]
  +
should be loaded in order to compile the [[C++]] 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++)
  +
//using namespace std;
  +
#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"
  +
  +
DB Q=2;
  +
z_type T(z_type z){ return z*(1.+z*z*(1.+z*Q));}
  +
  +
#include "nembran.cin"
  +
  +
z_type S=nembran(Q);
  +
  +
z_type F(z_type z){ z_type a,b,r,R ; a=Q-z; b=1. + 4.*Q*z; r = 27.*a*a + 4.*b*b*b ;
  +
DB x,y; x=Re(z); y=Im(z);
  +
if( x<=0 ) return -sqrt(r);
  +
if( y > Im(S) ) return I*sqrt(-r);
  +
if( y < - Im(S) ) return -I*sqrt(-r);
  +
if( x*Im(S)< fabs(y)*Re(S) ) return -sqrt(r);
  +
return sqrt(r);
  +
}
  +
  +
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
  +
int M=1001,M1=M+1;
  +
int N=501,N1=N+1;
  +
DB X[M1],Y[N1]; DB *g, *f, *w; // w is working array.
  +
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
char v[M1*N1]; // v is working array
  +
FILE *o;o=fopen("straro20a.eps","w"); ado(o,4002,4002);
  +
fprintf(o,"2001 2001 translate\n 1000 1000 scale\n");
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
DO(m,M1) X[m]=-2+.004*(m-.5);
  +
DO(n,N1) Y[n]=-2+.008*(n-.5);
  +
//for(n=0;n<N1;n++) { Y[n]=1.09*sinh((3./200.)*(n-200)); printf("%3d %9.6f\n",n,Y[n]); }
  +
for(m=-20;m<21;m+=5){M(.1*m,-2) L(.1*m,2) }
  +
for(n=-20;n<21;n+=5){M( -2,.1*n) L(2,.1*n)}
  +
fprintf(o,".001 W 0 0 0 RGB S\n");
  +
  +
//z=nembran(Q);
  +
x=Re(S);
  +
y=Im(S);
  +
M(x,-1)L(x,1)
  +
M(-1,y)L(1,y)
  +
M(-1,-y)L(1,-y)
  +
fprintf(o,".004 W .5 .5 0 RGB S\n");
  +
printf("x,y %18.15lf %18.15lf\n",x,y);
  +
  +
  +
DO(m,M1)DO(n,N1){ g[m*N1+n]=99999999;
  +
f[m*N1+n]=99999999;}
  +
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
  +
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)
  +
  +
// c=-1./z;
  +
// c=-(1.+3.*z*z)/(4.*z*z*z);
  +
  +
//c=Fru(z);
  +
c=F(z);
  +
  +
// p=abs(c-z)/(abs(c)+abs(z)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
  +
p=Re(c); q=Im(c); if(p>-100000 && p<100000 && q>-100000 && q<100000 ){ g[m*N1+n]=p;f[m*N1+n]=q;}
  +
}}
  +
fprintf(o,"1 setlinejoin 1 setlinecap\n");
  +
p=4;q=.5;
  +
  +
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,".001 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,".001 W .8 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,".001 W 0 0 .8 RGB S\n");
  +
for(m= 1;m<41;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".004 W .8 0 0 RGB S\n");
  +
for(m= 1;m<41;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".004 W 0 0 .8 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".003 W .6 0 .6 RGB S\n");
  +
for(m=-40;m<41;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".004 W 0 0 0 RGB S\n");
  +
// fprintf(o,"0 setlinejoin 0 setlinecap\n");
  +
// M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .02 W S\n");
  +
//#include "plofu.cin"
  +
  +
x=Re(S);
  +
y=Im(S);
  +
M(x,y)L(0,0)
  +
L(x,-y)
  +
fprintf(o,".01 W 1 1 .5 RGB 0 setlinecap S\n");
  +
printf("x,y %18.15lf %18.15lf\n",x,y);
  +
  +
fprintf(o,"showpage\n");
  +
fprintf(o,"%c%cTrailer\n",'%','%');
  +
fclose(o); free(f); free(g); free(w);
  +
  +
system("epstopdf straro20a.eps");
  +
system( "open straro20a.pdf"); //for macintosh
  +
getchar(); system("killall Preview"); // For macintosh
  +
}
  +
//
  +
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] generator of labels==
  +
  +
<poem><nomathjax><nowiki>
  +
\documentclass[12pt]{article}
  +
\usepackage{geometry}
  +
\paperwidth 4200pt
  +
\paperheight 4200pt
  +
\topmargin -100pt
  +
\oddsidemargin -80pt
  +
\textheight 4800px
  +
\parindent 0pt
  +
\usepackage{graphicx}
  +
\usepackage{rotating}
  +
\newcommand \rot {\begin{rotate}}
  +
\newcommand \ero {\end{rotate}}
  +
\newcommand \ing {\includegraphics}
  +
\newcommand \sx {\scalebox}
  +
\begin{document}
  +
\begin{picture}(4090,4050)
  +
%\put(20,14){\ing{nembran}}
  +
%\put(190,30){\ing{straro10a}}
  +
\put(190,30){\ing{straro20a}}
  +
\put(70,3940){\sx{14}{$y$}}
  +
\put(70,2970){\sx{14}{$1$}}
  +
\put(70,1970){\sx{14}{$0$}}
  +
\put(-30, 970){\sx{14}{$-1$}}
  +
\put(-30,-20){\sx{14}{$-2$}}
  +
\put( 50,-120){\sx{14}{$-2$}}
  +
\put(1050,-120){\sx{14}{$-1$}}
  +
\put(2170,-120){\sx{14}{$0$}}
  +
\put(3170,-120){\sx{14}{$1$}}
  +
\put(4110,-120){\sx{14}{$x$}}
  +
\put(330,2000){\sx{14}{\bf cut}}
  +
\put(1750,1994){\sx{14}{$v\!=\!0$}}
  +
\put(3330,2000){\sx{14}{$v\!=\!0$}}
  +
\put(2230,2060){\rot{69}\sx{14}{\bf cut}\ero}
  +
%
  +
\put(2180,3300){\rot{49}\sx{14}{$u\!=\!-40$}\ero}%
  +
\put(2450,3020){\rot{54}\sx{14}{$u\!=\!-10$}\ero}%
  +
\put(2700,3000){\rot{59}\sx{14}{$u\!=\!0$}\ero}%
  +
\put(2760,2800){\rot{68}\sx{14}{$u\!=\!10$}\ero}%
  +
\put(3180,2700){\rot{68}\sx{14}{$u\!=\!40$}\ero}%
  +
%
  +
\put(800, 3160){\rot{-72}\sx{14}{$v\!=\!-40$}\ero}%
  +
\put(1140, 3300){\rot{-64}\sx{14}{$v\!=\!-10$}\ero}%
  +
\put(1260, 3360){\rot{-64}\sx{14}{$v\!=\!0$}\ero}%
  +
\put(1360, 3420){\rot{-60}\sx{14}{$v\!=\!10$}\ero}%
  +
%
  +
\put(400, 1870){\rot{-5}\sx{14}{$u\!=\!-10$}\ero}
  +
\put(400, 1520){\rot{-7}\sx{14}{$u\!=\!-40$}\ero}
  +
%
  +
\put(760, 720){\rot{66}\sx{14}{$v\!=\!40$}\ero}
  +
\put(1090, 600){\rot{65}\sx{14}{$v\!=\!10$}\ero}
  +
\put(1200, 500){\rot{58}\sx{14}{$v\!=\!0$}\ero}
  +
\put(1320, 370){\rot{57}\sx{14}{$v\!=\!-10$}\ero}
  +
%
  +
\put(3300,2560){\rot{-14}\sx{14}{$v\!=\!40$}\ero}
  +
\put(3280,2160){\rot{-2}\sx{14}{$v\!=\!10$}\ero}
  +
\put(3190,1840){\rot{2}\sx{14}{$v\!=\!-10$}\ero}
  +
\put(3300, 1440){\rot{10}\sx{14}{$v\!=\!-40$}\ero}
  +
\end{picture}
  +
\end{document}
  +
  +
</nowiki></nomathjax></poem>

Latest revision as of 08:52, 1 December 2018

Complex map of function StraRo:

$u\!+\!\mathrm i v=\mathrm{Straro}_q(x\!+\!\mathrm i y)$

For positive $q$ and complez $z$, function $\mathrm{Straro}_q(z)$ is such that

$\mathrm{StraRo}_q(z)^2=27(z\!-\!q)^2+4(1\!+\!4qz)^3$

where $q\!>\!0$ and the cut lines are chosen in the special way:

From the negative branch point, there is cut line along the negative part of the real axis to $-\infty$.

Also, from each complex branch point, there is straight line to zero.

In the rest of the complex plane, for positive $q$, function $\mathrm{StraRo}$ is holomorphic.

The figure shows the map for $q\!=\!2$.

The special choice of the cuts above is used to define function ArcNem (with the same structure of cut lines), and it is used for evaluation of function AuNem, which is Abel function for the Nemtsov function. This function is proposed al illustration of the exotic iterates of the Nemtsov function constructed as zero.

The upper complex branch point is expressed with function NemBran,

$x_0\!+\!\mathrm i y_0=\mathrm{NemBran}(q)$

The complex brach points are shown with additional grid lines $x=x_0$ and $y=\pm y_0$


For $q=2$, the evaluation suggests

$x_0=0.089910637427530 ~$, $~ y_0= 0.296978217844034$


References


C++ generator of map

Files ado.cin, conto.cin, nembran.cin should be loaded 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 std::complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
#include "conto.cin"

DB Q=2;
z_type T(z_type z){ return z*(1.+z*z*(1.+z*Q));}

#include "nembran.cin"

z_type S=nembran(Q);

z_type F(z_type z){ z_type a,b,r,R ; a=Q-z; b=1. + 4.*Q*z; r = 27.*a*a + 4.*b*b*b ;
DB x,y; x=Re(z); y=Im(z);
if( x<=0 ) return -sqrt(r);
if( y > Im(S) ) return I*sqrt(-r);
if( y < - Im(S) ) return -I*sqrt(-r);
if( x*Im(S)< fabs(y)*Re(S) ) return -sqrt(r);
return sqrt(r);
}


int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
int M=1001,M1=M+1;
int N=501,N1=N+1;
DB X[M1],Y[N1]; DB *g, *f, *w; // w is working array.
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
char v[M1*N1]; // v is working array
FILE *o;o=fopen("straro20a.eps","w"); ado(o,4002,4002);
fprintf(o,"2001 2001 translate\n 1000 1000 scale\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
DO(m,M1) X[m]=-2+.004*(m-.5);
DO(n,N1) Y[n]=-2+.008*(n-.5);
//for(n=0;n<N1;n++) { Y[n]=1.09*sinh((3./200.)*(n-200)); printf("%3d %9.6f\n",n,Y[n]); }
for(m=-20;m<21;m+=5){M(.1*m,-2) L(.1*m,2) }
for(n=-20;n<21;n+=5){M( -2,.1*n) L(2,.1*n)}
 fprintf(o,".001 W 0 0 0 RGB S\n");

//z=nembran(Q);
x=Re(S);
y=Im(S);
M(x,-1)L(x,1)
M(-1,y)L(1,y)
M(-1,-y)L(1,-y)
fprintf(o,".004 W .5 .5 0 RGB S\n");
printf("x,y %18.15lf %18.15lf\n",x,y);


DO(m,M1)DO(n,N1){ g[m*N1+n]=99999999;
                  f[m*N1+n]=99999999;}
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)

// c=-1./z;
// c=-(1.+3.*z*z)/(4.*z*z*z);

//c=Fru(z);
c=F(z);

// p=abs(c-z)/(abs(c)+abs(z)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
 p=Re(c); q=Im(c); if(p>-100000 && p<100000 && q>-100000 && q<100000 ){ g[m*N1+n]=p;f[m*N1+n]=q;}
       }}
fprintf(o,"1 setlinejoin 1 setlinecap\n");
 p=4;q=.5;

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,".001 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,".001 W .8 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,".001 W 0 0 .8 RGB S\n");
for(m= 1;m<41;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".004 W .8 0 0 RGB S\n");
for(m= 1;m<41;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".004 W 0 0 .8 RGB S\n");
               conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".003 W .6 0 .6 RGB S\n");
for(m=-40;m<41;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".004 W 0 0 0 RGB S\n");
// fprintf(o,"0 setlinejoin 0 setlinecap\n");
// M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .02 W S\n");
//#include "plofu.cin"

x=Re(S);
y=Im(S);
M(x,y)L(0,0)
L(x,-y)
fprintf(o,".01 W 1 1 .5 RGB 0 setlinecap S\n");
printf("x,y %18.15lf %18.15lf\n",x,y);

fprintf(o,"showpage\n");
fprintf(o,"%c%cTrailer\n",'%','%');
fclose(o); free(f); free(g); free(w);

      system("epstopdf straro20a.eps");
      system( "open straro20a.pdf"); //for macintosh
      getchar(); system("killall Preview"); // For macintosh
}
//

Latex generator of labels


\documentclass[12pt]{article}
\usepackage{geometry}
\paperwidth 4200pt
\paperheight 4200pt
\topmargin -100pt
\oddsidemargin -80pt
\textheight 4800px
\parindent 0pt
\usepackage{graphicx}
\usepackage{rotating}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\newcommand \ing {\includegraphics}
\newcommand \sx {\scalebox}
\begin{document}
\begin{picture}(4090,4050)
%\put(20,14){\ing{nembran}}
%\put(190,30){\ing{straro10a}}
\put(190,30){\ing{straro20a}}
\put(70,3940){\sx{14}{$y$}}
\put(70,2970){\sx{14}{$1$}}
\put(70,1970){\sx{14}{$0$}}
\put(-30, 970){\sx{14}{$-1$}}
\put(-30,-20){\sx{14}{$-2$}}
\put( 50,-120){\sx{14}{$-2$}}
\put(1050,-120){\sx{14}{$-1$}}
\put(2170,-120){\sx{14}{$0$}}
\put(3170,-120){\sx{14}{$1$}}
\put(4110,-120){\sx{14}{$x$}}
\put(330,2000){\sx{14}{\bf cut}}
\put(1750,1994){\sx{14}{$v\!=\!0$}}
\put(3330,2000){\sx{14}{$v\!=\!0$}}
\put(2230,2060){\rot{69}\sx{14}{\bf cut}\ero}
%
\put(2180,3300){\rot{49}\sx{14}{$u\!=\!-40$}\ero}%
\put(2450,3020){\rot{54}\sx{14}{$u\!=\!-10$}\ero}%
\put(2700,3000){\rot{59}\sx{14}{$u\!=\!0$}\ero}%
\put(2760,2800){\rot{68}\sx{14}{$u\!=\!10$}\ero}%
\put(3180,2700){\rot{68}\sx{14}{$u\!=\!40$}\ero}%
%
\put(800, 3160){\rot{-72}\sx{14}{$v\!=\!-40$}\ero}%
\put(1140, 3300){\rot{-64}\sx{14}{$v\!=\!-10$}\ero}%
\put(1260, 3360){\rot{-64}\sx{14}{$v\!=\!0$}\ero}%
\put(1360, 3420){\rot{-60}\sx{14}{$v\!=\!10$}\ero}%
%
\put(400, 1870){\rot{-5}\sx{14}{$u\!=\!-10$}\ero}
\put(400, 1520){\rot{-7}\sx{14}{$u\!=\!-40$}\ero}
%
\put(760, 720){\rot{66}\sx{14}{$v\!=\!40$}\ero}
\put(1090, 600){\rot{65}\sx{14}{$v\!=\!10$}\ero}
\put(1200, 500){\rot{58}\sx{14}{$v\!=\!0$}\ero}
\put(1320, 370){\rot{57}\sx{14}{$v\!=\!-10$}\ero}
%
\put(3300,2560){\rot{-14}\sx{14}{$v\!=\!40$}\ero}
\put(3280,2160){\rot{-2}\sx{14}{$v\!=\!10$}\ero}
\put(3190,1840){\rot{2}\sx{14}{$v\!=\!-10$}\ero}
\put(3300, 1440){\rot{10}\sx{14}{$v\!=\!-40$}\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 20182,905 × 2,905 (1.28 MB)Maintenance script (talk | contribs)Importing image file

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:Straro20at.jpg&oldid=21688"