Difference between revisions of "File:Ausinsusinmapt50.jpg"

From TORI
Jump to navigation Jump to search
(Importing image file)
 
 
Line 1: Line 1:
  +
Range, where the Abel sin [[AuSin]] is [[inverse function]] of super sin [[SuSin]].
−
Importing image file
  
  +
  +
[[Range of validity]] of relation
  +
  +
$~\mathrm{SuSin}\big(\mathrm{AuSin}(z)\big)\!=\!z~$
  +
  +
in the complex plane $~z\!=\!x\!+\!\mathrm i y~$.
  +
  +
The region of validity of the equation is shaded with the rectangular grid.
  +
  +
The boundary of this range follows the lines $\Im\big(\mathrm{AuSin}(z)\Big)\!=\!0~$;
  +
these lines are also drawn.
  +
  +
==[[C++]] generator of rectangualr grid==
  +
<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 complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
#include "conto.cin"
  +
#include "arcsin.cin"
  +
#include"susin.cin"
  +
#include"ausin.cin"
  +
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
  +
printf("1.-ausin(1)=%19.17f , 1.-ausin(z_type(1.,1.e-14))=%19.17f ;1.- susun(1.)=%19.17f, 1.-susun(z_type(1.,1.e-14))=%19.17f \n",
  +
1.-Re(ausin(1.)), 1.-Re(ausin(z_type(1.,1.e-14))), 1.-Re(susin(1.)), 1.-Re(susin(z_type(1.,1.e-14))) );
  +
// getchar();
  +
  +
int M=548,M1=M+1;
  +
int N=421,N1=N+1;
  +
DB X[M1],Y[N1];
  +
DB *g, *f, *w; // w is working array.
  +
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
char v[M1*N1]; // v is working array
  +
FILE *o;o=fopen("ausinsusinmap1.eps","w"); ado(o,332,402);
  +
fprintf(o,"1 201 translate\n 100 100 scale\n");
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
DB e;
  +
e=1.4/sinh(1.);
  +
DO(m,M1) { t=(m+.5)/400.; X[m]=e*sinh(1.*t);}
  +
e=1.4/sinh(1.);
  +
DO(n,N1) { t=(n-210.5)/210.; Y[n]=e*sinh(1.*t);}
  +
for(m=0;m<4;m+=1){M(m,-2) L(m,2) }
  +
for(n=-2;n<3;n+=1){M(0,n) L(M_PI,n)}
  +
fprintf(o,".002 W 0 0 0 RGB S\n");
  +
M(M_PI/2,-2)L(M_PI/2,2)
  +
M(M_PI,-2)L(M_PI,2)
  +
fprintf(o,".001 W 0 0 0 RGB S\n");
  +
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
  +
f[m*N1+n]=9999;}
  +
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
  +
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)
  +
// c=arcsin(z);
  +
// c=sqrt(3./z);
  +
// c=susin(z);
  +
c=susin(ausin(z));
  +
// d=susin(c);
  +
// d=arcsin(su0(z+1.));
  +
// p=abs(z-d)/(abs(z)+abs(d)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
  +
p=Re(c); q=Im(c); if(p>-101 && p<101 && q>-101 && q<101 && abs(c-z)<.01){ g[m*N1+n]=p;f[m*N1+n]=q;}
  +
}}
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
p=200.;q=1.;
  +
/*
  +
p=9;q=.16;
  +
conto(o,g,w,v,X,Y,M,N,(15.6 ),-p,p); fprintf(o,".001 W .4 1 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(15. ),-p,p); fprintf(o,".005 W 1 0 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(14. ),-p,p); fprintf(o,".002 W .2 .2 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(13. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(12. ),-p,p); fprintf(o,".006 W 0 0 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(11. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(10. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (9. ),-p,p); fprintf(o,".006 W 0 1 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (8. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (7. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (6. ),-p,p); fprintf(o,".006 W 0 .5 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (5. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (4. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (3. ),-p,p); fprintf(o,".006 W 1 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (2. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (1. ),-p,p); fprintf(o,".004 W .5 0 0 RGB S\n");
  +
*/
  +
