File:Frac2zt.jpg

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define DO(x,y) for(x=0;x<y;x++)
#define DB double

#include"ado.cin"

DB A=2.;
DB T(DB n,DB x){return A*tan( (M_PI/2)*n+ atan(x/A));}

main(){ FILE *o; int m,n,k; DB x,y,t;
o=fopen("frac2z.eps","w");
ado(o,1002,1002);
#define M(x,y) fprintf(o,"%7.4f %7.4f M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%7.4f %7.4f L\n",0.+x,0.+y);
fprintf(o,"501 501 translate 100 100 scale 2 setlinecap\n");
for(n=-5;n<6;n++) { M(-5,n)L(5,n)}
for(m=-5;m<6;m++) { M(m,-5)L(m,5)}
fprintf(o,".005 W S\n");
M(-5,0)L(5,0)
M(0,-5)L(0,5)
fprintf(o,".02 W S\n");
n=0;DO(m,1001){x=-5.+.01*(m-.5);y=T(2.,x);if(y>-5&&y<5){ if(n==0){M(x,y) n=1;}else L(x,y)} else n=0;} fprintf(o,".03 W 0 1 0 RGB S\n");
n=0;DO(m,1001){x=-5.+.01*(m-.5);y=T(1.,x);if(y>-5&&y<5){ if(n==0){M(x,y) n=1;}else L(x,y)} else n=0;} fprintf(o,".03 W 0 1 0 RGB S\n");
DO(k,20){ t=-2.+.1*k; n=0;DO(m,1001){x=-5.+.01*(m-.5);y=T(t,x);if(y>-7.2&&y<7.2){ if(n==0){M(x,y) n=1;}else L(x,y)} else n=0;} fprintf(o,".01 W 0 0 0 RGB S\n");}
fprintf(o,"showpage\n"); fprintf(o,"%c%cTrailer\n",'%','%');
fclose(o);
system("epstopdf frac2z.eps");
system( "open frac2z.pdf");
}

<nomathjax>
% File Frac2z.pdf should be generated with the code above in order to compile the [[Latex]] document below.


</nomathjax>