Difference between revisions of "File:Nembrant.jpg"

From TORI
Jump to: navigation, search
(Importing image file)
 
 
Line 1: Line 1:
  +
[[Parametric plot]] of the branchpoint of function [[ArqNem]]$_q$,
Importing image file
 
  +
The imaginary part of the branch point versus the real part for positive values of parameter $q$:
  +
  +
$x=\Re(\mathrm{NemBran}(q))$
  +
  +
$y=\Im(\mathrm{NemBran}(q))$
  +
  +
where $z=\mathrm{NemBran}(q)$ is branchpoint of function
  +
$z\mapsto\mathrm{ArqNem}_q(z)$
  +
  +
Function [[NemBran]] is expressed through function [[NemBra]],
  +
that returns solution $z=\mathrm{NemBra}(q)$ of equation
  +
  +
$\mathrm{Nem}_q^{\,\prime}(z)=0$
  +
  +
The representation has the following form:
  +
  +
$\mathrm{NemBran}(q)=\mathrm{Nem}_q(\mathrm{NemBra}_q(z))$
  +
  +
While the moguls of the brachpoint is smaller than unity, and the linear term in the representation od the [[Nemtsov Function]] dominates. For this reason, the parametric plot of function [[NemBran]] looks similar to that of function [[NemBra]]. I remind the [[Nemtsov Function]] with parameter $q$ is defiend with
  +
  +
$\mathrm{Nem}_q(z)=z+z^3+q z^4$
  +
  +
Position of the branch point returned with function [[NemBra]] is the same for all the three versions of the inverse of the [[Nemtsov Function]]; not only for [[ArqNem]], but also for [[ArcNem]] and [[ArkNem]].
  +
  +
==Usage==
  +
  +
Function [[NemBran]] is necessary for evaluation of [[ArqNem]]; this function, in its turn, is necessary for the extension of the [[Abel function]] of the [[Nemtsov function]], denoted with [[AuNem]], to the complex plane.
  +
Function [[AuNem]] should be described in the updated version of the book [[Superfunction]] in the example of [[exotic iterate]] of a transfer function $T$ such that
  +
$T(0)\!=\!0$
  +
$T'(0)\!=\!1$
  +
$T''(0)\!=\!0$.
  +
  +
In such a way, the image Nembrant.jpg is prepared for the Book [[Superfunctions]].
  +
  +
One example of such function is $T=\sin$, bit function [[sin]] is symmetric, $\sin(-z)=-\sin(z)$.
  +
[[Superfunction]] and [[Abel Function]] of [[sin]] are denoted with [[SuSin]] and [[AuSin]].....
  +
The corresponding iterates of [[sin]] are reported in 2014,
  +
<ref>
  +
http://www.pphmj.com/references/8246.htm <br>
  +
http://mizugadro.mydns.jp/PAPERS/2014susin.pdf D.Kouznetsov. Super sin. Far East Jourmal of Mathematical Science, v.85, No.2, 2014, pages 219-238.
  +
</ref>
  +
  +
The example with [[Nemtsov function]] is important to show, that the symmetry is not essential for the construction of [[exotic iterate]]s; they can be constructed for asymmetric functions too.
  +
  +
==References==
  +
<references/>
  +
  +
==[[C++]] generator of the plot==
  +
  +
<poem><nomathjax><nowiki>
  +
#include <stdio.h>
  +
#include <math.h>
  +
#include <stdlib.h>
  +
#include <complex>
  +
typedef std::complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
#define DB double
  +
#define DO(x,y) for(x=0;x<y;x++)
  +
  +
void ado(FILE *O, int X, int Y)
  +
{ fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%');
  +
fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y);
  +
fprintf(O,"/M {moveto} bind def\n");
  +
fprintf(O,"/L {lineto} bind def\n");
  +
fprintf(O,"/S {stroke} bind def\n");
  +
fprintf(O,"/s {show newpath} bind def\n");
  +
fprintf(O,"/C {closepath} bind def\n");
  +
fprintf(O,"/F {fill} bind def\n");
  +
fprintf(O,"/o {.002 0 360 arc C S} bind def\n");
  +
fprintf(O,"/times-Roman findfont 20 scalefont setfont\n");
  +
fprintf(O,"/W {setlinewidth} bind def\n");
  +
fprintf(O,"/RGB {setrgbcolor} bind def\n");}
  +
//#include "ado.cin"
  +
  +
z_type nem(DB q,z_type z){ return z*(1.+z*z*(1.+z*q)); }
  +
