File:Sunept.jpg

Explicit plot of function SuNe:

$y=\mathrm{SuNe}_q(x)$

for $q=-1$, $q=-0.5$, $q=0$, $q=0.5$, $q=1$, $q=2$, $q=3$.

C++ generator of curves
//using namespace std; typedef std::complex z_type; //#include "conto.cin"
 * 1) include 
 * 2) include 
 * 3) include 
 * 4) define DB double
 * 5) define DO(x,y) for(x=0;x<y;x++)
 * 1) include
 * 1) define Re(x) x.real
 * 2) define Im(x) x.imag
 * 3) define I z_type(0.,1.)


 * 1) include "ado.cin"
 * 2) define M(x,y) fprintf(o,"%8.4lf %8.4lf M\n",0.+x,0.+y);
 * 3) define L(x,y) fprintf(o,"%8.4lf %8.4lf L\n",0.+x,0.+y);

DB Q; z_type nem(z_type z){ return z*(1.+z*z*(1.+z*Q)); } z_type nem1(z_type z){ return 1.+z*z*(3.+z*(4.*Q)); } // WARNING: Q is global!


 * 1) include "sune.cin"

DB SUNo=0;

z_type sunem(z_type z){ return sune(z + SUNo);}

int main{ int Max; int j,k,m,n; DB x,y, p,q, t; z_type z,c,d; // DB rr,ti; //FILE *o;o=fopen("suneplo4.eps","w");ado(o,420,420); FILE *o;o=fopen("sunep.eps","w");ado(o,1120,1020); fprintf(o,"1010 10 translate\n 100 100 scale 2 setlinecap\n"); for(m=-10;m<1;m+=1){ M(m,0)L(m,10)} for(n=0;n<11;n+=1){ M(-10,n)L(0,n)} fprintf(o,".008 W 0 0 0 RGB 2 setlinecap S\n");

Q=-1.; K=Q*Q; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".005 W 0 0 0 RGB 1 setlinejoin S\n");

Q=-.5; K=Q*Q; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".005 W 0 0 0 RGB 1 setlinejoin S\n");

Q=0; K=Q*Q; co; printf("Q=%9.4lf\n",Q); //DO(n,50){ y=Re(sunem(0)); SUNo-= y-1.; printf("%19.16lf %19.16lf\n", SUNo,y);} // getchar; DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".01 W 0 0 0 RGB 1 setlinejoin S\n");

Q=.5; K=Q*Q; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".01 W .7 0 .6 RGB 1 setlinejoin S\n");

Q=1; K=Q*Q; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".012 W .8 0 0 RGB 1 setlinejoin S\n");

Q=2; K=Q*Q; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".012 W 0 .6 0 RGB 1 setlinejoin S\n");

Q=3; K=Q*Q; SUNo=-3.; co; printf("Q=%9.4lf\n",Q); DO(n,1112){x=-10.1+.01*n; y=10.*Re(sune(10.*x)); if(n==0) M(x,y) else L(x,y); printf("%6.1lf %18.15lf\n", x,y); if(y>10.1) break;} fprintf(o,".012 W 0 0 .9 RGB 1 setlinejoin S\n");

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

system("epstopdf sunep.eps"); system( "open sunep.pdf"); //mac return 0; }

Latex generator of labels
\documentclass{mcom-l} % \documentclass[12pt]{article} \usepackage{graphics} \paperwidth 1098pt \paperheight 1010pt \usepackage{geometry} \usepackage{rotating} \textwidth 2560pt \textheight 2260pt \topmargin -98pt \oddsidemargin -84pt \parindent 0pt \pagestyle{empty} \newcommand \ing {\includegraphics} \newcommand \sx {\scalebox} \newcommand \rot {\begin{rotate}} \newcommand \ero {\end{rotate}} \begin{document}

\begin{picture}(600,1020) \put(4,0){\ing{sunep}} \put(1044,980){\sx{3.4}{$y$}} \put(1036,900){\sx{3.4}{$0.9$}} \put(1036,800){\sx{3.4}{$0.8$}} \put(1034,700){\sx{3.4}{$0.7$}} \put(1034,600){\sx{3.4}{$0.6$}} \put(1034,500){\sx{3.4}{$0.5$}} \put(1034,400){\sx{3.4}{$0.4$}} \put(1034,300){\sx{3.4}{$0.3$}} \put(1034,200){\sx{3.4}{$0.2$}} \put(1034,100){\sx{3.4}{$0.1$}} %\put(420, -7){\sx{3.3}{$0$}} \put( 82,14){\sx{3.3}{$-90$}} \put(182,14){\sx{3.3}{$-80$}} \put(282,14){\sx{3.3}{$-70$}} \put(382,14){\sx{3.3}{$-60$}} \put(482,14){\sx{3.3}{$-50$}} \put(582,14){\sx{3.3}{$-40$}} \put(682,14){\sx{3.3}{$-30$}} \put(782,14){\sx{3.3}{$-20$}} \put(882,14){\sx{3.3}{$-10$}} \put(1010,14){\sx{3.3}{$0$}} %\put(708,14){\sx{3.3}{$1$}} \put(1090,14){\sx{3.4}{$x$}} \put(988,926){\sx{3.2}{\rot{88}$q\!=\!3$\ero}}% \put(772,194){\sx{3.2}{\rot{30}$q\!=\!-1$\ero}}% \put(882,209){\sx{3.2}{\rot{44}$q\!=\!0$\ero}}% \put(924,176){\sx{3.2}{\rot{44}$q\!=\!2$\ero}}% \end{picture} \end{document}