File:Straroplot.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"

DB Q=.5;

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

z_type F1(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 ;
return sqrt(r);
}

int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d, cu,cd;
FILE *o;o=fopen("straroplo.eps","w"); ado(o,604,724);
//FILE *o;o=fopen("05.eps","w"); ado(o,604,724);
fprintf(o,"302 122 translate\n 100 100 scale\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
for(m=-30;m<31;m+=10){M(.1*m,-1.2) L(.1*m,6) }
fprintf(o,".003 W 0 0 0 RGB S\n");

for(n=-12;n<61;n+=1){M(-3,.1*n) L(3,.1*n)}
fprintf(o,".0006 W 0 0 0 RGB S\n");

for(n=-10;n<61;n+=10){M(-3,.1*n) L(3,.1*n)}
fprintf(o,".01 W 0 0 0 RGB S 0 setlinecap\n");

Q=0.5; S=nembran(Q);
DO(m,301){x=.01*(m+.1); y=.1*Re(F(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".02 W 1 0 0 RGB S\n");

DO(m,301){x=-.01*(m); y=.1*Re(F(x+.001*I)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break;
//printf("%8.2lf %8.2lf\n",x,y);
} fprintf(o,".02 W 1 0 0 RGB S\n");

DO(m,301){x=-3.+.01*(m+.1); y=.1*Re(F1(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".001 W 0 0 0 RGB S\n");

Q=1.; S=nembran(Q);
DO(m,301){x=.01*(m+.1); y=.1*Re(F(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".02 W 0 .8 0 RGB S\n");
DO(m,301){x=-.01*(m); y=.1*Re(F(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".02 W 0 .8 0 RGB S\n");
DO(m,301){x=-3.+.01*(m+.1); y=.1*Re(F1(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".001 W 0 0 0 RGB S\n");

Q=2.; S=nembran(Q);
DO(m,301){x=.01*(m+.1); y=.1*Re(F(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".02 W 0 0 1 RGB S\n");
DO(m,301){x=-.01*(m); y=.1*Re(F(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".02 W 0 0 1 RGB S\n");
DO(m,301){x=-3.+.01*(m+.1); y=.1*Re(F1(x)); if(m==0) M(x,y) else L(x,y) ; if(y>8) break; } fprintf(o,".001 W 0 0 0 RGB S\n");

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

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

\documentclass[12pt]{article}
\usepackage{geometry}
\paperwidth 606pt
\paperheight 726pt
\topmargin -102pt
\oddsidemargin -72pt
\textheight 800px
\parindent 0pt
\usepackage{graphicx}
\usepackage{rotating}
\newcommand \ing {\includegraphics}
\newcommand \sx {\scalebox}
\begin{document}
\begin{picture}(170,718)
%\put(20,14){\ing{nembraplo}}
\put(277,706){\sx{3.3}{$y$}}
\put(260,612){\sx{3.3}{$50$}}
\put(260,511){\sx{3.3}{$40$}}
\put(260,411){\sx{3.3}{$30$}}
\put(260,311){\sx{3.3}{$20$}}
\put(260,211){\sx{3.3}{$10$}}
\put(274,111){\sx{3.3}{$0$}}
\put(224, 11){\sx{3.3}{$-10$}}
\put( 74,94){\sx{3.3}{$-2$}}
\put(174,94){\sx{3.3}{$-1$}}
\put(294,94){\sx{3.3}{$0$}}
\put(394,94){\sx{3.3}{$1$}}
\put(494,94){\sx{3.3}{$2$}}
\put(584,95){\sx{3.33}{$x$}}
\put(0,0){\ing{straroplo}}
\end{picture}
\end{document}