*/
If you have a PH account, you can customize your PH profile.
Search
*/

View \xisrevfg.h

Calculus (Derivative:Tangent). v04

Submitted By: xhunga
Rating: (Not rated) (Rate It)

File Details
ZIP File Browser
/* f.h                    freeware                          [[Email Removed]] */

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void gplt_fgab(
double xmin,
double xmax,
double ymin,
double ymax,
double a,
double b,
  char feq[TXT_EQ],
  char geq[TXT_EQ],
double (*P_f)(double x),
double (*P_g)(double x)
)
{
FILE   *fp;
char    FA[FILENAME] = "a";
char    FB[FILENAME] = "b";
double  t;

        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,", ");
fprintf(fp,"%s",geq);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 1",FA);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 1",FB);


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

             fclose(fp);

        fp = fopen(FA,  "w");

t = ((*P_f)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a, t);
t = ((*P_g)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a, t);

             fclose(fp);

        fp = fopen(FB,  "w");

t = ((*P_f)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b, t);
t = ((*P_g)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b, t);

             fclose(fp);
}

/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void   gplt_SolidRevolDiskfg(
double xmin,
double xmax,
double ymin,
double ymax,
double a,
double b,
double r,
double step,
double (*P_f)(double x),
double (*P_g)(double x)
)
{
FILE   *fp;

char    FA[FILENAME] = "a";
char    FB[FILENAME] = "b";
char    FD[FILENAME] = "d";
char    FE[FILENAME] = "e";
char    FR[FILENAME] = "r";
char    FL[FILENAME] = "l";
char    FM[FILENAME] = "m";
char    FO[FILENAME] = "o";
char    FP[FILENAME] = "p";

double  t;
double  i;

        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");
t = ((*P_f)((a+b)/2.)) + 5. * r;
fprintf(fp," set label \"dx\" at %0.3f, %0.3f\n\n", ((a+b)/2.), t);
t = ((*P_f)((a+b)/2.));
fprintf(fp," set label \"f : outer radius\" at %0.3f, %0.3f\n\n",
             ((a+b)/2.)+2*r, t);
t = ((*P_g)((a+b)/2.));
fprintf(fp," set label \"g : inner radius\" at %0.3f, %0.3f\n\n",
             ((a+b)/2.)+2*r, t);
fprintf(fp," plot [%0.3f:%0.3f] [%0.3f:%0.3f] ",xmin,xmax,ymin,ymax);
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FA);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FB);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FD);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FE);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FL);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FM);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FO);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FP);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 9 pt 0",FR);
fprintf(fp,"\n\n");
fprintf(fp," reset");
             fclose(fp);


        fp = fopen(FA,  "w");
t = ((*P_f)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a,  t);
t = ((*P_g)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a,  t);
             fclose(fp);


        fp = fopen(FB,  "w");
t = ((*P_f)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  t);
t = ((*P_g)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  t);
             fclose(fp);


        fp = fopen(FD,  "w");
t = ((*P_f)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a, -t);
t = ((*P_g)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a, -t);
             fclose(fp);


        fp = fopen(FE,  "w");
t = ((*P_f)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  -t);
t = ((*P_g)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  -t);
             fclose(fp);


        fp = fopen(FR,  "w");
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.-(r), t);
t = (*P_g)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.-(r), t);
t = (*P_g)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.+(r), t);
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.+(r), t);
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.-(r), t);
             fclose(fp);


        fp = fopen(FL,  "w");
