/* f.h                    freeware                          [[Email Removed]] */

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void gplt_gh(
double xmin,
double xmax,
double ymin,
double ymax,
  char  geq[TXT_EQ],
  char  heq[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",geq);
fprintf(fp,", " );
fprintf(fp,"%s",heq);

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

 fclose(fp);
}

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void gplt_gxh(
double xmin,
double xmax,
double ymin,
double ymax,
  char  geq[TXT_EQ],
  char  heq[TXT_EQ],
double (*P_g)(double x),
double  x_0
)
{
double t;

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",geq);
fprintf(fp,", " );
fprintf(fp,"%s",heq);
fprintf(fp,", ");
fprintf(fp," \"x_0\" lt 3 ");

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

t = (*P_g)(x_0);
fprintf(fq," %0.6f   %0.6f", x_0, t);

 fclose(fp);
 fclose(fq);
}