Mathc initiation/Fichiers h : c47

Sommaire

Installer ce fichier dans votre répertoire de travail.

	kg_3d.h

/* --------------------------------- */
/* save as kg_3d.h                   */
/* --------------------------------- */
void G_3d(
CTRL_splot p,
double (*P_f)(double x, double y),
point2d step
)
{
FILE   *fp;
double   x;
double   y;

        fp = fopen("a_main.plt","w");
fprintf(fp,"reset\n");
fprintf(fp,"set xlabel   \"<x>\"\n");
fprintf(fp,"set ylabel   \"<y>\"\n");
fprintf(fp,"set zlabel   \"<z>\"\n");

fprintf(fp,"set    samples 40\n");
fprintf(fp,"set isosamples 40\n");
fprintf(fp,"set view %0.3f, %0.3f, %0.3f, %0.3f \n",
                    p.rot_x,p.rot_z,p.scale,p.scale_z);
fprintf(fp,"set xrange [%0.3f:%0.3f]\n",p.xmin,p.xmax);
fprintf(fp,"set yrange [%0.3f:%0.3f]\n",p.ymin,p.ymax);
fprintf(fp,"set zrange [%0.3f:%0.3f]\n",p.zmin,p.zmax);
fprintf(fp,"splot ");
fprintf(fp," \"a_ka\" with lines,\\\n");
fprintf(fp," \"a_kb\" with lines");

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

  fclose(fp);

  fp = fopen("a_ka","w");
   for(x = p.xmin; x <= p.xmax; x += step.x)
   {
    fprintf(fp,"\n\n");
    for(y = p.ymin; y <= p.ymax; y += step.y)
        fprintf(fp," %6.3f   %6.3f   %6.3f\n", x, y, ((*P_f)(x,y)) );
   }
  fclose(fp);


  fp = fopen("a_kb","w");
   for(y = p.ymin; y <= p.ymax; y += step.y)
   {
    fprintf(fp,"\n\n");
    for(x = p.xmin; x <= p.xmax; x += step.x)
        fprintf(fp," %6.3f   %6.3f   %6.3f\n", x, y, ((*P_f)(x,y)) );
   }
  fclose(fp);
}
/* --------------------------------- */
/* --------------------------------- */

Déclaration des fichiers h.

---