z_type nem1(DB q,z_type z){ return 1.+z*z*(3.+z*(4.*q)); }
  +
z_type nembra0(DB q){
  +
return -0.5773502691896258*I+
  +
q*(2./9+
  +
q*(0.2138334330331947*I+
  +
q*(-0.2633744855967078 +
  +
q*(-0.3658927631901332*I+
  +
q*(0.5462581923487273 +
  +
q*(0.8556857213229570*I+
  +
q*(-1.387322393266609
  +
))))))) ;}
  +
  +
z_type nembrao(DB q){ z_type x,y,z,s; x=conj(pow(z_type(-.25/q,0.),1./3.)); y=x*x; z=y*y;
  +
s=1.+y*(1.+y*(1.+y*(2./3.+z*(-2./3.+y*(-7./9.+z*(11./9.+y*(130./91.)
  +
) ) ) ) ) );
  +
return s*x;}
  +
  +
z_type nembra(DB q){
  +
if(fabs(q)<.021) return nembra0(q);
  +
if(fabs(q) >20.) return nembrao(q);
  +
z_type Q,v,V; Q=q*q; v=-1.-8.*Q+4.*sqrt(Q+4.*Q*Q); V=pow(v,1./3.); return (.25/q)*(-1.+1./V+V); }
  +
  +
z_type nembran(DB q){ return nem(q,nembra(q));}
  +
  +
