File:Nembrant.jpg

From TORI
Revision as of 08:44, 1 December 2018 by Maintenance script (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Original file(361 × 871 pixels, file size: 65 KB, MIME type: image/jpeg)

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

The following 3 pages use this file:

Metadata