for(m=-4;m<4;m++)for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q,q);fprintf(o,".006 W 0 .6 0 RGB S\n");
  +
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q,q);fprintf(o,".006 W .9 0 0 RGB S\n");
  +
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q,q);fprintf(o,".006 W 0 0 .9 RGB S\n");
  +
for(m= 1;m<8;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".009 W .8 0 0 RGB S\n");
  +
for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".009 W 0 0 .8 RGB S\n");
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".009 W .5 0 .5 RGB S\n");
  +
for(m=-8;m<9;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".009 W 0 0 0 RGB S\n");
  +
/*
  +
fprintf(o,"0 setlinejoin 0 setlinecap\n");
  +
M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .01 W S\n");
  +
DO(n,51) {x=-.2*n; M(x-.01,0) L(x-.09,0) } fprintf(o,"0 0 0 RGB .03 W S\n");
  +
//#include "plofu.cin"
  +
*/
  +
  +
fprintf(o,"showpage\n");
  +
fprintf(o,"%c%cTrailer\n",'%','%');
  +
fclose(o); free(f); free(g); free(w);
  +
system("epstopdf ausinsusinmap1.eps");
  +
system( "open ausinsusinmap1.pdf"); //for macintosh
  +
getchar(); system("killall Preview"); // For macintosh
  +
}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[C++]] generator of boundaries==
  +
<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 complex<double> z_type;
  +
#define Re(x) x.real()
  +
#define Im(x) x.imag()
  +
#define I z_type(0.,1.)
  +
#include "conto.cin"
  +
#include "arcsin.cin"
  +
#include"susin.cin"
  +
#include"ausin.cin"
  +
  +
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
  +
printf("1.-ausin(1)=%19.17f , 1.-ausin(z_type(1.,1.e-14))=%19.17f ;1.- susun(1.)=%19.17f, 1.-susun(z_type(1.,1.e-14))=%19.17f \n",
  +
1.-Re(ausin(1.)), 1.-Re(ausin(z_type(1.,1.e-14))), 1.-Re(susin(1.)), 1.-Re(susin(z_type(1.,1.e-14))) );
  +
// getchar();
  +
  +
int M=548,M1=M+1;
  +
int N=421,N1=N+1;
  +
DB X[M1],Y[N1];
  +
DB *g, *f, *w; // w is working array.
  +
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
  +
char v[M1*N1]; // v is working array
  +
FILE *o;o=fopen("ausinsusinmap2.eps","w"); ado(o,332,402);
  +
fprintf(o,"1 201 translate\n 100 100 scale\n");
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
DB e;
  +
e=2./sinh(1.);
  +
DO(m,M1) { t=(m+.5)/400.; X[m]=e*sinh(1.*t);}
  +
e=2./sinh(1.);
  +
DO(n,N1) { t=(n-210.5)/210.; Y[n]=e*sinh(1.*t);}
  +
for(m=0;m<4;m+=1){M(m,-2) L(m,2) }
  +
for(n=-2;n<3;n+=1){M(0,n) L(M_PI,n)}
  +
fprintf(o,".0004 W 0 0 0 RGB S\n");
  +
M(M_PI/2,-2)L(M_PI/2,2)
  +
M(M_PI,-2)L(M_PI,2)
  +
fprintf(o,".0003 W 0 0 0 RGB S\n");
  +
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
  +
f[m*N1+n]=9999;}
  +
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
  +
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)
  +
// c=arcsin(z);
  +
// c=sqrt(3./z);
  +
c=ausin(z);
  +
// c=susin(ausin(z));
  +
// d=susin(c);
  +
// d=arcsin(su0(z+1.));
  +
// p=abs(z-d)/(abs(z)+abs(d)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
  +
p=Re(c); q=Im(c); if(p>-101 && p<101 && q>-101 && q<101 ){ g[m*N1+n]=p;f[m*N1+n]=q;}
  +
}}
  +
fprintf(o,"1 setlinejoin 2 setlinecap\n");
  +
p=200.;q=1.;
  +
