xbfplt.h

c源码

/* f.h                    freeware                          [e-mail] */


/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void gplt_f(
double xmin,
double xmax,
double ymin,
double ymax,
  char  feq[TXT_EQ]
)
{
FILE   *fp;

        fp = fopen("f_Df.plt","w");

fprintf(fp,"#\n");
fprintf(fp,"# Worksheet for gp371w32.exe by xhunga bernard 06000 Nice France\n");
fprintf(fp,"#\n\n");
fprintf(fp," set zeroaxis lt 8\n");
fprintf(fp," set grid\n\n");
fprintf(fp," plot [%0.3f:%0.3f] [%0.3f:%0.3f] ",xmin,xmax,ymin,ymax);
fprintf(fp,"%s",feq);

fprintf(fp,"\n\n");
fprintf(fp," reset");

 fclose(fp);
}


/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void gplt_root
(
double xmin,
double xmax,
double ymin,
double ymax,
  char  eqfunct[TXT_EQ],
double (*P_f)(double x),
double  root
)
{
double r;
FILE   *fp;
FILE   *fq;

        fp = fopen("f_Df.plt","w");
        fq = fopen(     "x_0","w");

fprintf(fp,"#\n");
fprintf(fp,"# Worksheet for gp371w32.exe by xhunga bernard 06000 Nice France\n");
fprintf(fp,"#\n\n");
fprintf(fp," set zeroaxis lt 8\n");
fprintf(fp," set grid\n\n");
fprintf(fp," plot [%0.3f:%0.3f] [%0.3f:%0.3f] ",xmin,xmax,ymin,ymax);
fprintf(fp,"%s",eqfunct);
fprintf(fp,", ");
fprintf(fp," \"x_0\" lt 3 ");

fprintf(fp,"\n\n");
fprintf(fp," reset");
    r =    (*P_f)(root);
fprintf(fq," %0.6f   %0.6f", root, r);

 fclose(fp);
 fclose(fq);
}