i = a;
do{
   t = ((*P_f)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FM,  "w");
i = a;
do{
   t = ((*P_g)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FO,  "w");
i = a;
do{
   t = ((*P_f)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, -t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FP,  "w");
i = a;
do{
   t = ((*P_g)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, -t);
   i += step;
}while(i <=  b);
             fclose(fp);
}


/* ------------------------------------- FUNCTION --------------------------  */
/* Do    :                                                                    */
/*                                                                            */
/* -------------------------------------------------------------------------- */
void   gplt_SolidRevolCylindShellfg(
double xmin,
double xmax,
double ymin,
double ymax,
double a,
double b,
double r,
double step,
double (*P_f)(double x),
double (*P_g)(double x)
)
{
FILE   *fp;
char    FA[FILENAME] = "a";
char    FB[FILENAME] = "b";
char    FD[FILENAME] = "d";
char    FE[FILENAME] = "e";
char    FL[FILENAME] = "l";
char    FM[FILENAME] = "m";
char    FO[FILENAME] = "o";
char    FP[FILENAME] = "p";
char    FR[FILENAME] = "r";
char    FX[FILENAME] = "x";

double  t;
double  i;

        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");
t = ((*P_f)((a+b)/2.)) + 2. * r;
fprintf(fp," set label \"dx\" at %0.3f, %0.3f\n\n", ((a+b)/2.), t);
t = ((*P_g)((a+b)/2.)) - 2. * r;
fprintf(fp," set label \"average radius\" at %0.3f, %0.3f \n\n", ((a+b)/2.), t);
t = (((*P_f)((a+b)/2.)) + ((*P_g)((a+b)/2.))) / 2.;
fprintf(fp," set label \"altitude : f-g\" at %0.3f, %0.3f \n\n", ((a+b)/2.)+2*r, t);
fprintf(fp," plot [%0.3f:%0.3f] [%0.3f:%0.3f] ",xmin,xmax,ymin,ymax);
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FA);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FB);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FD);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FE);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FL);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FM);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FO);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 0",FP);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 9 pt 0",FR);
fprintf(fp,", ");
fprintf(fp," \"%s\" with linesp lt 3 pt 1",FX);
fprintf(fp,"\n\n");
fprintf(fp," reset");
             fclose(fp);


        fp = fopen(FA,  "w");
t = ((*P_f)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a,  t);
t = ((*P_g)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", a,  t);
             fclose(fp);


        fp = fopen(FB,  "w");
t = ((*P_f)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  t);
t = ((*P_g)(b)) ;
fprintf(fp," %0.6f   %0.6f\n", b,  t);
             fclose(fp);


        fp = fopen(FD,  "w");
t = ((*P_f)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", -a,  t);
t = ((*P_g)(a)) ;
fprintf(fp," %0.6f   %0.6f\n", -a,  t);
             fclose(fp);


        fp = fopen(FE,  "w");
t = ((*P_f)(b)) ;
fprintf(fp," %0.6f   %0.6f\n",-b,  t);
t = ((*P_g)(b)) ;
fprintf(fp," %0.6f   %0.6f\n",-b,  t);
             fclose(fp);


        fp = fopen(FR,  "w");
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.- (r ),   t);
t = (*P_g)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.- (r ),   t);
t = (*P_g)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.+ (r ),   t);
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.+ (r ),   t);
t = (*P_f)( (a+b)/2.) ;
fprintf(fp," %0.6f   %0.6f\n", (a+b)/2.- (r ),   t);
             fclose(fp);


        fp = fopen(FL,  "w");
i = a;
do{
   t = ((*P_f)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FM,  "w");
i = a;
do{
   t = ((*P_g)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", i, t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FO,  "w");
i = a;
do{
   t = ((*P_f)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", -i, t);
   i += step;
}while(i <=  b);
             fclose(fp);


        fp = fopen(FP,  "w");
i = a;
do{
   t = ((*P_g)(i)) ;
   fprintf(fp," %0.6f   %0.6f\n", -i, t);
   i += step;
}while(i <=  b);
             fclose(fp);

         fp = fopen(FX,  "w");
t = ((*P_g)((a+b)/2.)) ;
fprintf(fp," %0.6f   %0.6f\n", ((a+b)/2.),  t);
             fclose(fp);
}

corner
Home About PH Privacy Statement Contact Us Advertise Here
© 1996-2008 CommunityHeaven LLC. All rights reserved. Reproduction in whole or in part, in any form or medium without express written permission is prohibited.
Violators of this policy may be subject to legal action. Please read our Terms Of Use and Privacy Statement for more information.
North American business development: Nicolai Wadstrom. Publisher: Lars Hagelin.