File:Tet10bxr.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 std::complex<double> z_type;
#define Re(x) x.real()
#define Im(x) x.imag()
#define I z_type(0.,1.)
//b=10
#include "f4ten.cin"
#include "fit1.cin"
void ado(FILE *O, int X, int Y)
{ fprintf(O,"%c!PS-Adobe-2.0 EPSF-2.0\n",'%');
       fprintf(O,"%c%cBoundingBox: 0 0 %d %d\n",'%','%',X,Y);
       fprintf(O,"/M {moveto} bind def\n");
       fprintf(O,"/L {lineto} bind def\n");
       fprintf(O,"/S {stroke} bind def\n");
       fprintf(O,"/s {show newpath} bind def\n");
       fprintf(O,"/C {closepath} bind def\n");
       fprintf(O,"/F {fill} bind def\n");
       fprintf(O,"/o {.1 0 360 arc C S} bind def\n");
       fprintf(O,"/times-Roman findfont 20 scalefont setfont\n");
       fprintf(O,"/W {setlinewidth} bind def\n");
       fprintf(O,"/RGB {setrgbcolor} bind def\n");}
 //#include "ado.cin"
#define M(x,y) fprintf(o,"%6.4f %6.4f M\n",0.+x,0.+y);
#define L(x,y) fprintf(o,"%6.4f %6.4f L\n",0.+x,0.+y);
#define o(x,y) fprintf(o,"%6.4f %6.4f o\n",0.+x,0.+y);
int main(){ int j,k,m,n; DB p,q,t1,t3,u,v,w,x,y; z_type z,c,d;
FILE *o;o=fopen("tet10bx.eps","w");ado(o,134,124);
fprintf(o,"22 22 translate\n 10 10 scale\n");
fprintf(o,"2 setlinecap\n");
for(m=-2;m<12;m++){if(m!=0){M(m,-2)L(m,10)}}
for(n= -2;n<11;n++){if(n!=0){M(-2,n)L(11,n)}} fprintf(o,".006 W 0 0 0 RGB S\n");
M(-2,0)L(11.1,0) M(0, -2)L(0,10.1) fprintf(o,".03 W 0 0 0 RGB S\n");
M(0,M_E)L(11.,M_E) fprintf(o,".006 W 0 0 0 RGB S\n");
fprintf(o,"1 setlinejoin 1 setlinecap\n");
//DO(m,400){x=-1.99+.008*m;y=Re(FIT1(log(10.),x)); if(y>10.3) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".03 W .5 0 .5 RGB S\n");
DO(m,400){x=-1.80+.02*m; y=Re(FIT1(log(2.),x)); if(y>10.3) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".03 W .5 0 .5 RGB S\n");
DO(m,400){x=-1.78+.02*m; y=Re(FIT1(log(1.9),x)); if(y>10.33) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.77+.02*m; y=Re(FIT1(log(1.8),x)); if(y>10.33) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.74+.02*m; y=Re(FIT1(log(1.7),x)); if(y>10.3) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.72+.03*m; y=Re(FIT1(log(1.6),x)); if(y>10.3) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.68+.04*m; y=Re(FIT1(log(1.5),x)); if(y>10.3) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.64+.04*m; y=Re(FIT1(log(1.4),x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.58+.04*m; y=Re(FIT1(log(1.3),x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.52+.04*m; y=Re(FIT1(log(1.2),x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
DO(m,400){x=-1.42+.04*m; y=Re(FIT1(log(1.1),x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 .5 0 RGB S\n");
//
DO(m,400){x=-1.65+.04*m; y=Re(FIT1(1./M_E,x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".03 W 0 0 .7 RGB S\n");
DO(m,400){x=-1.64+.04*m; y=Re(FIT1(log(sqrt(2.)),x)); if(x>11.1) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".03 W .8 0 0 RGB S\n");
DO(m,400){x=-1.873+.01*m; y=Re(FIT1(1.,x)); if(y>11) break; if(m==0)M(x,y)else L(x,y)} fprintf(o,".02 W 0 0 0 RGB S\n");
//
DO(m,401){x=-1.987+.01*m;y=Re(F4TEN(x));if(m==0)M(x,y)else L(x,y)}fprintf(o,".02 W .5 0 .5 RGB S\n");
//DO(m,39){x=-1.89+.098*m; y=Re(FSEXP(x));if(m==0)M(x,y)else L(x,y)} fprintf(o,".01 W 0 0 0 RGB S\n");
//DO(m,48){x=-1.82+.0973*m; y=Re(F2(x)); if(m==0)M(x,y)else L(x,y)} fprintf(o,".01 W 0 0 0 RGB S\n");
//DO(m,126){x=-1.72+.1*m; y=Re(F15(x)); if(m==0)M(x,y)else L(x,y)} fprintf(o,".01 W 0 0 0 RGB S\n");
//DO(m,130){x=-1.65+.1*m; y=Re(E1ETF(x)); if(m==0)M(x,y)else L(x,y)} fprintf(o,".01 W 0 0 0 RGB S\n");
//DO(m,130){x=-1.64+.1*m; y=Re(f21E(x)); if(m==0)M(x,y)else L(x,y)} fprintf(o,".01 W 0 0 0 RGB S\n");
M(-1.998,-2)L(-1.992,-.01)L(-1,0)L(-.01,.01)L(0,1)L(.012,10.1) fprintf(o,".02 W .4 0 .8 RGB S\n");
M(-1.01,-2)L(-1,0)L(-.99,.99)L(0,1)L(10.97,1.01) fprintf(o,".02 W .2 .4 0 RGB S\n");
fprintf(o,"showpage\n%cTrailer",'%'); fclose(o);
       system("epstopdf tet10bx.eps");
       system( "open tet10bx.pdf"); //mac
// system( "xpdf tet10bx.pdf"); // linux
getchar(); system("killall Preview");// mac
}