Mathc initiation/Fichiers h : 812
Installer ce fichier dans votre répertoire de travail.
kg_3d.h |
---|
/* --------------------------------- */
/* save as kg_3d.h */
/* --------------------------------- */
void G_3d_eq(
CTRL_splot p,
char feq[],
double (*P_f)(double x, double y),
pt2d Q
)
{
FILE *fp;
fp = fopen("a_main.plt","w");
fprintf(fp,"reset\n");
fprintf(fp,"set samples 100\n");
fprintf(fp,"set isosamples 100\n");
fprintf(fp,"set hidden3d\n");
fprintf(fp,"set dgrid3d\n");
fprintf(fp,"set contour both\n");
fprintf(fp,"set cntrparam levels disc %6.3f\n",(*P_f)(Q.x,Q.y));
fprintf(fp,"set xlabel \"X axis\"\n");
fprintf(fp,"set ylabel \"Y axis\"\n");
fprintf(fp,"set zlabel \"Z axis\" offset 1, 0\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\" lt 7 pt 7 ps 3,\\\n");
fprintf(fp,"%s lt 9",feq);
fprintf(fp,"\n\n");
fclose(fp);
fp = fopen("a_ka","w");
fprintf(fp," %6.3f %6.3f %6.3f\n",
Q.x, Q.y, ((*P_f)(Q.x,Q.y)));
fclose(fp);
}
/* --------------------------------- */
/* --------------------------------- */