/*
  +
p=9;q=.16;
  +
conto(o,g,w,v,X,Y,M,N,(15.6 ),-p,p); fprintf(o,".001 W .4 1 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(15. ),-p,p); fprintf(o,".005 W 1 0 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(14. ),-p,p); fprintf(o,".002 W .2 .2 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(13. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(12. ),-p,p); fprintf(o,".006 W 0 0 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(11. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N,(10. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (9. ),-p,p); fprintf(o,".006 W 0 1 1 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (8. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (7. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (6. ),-p,p); fprintf(o,".006 W 0 .5 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (5. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (4. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (3. ),-p,p); fprintf(o,".006 W 1 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (2. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
  +
conto(o,g,w,v,X,Y,M,N, (1. ),-p,p); fprintf(o,".004 W .5 0 0 RGB S\n");
  +
*/
  +
/*
  +
for(m=-4;m<4;m++)for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q,q);fprintf(o,".0016 W 0 .6 0 RGB S\n");
  +
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q,q);fprintf(o,".0016 W .9 0 0 RGB S\n");
  +
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q,q);fprintf(o,".0016 W 0 0 .9 RGB S\n");
  +
for(m= 1;m<8;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".004 W .8 0 0 RGB S\n");
  +
for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".004 W 0 0 .8 RGB S\n");
  +
*/
  +
conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".012 W .5 0 .5 RGB S\n");
  +
/*
  +
for(m=-8;m<9;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".004 W 0 0 0 RGB S\n");
  +
fprintf(o,"0 setlinejoin 0 setlinecap\n");
  +
M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .01 W S\n");
  +
DO(n,51) {x=-.2*n; M(x-.01,0) L(x-.09,0) } fprintf(o,"0 0 0 RGB .03 W S\n");
  +
//#include "plofu.cin"
  +
*/
  +
  +
fprintf(o,"showpage\n");
  +
fprintf(o,"%c%cTrailer\n",'%','%');
  +
fclose(o); free(f); free(g); free(w);
  +
system("epstopdf ausinsusinmap2.eps");
  +
system( "open ausinsusinmap2.pdf"); //for macintosh
  +
getchar(); system("killall Preview"); // For macintosh
  +
}
  +
</nowiki></nomathjax></poem>
  +
  +
==[[Latex]] generator of labels==
  +
<poem><nomathjax><nowiki>
  +
  +
\documentclass[12pt]{article}
  +
\usepackage{geometry}
  +
\usepackage{graphics}
  +
\usepackage{rotating}
  +
\paperwidth 3288pt
  +
\paperheight 4150pt
  +
\topmargin -100pt
  +
\oddsidemargin -72pt
  +
\textwidth 3200pt
  +
\textheight 4200pt
  +
\newcommand \sx {\scalebox}
  +
\newcommand \rot {\begin{rotate}}
  +
\newcommand \ero {\end{rotate}}
  +
\pagestyle{empty}
  +
\begin{document}
  +
\sx{10}{\begin{picture}(328,412)
  +
\put(10,9){\includegraphics{ausinsusinmap1}}
  +
\put(10,9){\includegraphics{ausinsusinmap2}}
  +
\put(2,406){\sx{1.2}{$y$}}
  +
\put(2,306){\sx{1.2}{$1$}}
  +
\put(2,206){\sx{1.2}{$0$}}
  +
\put(-7,106){\sx{1.2}{$-1$}}
  +
%\put(-7, 06){\sx{1.2}{$-2$}}
  +
\put(9,-1){\sx{1.2}{$0$}}
  +
\put(109,-1){\sx{1.2}{$1$}}
  +
\put(157,-1){\sx{1.2}{$\pi/2$}}
  +
\put(209,-1){\sx{1.2}{$2$}}
  +
\put(309,-1){\sx{1.2}{$3$}}
  +
\put(319,-1){\sx{1.2}{$x$}}
  +
\end{picture}}
  +
\end{document}
  +
</nowiki></nomathjax></poem>
  +
  +
  +
[[Category:AuSin]]
  +
[[Category:Book]]
  +
