File:Sunept.jpg

#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"

#include "ado.cin"
#define M(x,y) fprintf(o,"%8.4lf %8.4lf M\n",0.+x,0.+y);
#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!

#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;
}

\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}