int main(){ int m,n; z_type z,c,d,e; DB x,y,q,u,v;
  +
  +
DO(n,41){ q=-.04+.01*n; c=nembra0(q); d=nem(q,c); e=nem1(q,c);
  +
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}
  +
printf("\n");
  +
DO(n,41){ q=.1*n; c=nembra(q); d=nem(q,c); e=nem1(q,c);
  +
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}
  +
printf("\n");
  +
DO(n,41){ q=1.+1.*n; c=nembrao(q); d=nem(q,c); e=nem1(q,c);
  +
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}
  +
  +
//FILE *o; o=fopen("nembraplo.eps","w"); ado(o,154,604);
  +
FILE *o; o=fopen("nembran.eps","w"); ado(o,154,404);
  +
fprintf(o,"2 2 translate 1000 1000 scale 2 setlinecap\n");
  +
#define M(x,y) fprintf(o,"%8.6f %8.6f M\n",(0.+x),(0.+y));
  +
#define L(x,y) fprintf(o,"%8.6f %8.6f L\n",(0.+x),(0.+y));
  +
#define o(x,y) fprintf(o,"%8.6f %8.6f o\n",(0.+x),(0.+y));
  +
  +
for(n=0;n<41;n+=5){ M(0,.01*n) L(.15,.01*n) }
  +
for(n=0;n<31;n+=5){ M(.01*n,0) L(.01*n,.4) }
  +
fprintf(o,".0006 W S 1 setlinejoin\n");
  +
fprintf(o,"1 0 0 RGB .001 W\n");
  +
  +
for(n= 1;n<10;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
//for(n=11;n<20;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
//for(n=21;n<30;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
fprintf(o,"0 .5 0 RGB\n");
  +
for(n= 1;n<10;n++){ q=1.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
fprintf(o,"0 0 1 RGB\n");
  +
for(n= 1;n<10;n++){ q=10.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
fprintf(o,".5 0 .8 RGB\n");
  +
for(n= 1;n<10;n++){ q=100.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
fprintf(o,"0 0 0 RGB\n");
  +
for(n= 0;n<11;n+=10){ q=100.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
  +
q=10000.; c=nembran(q); x=Re(c); y=-Im(c); o(x,y);
  +
  +
for(n= 0;n<400;n++){ q=-.006+.04*(n+.01*n*n+.001*n*n*n); c=nembran(q); x=Re(c); y=-Im(c); if(n==0) M(x,y) else L(x,y);
  +
printf("%9.3lf %16.6lf %16.6lf\n",q,x,y);
  +
}
  +
L(0,0)
  +
fprintf(o,"0 0 0 RGB .001 W 0 setlinecap S\n");
  +
  +
fprintf(o,"showpage\n%c%cTrailer\n",'%','%');
  +
fclose(o);
  +
  +
system("epstopdf nembran.eps");
  +
system("open nembran.pdf");
  +
getchar();
  +
return 0;}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] generator of labels==
  +
  +
<poem><nomathjax><nowiki>
  +
\documentclass[12pt]{article}
  +
\usepackage{geometry}
  +
\paperwidth 174pt
  +
\paperheight 420pt
  +
\topmargin -100pt
  +
\oddsidemargin -72pt
  +
\textheight 800px
  +
\parindent 0pt
  +
\usepackage{graphicx}
  +
\usepackage{rotating}
  +
\newcommand \ing {\includegraphics}
  +
\newcommand \sx {\scalebox}
  +
\begin{document}
  +
\begin{picture}(170,410)
  +
%\put(20,14){\ing{nembraplo}}
  +
\put(20,14){\ing{nembran}}
  +
\put(8,411){\sx{1.4}{$y$}}
  +
%\put(0,512){\sx{1.4}{$0.5$}}
  +
%\put(0,412){\sx{1.4}{$0.4$}}
  +
\put(0,312){\sx{1.4}{$0.3$}}
  +
\put(0,212){\sx{1.4}{$0.2$}}
  +
\put(0,112){\sx{1.4}{$0.1$}}
  +
\put(8,12){\sx{1.4}{$0$}}
  +
\put(20,1){\sx{1.4}{$0$}}
  +
\put(60,1){\sx{1.4}{$0.05$}}
  +
\put(112,1){\sx{1.4}{$0.1$}}
  +
\put(162,2){\sx{1.4}{$x$}}
  +
%\put(18,601){\sx{1.1}{$q\!=\!0$}}
  +
\put(0,389){\sx{1.1}{$q\!=\!0$}}
  +
%\put(63,593){\sx{1.1}{$q\!=\!0.2$}}
  +
\put(43,404){\sx{1.1}{$q\!=\!0.2$}}
  +
%\put(97,574){\sx{1.1}{$q\!=\!0.4$}}
  +
\put(62,394){\sx{1.1}{$q\!=\!0.4$}}
  +
\put(78,384){\sx{1.1}{$q\!=\!0.6$}}
  +
\put(88,372){\sx{1.1}{$q\!=\!0.8$}}
  +
\put(97,356){\sx{1.1}{$q\!=\!1$}}
  +
\put(116,312){\sx{1.1}{$q\!=\!2$}}
  +
\put(120,284){\sx{1.1}{$q\!=\!3$}}
  +
\put(122,263){\sx{1.1}{$q\!=\!4$}}
  +
\put(122,250){\sx{1.1}{$q\!=\!5$}}
  +
\put(112,200){\sx{1.1}{$q\!=\!10$}}
  +
\put(104,164){\sx{1.1}{$q\!=\!20$}}
  +
\put(94,145){\sx{1.1}{$q\!=\!30$}}
  +
\put( 74,99){\sx{1.1}{$q\!=\!100$}}
  +
\put( 49,52){\sx{1.1}{$q\!=\!1000$}}
  +
\put( 36,32){\sx{1.1}{$q\!=\!10000$}}
  +
\end{picture}
  +
\end{document}
  +
</nowiki></nomathjax></poem>
  +
  +
[[Category:Book]]
  +
[[Category:BookPlot]]
  +
[[Category:C++]]
  +
[[Category:Latex]]
  +
[[Category:NemBrant]]
  +
[[Category:ArqNem]]
  +
[[Category:SuNem]]
  +
[[Category:AuNem]]
  +
[[Category:Nemtsov function]]
  +
[[Category:Parametric plot]]

Latest revision as of 08:44, 1 December 2018

Parametric plot of the branchpoint of function ArqNem$_q$, The imaginary part of the branch point versus the real part for positive values of parameter $q$:

$x=\Re(\mathrm{NemBran}(q))$

$y=\Im(\mathrm{NemBran}(q))$

where $z=\mathrm{NemBran}(q)$ is branchpoint of function $z\mapsto\mathrm{ArqNem}_q(z)$

Function NemBran is expressed through function NemBra, that returns solution $z=\mathrm{NemBra}(q)$ of equation

$\mathrm{Nem}_q^{\,\prime}(z)=0$

The representation has the following form:

$\mathrm{NemBran}(q)=\mathrm{Nem}_q(\mathrm{NemBra}_q(z))$

While the moguls of the brachpoint is smaller than unity, and the linear term in the representation od the Nemtsov Function dominates. For this reason, the parametric plot of function NemBran looks similar to that of function NemBra. I remind the Nemtsov Function with parameter $q$ is defiend with

$\mathrm{Nem}_q(z)=z+z^3+q z^4$

Position of the branch point returned with function NemBra is the same for all the three versions of the inverse of the Nemtsov Function; not only for ArqNem, but also for ArcNem and ArkNem.

Usage

Function NemBran is necessary for evaluation of ArqNem; this function, in its turn, is necessary for the extension of the Abel function of the Nemtsov function, denoted with AuNem, to the complex plane. Function AuNem should be described in the updated version of the book Superfunction in the example of exotic iterate of a transfer function $T$ such that $T(0)\!=\!0$ $T'(0)\!=\!1$ $T(0)\!=\!0$.

In such a way, the image Nembrant.jpg is prepared for the Book Superfunctions.

One example of such function is $T=\sin$, bit function sin is symmetric, $\sin(-z)=-\sin(z)$. Superfunction and Abel Function of sin are denoted with SuSin and AuSin..... The corresponding iterates of sin are reported in 2014, [1]

The example with Nemtsov function is important to show, that the symmetry is not essential for the construction of exotic iterates; they can be constructed for asymmetric functions too.

References

  1. http://www.pphmj.com/references/8246.htm
    http://mizugadro.mydns.jp/PAPERS/2014susin.pdf D.Kouznetsov. Super sin. Far East Jourmal of Mathematical Science, v.85, No.2, 2014, pages 219-238.

C++ generator of the plot


#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <complex>
typedef std::complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
#define DB double
#define DO(x,y) for(x=0;x<y;x++)

void ado(FILE *O, int X, int Y)
{ fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%');
        fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y);
        fprintf(O,"/M {moveto} bind def\n");
        fprintf(O,"/L {lineto} bind def\n");
        fprintf(O,"/S {stroke} bind def\n");
        fprintf(O,"/s {show newpath} bind def\n");
        fprintf(O,"/C {closepath} bind def\n");
        fprintf(O,"/F {fill} bind def\n");
        fprintf(O,"/o {.002 0 360 arc C S} bind def\n");
        fprintf(O,"/times-Roman findfont 20 scalefont setfont\n");
        fprintf(O,"/W {setlinewidth} bind def\n");
        fprintf(O,"/RGB {setrgbcolor} bind def\n");}
//#include "ado.cin"

z_type nem(DB q,z_type z){ return z*(1.+z*z*(1.+z*q)); }
z_type nem1(DB q,z_type z){ return 1.+z*z*(3.+z*(4.*q)); }
z_type nembra0(DB q){
return -0.5773502691896258*I+
        q*(2./9+
        q*(0.2138334330331947*I+
        q*(-0.2633744855967078 +
        q*(-0.3658927631901332*I+
        q*(0.5462581923487273 +
        q*(0.8556857213229570*I+
        q*(-1.387322393266609
        ))))))) ;}

z_type nembrao(DB q){ z_type x,y,z,s; x=conj(pow(z_type(-.25/q,0.),1./3.)); y=x*x; z=y*y;
s=1.+y*(1.+y*(1.+y*(2./3.+z*(-2./3.+y*(-7./9.+z*(11./9.+y*(130./91.)
       ) ) ) ) ) );
return s*x;}

z_type nembra(DB q){
if(fabs(q)<.021) return nembra0(q);
if(fabs(q) >20.) return nembrao(q);
z_type Q,v,V; Q=q*q; v=-1.-8.*Q+4.*sqrt(Q+4.*Q*Q); V=pow(v,1./3.); return (.25/q)*(-1.+1./V+V); }

z_type nembran(DB q){ return nem(q,nembra(q));}

int main(){ int m,n; z_type z,c,d,e; DB x,y,q,u,v;

DO(n,41){ q=-.04+.01*n; c=nembra0(q); d=nem(q,c); e=nem1(q,c);
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}
printf("\n");
DO(n,41){ q=.1*n; c=nembra(q); d=nem(q,c); e=nem1(q,c);
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}
printf("\n");
DO(n,41){ q=1.+1.*n; c=nembrao(q); d=nem(q,c); e=nem1(q,c);
printf("%6.3lf %9.5lf %8.5lf %9.5lf %8.5lf %19.15lf %18.15lf\n",q,Re(c),Im(c), Re(d),Im(d), Re(e),Im(e) );}

//FILE *o; o=fopen("nembraplo.eps","w"); ado(o,154,604);
FILE *o; o=fopen("nembran.eps","w"); ado(o,154,404);
fprintf(o,"2 2 translate 1000 1000 scale 2 setlinecap\n");
#define M(x,y) fprintf(o,"%8.6f %8.6f M\n",(0.+x),(0.+y));
#define L(x,y) fprintf(o,"%8.6f %8.6f L\n",(0.+x),(0.+y));
#define o(x,y) fprintf(o,"%8.6f %8.6f o\n",(0.+x),(0.+y));

for(n=0;n<41;n+=5){ M(0,.01*n) L(.15,.01*n) }
for(n=0;n<31;n+=5){ M(.01*n,0) L(.01*n,.4) }
fprintf(o,".0006 W S 1 setlinejoin\n");
fprintf(o,"1 0 0 RGB .001 W\n");

for(n= 1;n<10;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
//for(n=11;n<20;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
//for(n=21;n<30;n++){ q=.1*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
fprintf(o,"0 .5 0 RGB\n");
for(n= 1;n<10;n++){ q=1.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
fprintf(o,"0 0 1 RGB\n");
for(n= 1;n<10;n++){ q=10.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
fprintf(o,".5 0 .8 RGB\n");
for(n= 1;n<10;n++){ q=100.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
fprintf(o,"0 0 0 RGB\n");
for(n= 0;n<11;n+=10){ q=100.*n; c=nembran(q); x=Re(c); y=-Im(c); o(x,y); }
                      q=10000.; c=nembran(q); x=Re(c); y=-Im(c); o(x,y);

for(n= 0;n<400;n++){ q=-.006+.04*(n+.01*n*n+.001*n*n*n); c=nembran(q); x=Re(c); y=-Im(c); if(n==0) M(x,y) else L(x,y);
printf("%9.3lf %16.6lf %16.6lf\n",q,x,y);
}
L(0,0)
fprintf(o,"0 0 0 RGB .001 W 0 setlinecap S\n");

fprintf(o,"showpage\n%c%cTrailer\n",'%','%');
fclose(o);

system("epstopdf nembran.eps");
system("open nembran.pdf");
getchar();
return 0;}

Latex generator of labels


\documentclass[12pt]{article}
\usepackage{geometry}
\paperwidth 174pt
\paperheight 420pt
\topmargin -100pt
\oddsidemargin -72pt
\textheight 800px
\parindent 0pt
\usepackage{graphicx}
\usepackage{rotating}
\newcommand \ing {\includegraphics}
\newcommand \sx {\scalebox}
\begin{document}
\begin{picture}(170,410)
%\put(20,14){\ing{nembraplo}}
\put(20,14){\ing{nembran}}
\put(8,411){\sx{1.4}{$y$}}
%\put(0,512){\sx{1.4}{$0.5$}}
%\put(0,412){\sx{1.4}{$0.4$}}
\put(0,312){\sx{1.4}{$0.3$}}
\put(0,212){\sx{1.4}{$0.2$}}
\put(0,112){\sx{1.4}{$0.1$}}
\put(8,12){\sx{1.4}{$0$}}
\put(20,1){\sx{1.4}{$0$}}
\put(60,1){\sx{1.4}{$0.05$}}
\put(112,1){\sx{1.4}{$0.1$}}
\put(162,2){\sx{1.4}{$x$}}
%\put(18,601){\sx{1.1}{$q\!=\!0$}}
\put(0,389){\sx{1.1}{$q\!=\!0$}}
%\put(63,593){\sx{1.1}{$q\!=\!0.2$}}
\put(43,404){\sx{1.1}{$q\!=\!0.2$}}
%\put(97,574){\sx{1.1}{$q\!=\!0.4$}}
\put(62,394){\sx{1.1}{$q\!=\!0.4$}}
\put(78,384){\sx{1.1}{$q\!=\!0.6$}}
\put(88,372){\sx{1.1}{$q\!=\!0.8$}}
\put(97,356){\sx{1.1}{$q\!=\!1$}}
\put(116,312){\sx{1.1}{$q\!=\!2$}}
\put(120,284){\sx{1.1}{$q\!=\!3$}}
\put(122,263){\sx{1.1}{$q\!=\!4$}}
\put(122,250){\sx{1.1}{$q\!=\!5$}}
\put(112,200){\sx{1.1}{$q\!=\!10$}}
\put(104,164){\sx{1.1}{$q\!=\!20$}}
\put(94,145){\sx{1.1}{$q\!=\!30$}}
\put( 74,99){\sx{1.1}{$q\!=\!100$}}
\put( 49,52){\sx{1.1}{$q\!=\!1000$}}
\put( 36,32){\sx{1.1}{$q\!=\!10000$}}
\end{picture}
\end{document}

File history

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

Date/TimeThumbnailDimensionsUserComment
current06:13, 1 December 2018Thumbnail for version as of 06:13, 1 December 2018361 × 871 (65 KB)Maintenance script (talk | contribs)Importing image file
  • You cannot overwrite this file.

The following 3 pages link to this file:

Metadata