[[Category:BookMap]]
  +
[[Category:SuSin]]
  +
[[Category:Inverse function]]
  +
[[Category:Complex map]]
  +
[[Category:C++]]
  +
[[Category:Latex]]
  +
[[Category:Book]]

Latest revision as of 08:30, 1 December 2018

Range, where the Abel sin AuSin is inverse function of super sin SuSin.

Range of validity of relation

$~\mathrm{SuSin}\big(\mathrm{AuSin}(z)\big)\!=\!z~$

in the complex plane $~z\!=\!x\!+\!\mathrm i y~$.

The region of validity of the equation is shaded with the rectangular grid.

The boundary of this range follows the lines $\Im\big(\mathrm{AuSin}(z)\Big)\!=\!0~$; these lines are also drawn.

C++ generator of rectangualr grid


#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 complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
#include "conto.cin"
#include "arcsin.cin"
#include"susin.cin"
#include"ausin.cin"

int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
printf("1.-ausin(1)=%19.17f , 1.-ausin(z_type(1.,1.e-14))=%19.17f ;1.- susun(1.)=%19.17f, 1.-susun(z_type(1.,1.e-14))=%19.17f \n",
     1.-Re(ausin(1.)), 1.-Re(ausin(z_type(1.,1.e-14))), 1.-Re(susin(1.)), 1.-Re(susin(z_type(1.,1.e-14))) );
// getchar();

int M=548,M1=M+1;
int N=421,N1=N+1;
DB X[M1],Y[N1];
DB *g, *f, *w; // w is working array.
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
char v[M1*N1]; // v is working array
FILE *o;o=fopen("ausinsusinmap1.eps","w"); ado(o,332,402);
fprintf(o,"1 201 translate\n 100 100 scale\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
DB e;
e=1.4/sinh(1.);
DO(m,M1) { t=(m+.5)/400.; X[m]=e*sinh(1.*t);}
e=1.4/sinh(1.);
DO(n,N1) { t=(n-210.5)/210.; Y[n]=e*sinh(1.*t);}
for(m=0;m<4;m+=1){M(m,-2) L(m,2) }
for(n=-2;n<3;n+=1){M(0,n) L(M_PI,n)}
 fprintf(o,".002 W 0 0 0 RGB S\n");
M(M_PI/2,-2)L(M_PI/2,2)
M(M_PI,-2)L(M_PI,2)
 fprintf(o,".001 W 0 0 0 RGB S\n");
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
                       f[m*N1+n]=9999;}
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)
// c=arcsin(z);
// c=sqrt(3./z);
// c=susin(z);
   c=susin(ausin(z));
// d=susin(c);
// d=arcsin(su0(z+1.));
// p=abs(z-d)/(abs(z)+abs(d)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
 p=Re(c); q=Im(c); if(p>-101 && p<101 && q>-101 && q<101 && abs(c-z)<.01){ g[m*N1+n]=p;f[m*N1+n]=q;}
       }}
fprintf(o,"1 setlinejoin 2 setlinecap\n");
 p=200.;q=1.;
/*
  p=9;q=.16;
 conto(o,g,w,v,X,Y,M,N,(15.6 ),-p,p); fprintf(o,".001 W .4 1 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(15. ),-p,p); fprintf(o,".005 W 1 0 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(14. ),-p,p); fprintf(o,".002 W .2 .2 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(13. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(12. ),-p,p); fprintf(o,".006 W 0 0 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(11. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(10. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (9. ),-p,p); fprintf(o,".006 W 0 1 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (8. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (7. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (6. ),-p,p); fprintf(o,".006 W 0 .5 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (5. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (4. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (3. ),-p,p); fprintf(o,".006 W 1 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (2. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (1. ),-p,p); fprintf(o,".004 W .5 0 0 RGB S\n");
*/
for(m=-4;m<4;m++)for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q,q);fprintf(o,".006 W 0 .6 0 RGB S\n");
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q,q);fprintf(o,".006 W .9 0 0 RGB S\n");
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q,q);fprintf(o,".006 W 0 0 .9 RGB S\n");
for(m= 1;m<8;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".009 W .8 0 0 RGB S\n");
for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".009 W 0 0 .8 RGB S\n");
               conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".009 W .5 0 .5 RGB S\n");
