Difference between revisions of "File:Straroplot.jpg"

From TORI
Jump to navigation Jump to search
(Importing image file)
 
 
Line 1: Line 1:
  +
Esplicit plot of function [[StraRo]] for various vvalues of parameter:
−
Importing image file
  
  +
  +
$q\!=\!0.5$ , red
  +
  +
$q\!=\!1$ , green
  +
  +
$q\!=\!2$ , blue
  +
  +
Dependences $y=\mathrm{StraRo}_q(x)$ are shown with thick lines.
  +
These lines have jump at zero. Such a jump is unavoidable wile function [[StraRo]] is holomorphic above the upper branch point and below the lower branch point.
  +
However, the jump can be moved away from the real axis in the complex plane; the corresponding extensions are shown with thin lines.
  +
  +
==Refereneces==
  +
<references/>
  +
  +
[[Category:ArqNem]]
  +
[[Category:AuNem]]
  +
[[Category:Book]]
  +
[[Category:BookPlot]]
  +
[[Category:C++]]
  +
[[Category:Latex]]
  +
[[Category:Nemtsov Function]]
  +
  +
==[[C++]] generator of curves==
  +
  +
<poem><nomathjax><nowiki>
  +
  +
#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
  +
}
  +
  +
</nowiki></nomathjax></poem>
  +
==[[Latex]] generator of labels==
  +
  +
<poem><nomathjax><nowiki>
  +
\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}
  +
</nowiki></nomathjax></poem>

Latest revision as of 08:52, 1 December 2018

Esplicit plot of function StraRo for various vvalues of parameter:

$q\!=\!0.5$ , red

$q\!=\!1$ , green

$q\!=\!2$ , blue

Dependences $y=\mathrm{StraRo}_q(x)$ are shown with thick lines. These lines have jump at zero. Such a jump is unavoidable wile function StraRo is holomorphic above the upper branch point and below the lower branch point. However, the jump can be moved away from the real axis in the complex plane; the corresponding extensions are shown with thin lines.

Refereneces


C++ generator of curves



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

Latex generator of labels


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

File history

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

Date/TimeThumbnailDimensionsUserComment
current06:14, 1 December 2018Thumbnail for version as of 06:14, 1 December 20181,257 × 1,506 (265 KB)Maintenance script (talk | contribs)Importing image file

You cannot overwrite this file.

The following page uses this file:

Metadata

This file contains additional information, probably added from the digital camera or scanner used to create or digitize it.

If the file has been modified from its original state, some details may not fully reflect the modified file.

Retrieved from "https://mizugadro.mydns.jp/t/index.php?title=File:Straroplot.jpg&oldid=21689"