File:Logi5T1500.png

#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 "ado.cin"
#include "efjh.cin"
int main(){ int j,k,m,n; DB x,y, p,q, t; z_type z,c,d;
FILE *o;o=fopen("logi6.eps","w");ado(o,164,24);
fprintf(o,"62 2 translate\n 20 20 scale\n");
#define M(x,y) fprintf(o,"%6.3f %6.3f M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%6.3f %6.3f L\n",0.+x,0.+y);

fprintf(o,"1 setlinejoin 2 setlinecap\n");
for(m=-3;m<6;m++){if(m==0){M(m,-.04)L(m,1.06)} else{M(m,0)L(m,1)}}
for(n=0;n<2;n++){ M( -3,n)L(5,n)}
fprintf(o,".008 W 0 0 0 RGB S\n");

maq(3.4);
DO(m,1004) { x=-3.+8.*sqrt(.001*m); y=Re(F(x)); if(m==0)M(x,y) else L(x,y);}
fprintf(o,".01 W 0 .7 0 RGB S\n");

maq(3.);
DO(m,1004) { x=-3.+8.*sqrt(.001*m); y=Re(F(x)); if(m==0)M(x,y) else L(x,y);}
fprintf(o,".015 W 1 0 0 RGB [.03 .04] 0 setdash S\n");

fprintf(o,"1 setlinejoin 1 setlinecap\n");

maq(3.8);
DO(m,1004) { x=-3.+8.*sqrt(.001*m); y=Re(F(x)); if(y>-2) { if(m==0)M(x,y) else L(x,y);} }
fprintf(o,".015 W 0 0 1 RGB [.001 .025] 0 setdash S\n");

fprintf(o,"showpage\n%c%cTrailer",'%','%'); fclose(o);
        system("epstopdf logi6.eps");
        system( "open logi6.pdf");
        getchar(); system("killall Preview");
}