for(m=-8;m<9;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".009 W 0 0 0 RGB S\n");
/*
fprintf(o,"0 setlinejoin 0 setlinecap\n");
M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .01 W S\n");
DO(n,51) {x=-.2*n; M(x-.01,0) L(x-.09,0) } fprintf(o,"0 0 0 RGB .03 W S\n");
//#include "plofu.cin"
*/

fprintf(o,"showpage\n");
fprintf(o,"%c%cTrailer\n",'%','%');
fclose(o); free(f); free(g); free(w);
      system("epstopdf ausinsusinmap1.eps");
      system( "open ausinsusinmap1.pdf"); //for macintosh
      getchar(); system("killall Preview"); // For macintosh
}

C++ generator of boundaries



#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 complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
#include "conto.cin"
#include "arcsin.cin"
#include"susin.cin"
#include"ausin.cin"

int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
printf("1.-ausin(1)=%19.17f , 1.-ausin(z_type(1.,1.e-14))=%19.17f ;1.- susun(1.)=%19.17f, 1.-susun(z_type(1.,1.e-14))=%19.17f \n",
     1.-Re(ausin(1.)), 1.-Re(ausin(z_type(1.,1.e-14))), 1.-Re(susin(1.)), 1.-Re(susin(z_type(1.,1.e-14))) );
// getchar();

int M=548,M1=M+1;
int N=421,N1=N+1;
DB X[M1],Y[N1];
DB *g, *f, *w; // w is working array.
g=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
f=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
w=(DB *)malloc((size_t)((M1*N1)*sizeof(DB)));
char v[M1*N1]; // v is working array
FILE *o;o=fopen("ausinsusinmap2.eps","w"); ado(o,332,402);
fprintf(o,"1 201 translate\n 100 100 scale\n");
fprintf(o,"1 setlinejoin 2 setlinecap\n");
DB e;
e=2./sinh(1.);
DO(m,M1) { t=(m+.5)/400.; X[m]=e*sinh(1.*t);}
e=2./sinh(1.);
DO(n,N1) { t=(n-210.5)/210.; Y[n]=e*sinh(1.*t);}
for(m=0;m<4;m+=1){M(m,-2) L(m,2) }
for(n=-2;n<3;n+=1){M(0,n) L(M_PI,n)}
 fprintf(o,".0004 W 0 0 0 RGB S\n");
M(M_PI/2,-2)L(M_PI/2,2)
M(M_PI,-2)L(M_PI,2)
 fprintf(o,".0003 W 0 0 0 RGB S\n");
DO(m,M1)DO(n,N1){ g[m*N1+n]=9999;
                       f[m*N1+n]=9999;}
DO(m,M1){x=X[m]; if(m/10*10==m) printf("x=%6.3f\n",x);
DO(n,N1){y=Y[n]; z=z_type(x,y); //if(abs(z+2.)>.019)
// c=arcsin(z);
// c=sqrt(3./z);
 c=ausin(z);
// c=susin(ausin(z));
// d=susin(c);
// d=arcsin(su0(z+1.));
// p=abs(z-d)/(abs(z)+abs(d)); p=-log(p)/log(10.); if(p>0 && p<17) g[m*N1+n]=p;
 p=Re(c); q=Im(c); if(p>-101 && p<101 && q>-101 && q<101 ){ g[m*N1+n]=p;f[m*N1+n]=q;}
       }}
fprintf(o,"1 setlinejoin 2 setlinecap\n");
 p=200.;q=1.;
