#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#define NDMAX 1000
void main_part(char fnameR[],char fnameW[]);
void TeX_DRN_XJIKU(FILE *fout,double xmin,double xmax,double dx,
    double kxi,double kxf,double kyi,double kyf);
void TeX_DRN_YJIKU(FILE *fout,double ymin,double ymax,double dy,
    double kxi,double kxf,double kyi,double kyf);
void TeX_DRL_XJIKU(FILE *fout,double xmin,double xmax,
    double kxi,double kxf,double kyi,double kyf);
void TeX_DRL_YJIKU(FILE *fout,double ymin,double ymax,
    double kxi,double kxf,double kyi,double kyf);
/*---------------------------------------------------------------------------------*/

int main(int argc,char *argv[])
{
    char fnameR[50]; /*入力ファイル名*/
    char fnameW[50]; /*出力ファイル名*/

    strcpy(fnameR,argv[1]);
    strcpy(fnameW,argv[2]);

    main_part(fnameR,fnameW);

    printf("リターンキーを押してください");
    getchar();
    return 0;
}
/*---------------------------------------------------------------------------------*/

void main_part(char fnameR[],char fnameW[])
{
    /*作図用変数宣言*/
    double HSIZE;    /*画像横寸法（mm）*/
    double DHL;      /*画像左余白寸法：縦軸描画用（mm）*/
    double DHR;      /*画像右余白寸法（mm）*/
    double VSIZE;    /*画像縦寸法（mm）*/
    double DVL;      /*画像下余白寸法：横軸描画用（mm）*/
    double DVU;      /*画像上余白寸法（mm）*/
    char sxjiku[50]; /*ｘ軸名（単位含む）*/
    double xmin;     /*ｘ軸最小値*/
    double xmax;     /*ｘ軸最大値*/
    double dx;       /*ｘ軸増分*/
    char slogx[2];   /*ｘ軸の種類（L：対数，N：普通）*/
    char syjiku[50]; /*ｙ軸名（単位含む）*/
    double ymin;     /*ｙ軸最小値*/
    double ymax;     /*ｙ軸最大値*/
    double dy;       /*ｙ軸増分*/
    char slogy[2];   /*ｙ軸の種類（L：対数，N：普通）*/
    double kxi,kxf,kyi,kyf;
    double xx,yy;
    double kxx,kyy;
    int nd; /*データ組数−１
    int i;
    /*プロットデータ格納配列宣言*/
    double datax[NDMAX];
    double datay[NDMAX];
    /*入出力用変数宣言*/
    FILE *fin,*fout;
    int i,j,n,m;
    char dat[256];
    char sv[10];

    /*データ読み込み*/
    fin=fopen(fnameR,"r");
    fgets(dat,sizeof dat,fin);
    fgets(dat,sizeof dat,fin);
        strcpy(sxjiku,strtok(dat,","));
        strcpy(syjiku,strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        HSIZE=atof(strtok(dat,","));
        DHL  =atof(strtok(NULL,","));
        DHR  =atof(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        VSIZE=atof(strtok(dat,","));
        DVL  =atof(strtok(NULL,","));
        DVU  =atof(strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        xmin=atof(strtok(dat,","));
        xmax=atof(strtok(NULL,","));
        dx  =atof(strtok(NULL,","));
        strcpy(slogx,strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        ymin=atof(strtok(dat,","));
        ymax=atof(strtok(NULL,","));
        dy  =atof(strtok(NULL,","));
        strcpy(slogy,strtok(NULL,","));
    fgets(dat,sizeof dat,fin);
        nd=atof(strtok(dat,","))-1;
    for(i=0;i<=nd;i++){
        fgets(dat,sizeof dat,fin);
            datax[i]=atof(strtok(dat,","));
            datay[i]=atof(strtok(NULL,","));
    }
    fclose(fin);

    /*データ対数値処理*/
    if(strncmp(slogx,"L",1)==0){
        for(i=0;i<=nd;i++){
            if(datax[i]<pow(10.0,xmin)){datax[i]=pow(10.0,xmin);}
            datax[i]=log10(datax[i]);
        }
    }
    if(strncmp(slogy,"L",1)==0){
        for(i=0;i<=nd;i++){
            if(datay[i]<pow(10.0,ymin)){datay[i]=pow(10.0,ymin);}
            datay[i]=log10(datay[i]);
        }
    }

    /*プロット領域定義*/
    kxi=DHL;
    kxf=HSIZE-DHR;
    kyi=DVL;
    kyf=VSIZE-DVU;

    /*TeXファイル書き込み*/
    fout=fopen(fnameW,"w");
    fprintf(fout,"\\documentclass[a4paper,10pt]{jsarticle}\n");
    fprintf(fout,"\\input{jartfmt}\n");
    fprintf(fout,"\\usepackage[dvipdfm]{graphics,graphicx}\n");
    fprintf(fout,"\\begin{document}\n");
    fprintf(fout,"\\begin{center}\n");
    fprintf(fout,"\\unitlength=1mm\n");
    fprintf(fout,"\\begin{picture}(%.3f,%.3f)\n",HSIZE,VSIZE);
    /*ｘ−ｙ軸描画*/
    if(strncmp(slogx,"N",1)==0){TeX_DRN_XJIKU(fout,xmin,xmax,dx,kxi,kxf,kyi,kyf);}
    if(strncmp(slogy,"N",1)==0){TeX_DRN_YJIKU(fout,ymin,ymax,dy,kxi,kxf,kyi,kyf);}
    if(strncmp(slogx,"L",1)==0){TeX_DRL_XJIKU(fout,xmin,xmax,kxi,kxf,kyi,kyf);}
    if(strncmp(slogy,"L",1)==0){TeX_DRL_YJIKU(fout,ymin,ymax,kxi,kxf,kyi,kyf);}
    /*枠線描画*/
    fprintf(fout,"\\path(%.3f,%.3f)(%.3f,%.3f)\n",kxi,kyi,kxf,kyi);
    fprintf(fout,"\\path(%.3f,%.3f)(%.3f,%.3f)\n",kxf,kyi,kxf,kyf);
    fprintf(fout,"\\path(%.3f,%.3f)(%.3f,%.3f)\n",kxf,kyf,kxi,kyf);
    fprintf(fout,"\\path(%.3f,%.3f)(%.3f,%.3f)\n",kxi,kyf,kxi,kyi);
    /*ｘ軸名描画*/
    fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){%s}}\n",0.5*(kxi+kxf),kyi-8.0,sxjiku);
    /*ｙ軸名描画*/
    fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){\\rotatebox{90}{%s}}}\n",kxi-15.0,0.5*(kyi+kyf),syjiku);
    /*データプロット*/
    fprintf(fout,"\\newcommand{\\plotchar}{\\makebox(0,0){$\\bullet$}}\n");
    fprintf(fout,"\\begin{drawjoin}\n");
    for(i=0;i<=nd;i++){
        xx=datax[i];
        yy=datay[i];
        kxx=kxi+(xx-xmin)*(kxf-kxi)/(xmax-xmin);
        kyy=kyi+(yy-ymin)*(kyf-kyi)/(ymax-ymin);
        fprintf(fout,"\\jput(%.3f,%.3f){\\plotchar}\n",kxx,kyy);
    }
    fprintf(fout,"\\end{drawjoin}\n");
    fprintf(fout,"\\end{picture}\n");
    fprintf(fout,"\\end{center}\n");
    fprintf(fout,"\\end{document}\n");
    fclose(fout);
}
/*---------------------------------------------------------------------------------*/
void TeX_DRN_XJIKU(FILE *fout,double xmin,double xmax,double dx,
    double kxi,double kxf,double kyi,double kyf)
{
    int i,ix;
    double xx,kxx,wv;
    char sv[10];
    /*普通ｘ軸描画*/
    ix=(int)((xmax-xmin)/dx);
    for(i=0;i<=ix;i++){
        xx=xmin+(double)i*dx;
        kxx=kxi+(xx-xmin)*(kxf-kxi)/(xmax-xmin);
        wv=xmin+(double)i*dx;
        sprintf(sv,"%.0f",wv);
        if((int)(dx*1000.0)%10!=0)sprintf(sv,"%.3f",wv);
        if((int)(dx*1000.0)%10==0&&(int)(dx*100.0)%10!=0)sprintf(sv,"%.2f",wv);
        if((int)(dx*1000.0)%10==0&&(int)(dx*100.0)%10==0&&(int)(dx*10.0)%10!=0)sprintf(sv,"%.1f",wv);
        fprintf(fout,"\\dottedline{1}(%.3f,%.3f)(%.3f,%.3f)\n",kxx,kyi,kxx,kyf);
        fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){%s}}\n",kxx,kyi-3.0,sv);
    }
}
/*---------------------------------------------------------------------------------*/
void TeX_DRN_YJIKU(FILE *fout,double ymin,double ymax,double dy,
    double kxi,double kxf,double kyi,double kyf)
{
    int i,iy;
    double yy,kyy,wv;
    char sv[10];
    /*普通ｙ軸描画*/
    iy=(int)((ymax-ymin)/dy);
    for(i=0;i<=iy;i++){
        yy=ymin+(double)i*dy;
        kyy=kyi+(yy-ymin)*(kyf-kyi)/(ymax-ymin);
        wv=ymin+(double)i*dy;
        sprintf(sv,"%.0f",wv);
        if((int)(dy*1000.0)%10!=0)sprintf(sv,"%.3f",wv);
        if((int)(dy*1000.0)%10==0&&(int)(dy*100.0)%10!=0)sprintf(sv,"%.2f",wv);
        if((int)(dy*1000.0)%10==0&&(int)(dy*100.0)%10==0&&(int)(dy*10.0)%10!=0)sprintf(sv,"%.1f",wv);
        fprintf(fout,"\\dottedline{1}(%.3f,%.3f)(%.3f,%.3f)\n",kxi,kyy,kxf,kyy);
        fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){%s}}\n",kxi-6.0,kyy,sv);
    }
}
/*---------------------------------------------------------------------------------*/
void TeX_DRL_XJIKU(FILE *fout,double xmin,double xmax,
    double kxi,double kxf,double kyi,double kyf)
{
    int i,j;
    double xx,kxx;
    /*対数ｘ軸描画*/
    for(i=(int)xmin;i<=(int)xmax-1;i++){
        for(j=1;j<=9;j++){
            xx=log10((double)j*pow(10.0,(double)i));
            kxx=kxi+(xx-xmin)*(kxf-kxi)/(xmax-xmin);
            fprintf(fout,"\\dottedline{1}(%.3f,%.3f)(%.3f,%.3f)\n",kxx,kyi,kxx,kyf);
        }
    }
    for(i=(int)xmin;i<=(int)xmax;i++){
        xx=(double)i;
        kxx=kxi+(xx-xmin)*(kxf-kxi)/(xmax-xmin);
        fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){10$^{%d}$}}\n",kxx,kyi-3.0,i);
    }
}
/*---------------------------------------------------------------------------------*/
void TeX_DRL_YJIKU(FILE *fout,double ymin,double ymax,
    double kxi,double kxf,double kyi,double kyf)
{
    int i,j;
    double yy,kyy;
    /*対数ｙ軸描画*/
    for(i=(int)ymin;i<=(int)ymax-1;i++){
        for(j=1;j<=9;j++){
            yy=log10((double)j*pow(10.0,(double)i));
            kyy=kyi+(yy-ymin)*(kyf-kyi)/(ymax-ymin);
            fprintf(fout,"\\dottedline{1}(%.3f,%.3f)(%.3f,%.3f)\n",kxi,kyy,kxf,kyy);
        }
    }
    for(i=(int)ymin;i<=(int)ymax;i++){
        yy=(double)i;
        kyy=kyi+(yy-ymin)*(kyf-kyi)/(ymax-ymin);
        fprintf(fout,"\\put(%.3f,%.3f){\\makebox(0,0){10$^{%d}$}}\n",kxi-5.0,kyy,i);
    }
}
/*---------------------------------------------------------------------------------*/