/*
  p=9;q=.16;
 conto(o,g,w,v,X,Y,M,N,(15.6 ),-p,p); fprintf(o,".001 W .4 1 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(15. ),-p,p); fprintf(o,".005 W 1 0 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(14. ),-p,p); fprintf(o,".002 W .2 .2 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(13. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(12. ),-p,p); fprintf(o,".006 W 0 0 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(11. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N,(10. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (9. ),-p,p); fprintf(o,".006 W 0 1 1 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (8. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (7. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (6. ),-p,p); fprintf(o,".006 W 0 .5 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (5. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (4. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (3. ),-p,p); fprintf(o,".006 W 1 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (2. ),-p,p); fprintf(o,".002 W 0 0 0 RGB S\n");
 conto(o,g,w,v,X,Y,M,N, (1. ),-p,p); fprintf(o,".004 W .5 0 0 RGB S\n");
*/
/*
for(m=-4;m<4;m++)for(n=1;n<10;n+=1)conto(o,f,w,v,X,Y,M,N,(m+.1*n),-q,q);fprintf(o,".0016 W 0 .6 0 RGB S\n");
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N,-(m+.1*n),-q,q);fprintf(o,".0016 W .9 0 0 RGB S\n");
for(m=0;m<4;m++) for(n=1;n<10;n+=1)conto(o,g,w,v,X,Y,M,N, (m+.1*n),-q,q);fprintf(o,".0016 W 0 0 .9 RGB S\n");
for(m= 1;m<8;m++) conto(o,f,w,v,X,Y,M,N, (0.-m),-p,p);fprintf(o,".004 W .8 0 0 RGB S\n");
for(m= 1;m<9;m++) conto(o,f,w,v,X,Y,M,N, (0.+m),-p,p);fprintf(o,".004 W 0 0 .8 RGB S\n");
*/
               conto(o,f,w,v,X,Y,M,N, (0. ),-p,p); fprintf(o,".012 W .5 0 .5 RGB S\n");
/*
for(m=-8;m<9;m++)conto(o,g,w,v,X,Y,M,N,(0.+m),-p,p);fprintf(o,".004 W 0 0 0 RGB S\n");
fprintf(o,"0 setlinejoin 0 setlinecap\n");
M(-10,0)L(0,0) fprintf(o,"1 1 1 RGB .01 W S\n");
DO(n,51) {x=-.2*n; M(x-.01,0) L(x-.09,0) } fprintf(o,"0 0 0 RGB .03 W S\n");
//#include "plofu.cin"
*/

fprintf(o,"showpage\n");
fprintf(o,"%c%cTrailer\n",'%','%');
fclose(o); free(f); free(g); free(w);
      system("epstopdf ausinsusinmap2.eps");
      system( "open ausinsusinmap2.pdf"); //for macintosh
      getchar(); system("killall Preview"); // For macintosh
}

Latex generator of labels



\documentclass[12pt]{article}
\usepackage{geometry}
\usepackage{graphics}
\usepackage{rotating}
\paperwidth 3288pt
\paperheight 4150pt
\topmargin -100pt
\oddsidemargin -72pt
\textwidth 3200pt
\textheight 4200pt
\newcommand \sx {\scalebox}
\newcommand \rot {\begin{rotate}}
\newcommand \ero {\end{rotate}}
\pagestyle{empty}
\begin{document}
\sx{10}{\begin{picture}(328,412)
\put(10,9){\includegraphics{ausinsusinmap1}}
\put(10,9){\includegraphics{ausinsusinmap2}}
\put(2,406){\sx{1.2}{$y$}}
\put(2,306){\sx{1.2}{$1$}}
\put(2,206){\sx{1.2}{$0$}}
\put(-7,106){\sx{1.2}{$-1$}}
%\put(-7, 06){\sx{1.2}{$-2$}}
\put(9,-1){\sx{1.2}{$0$}}
\put(109,-1){\sx{1.2}{$1$}}
\put(157,-1){\sx{1.2}{$\pi/2$}}
\put(209,-1){\sx{1.2}{$2$}}
\put(309,-1){\sx{1.2}{$3$}}
\put(319,-1){\sx{1.2}{$x$}}
\end{picture}}
\end{document}

File history

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

Date/TimeThumbnailDimensionsUserComment
current06:10, 1 December 2018Thumbnail for version as of 06:10, 1 December 20182,274 × 2,871 (1,015 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:Ausinsusinmapt50.jpg&oldid=18292"