home *** CD-ROM | disk | FTP | other *** search
/ Club Amiga de Montreal - CAM / CAM_CD_1.iso / files / 344b.lha / plplot_v2.6 / examples / example08.c < prev    next >
Encoding:
C/C++ Source or Header  |  1990-01-27  |  1.3 KB  |  62 lines
  1. /* Demonstration of 3-d plotting */
  2.  
  3. #include "plplot.h"
  4. #include <math.h>
  5.  
  6. #define   NPTS   41
  7.  
  8. static int opt[] = {1,   2,   3,   3};
  9.  
  10. static float alt[] = {60.0,20.0,60.0,60.0};
  11.  
  12. static float az[] = {30.0,60.0,120.0,160.0};
  13.  
  14. static char *title[4] =  {
  15.    "#frPLPLOT Example 8 - Alt=60, Az=30, Opt=1",
  16.    "#frPLPLOT Example 8 - Alt=20, Az=60, Opt=2",
  17.    "#frPLPLOT Example 8 - Alt=60, Az=120, Opt=3",
  18.    "#frPLPLOT Example 8 - Alt=60, Az=160, Opt=3"
  19. };
  20.  
  21. static float x[NPTS], y[NPTS], z[NPTS][NPTS];
  22.  
  23. main()
  24. {
  25.    int i, j, k;
  26.  
  27.    float xx, yy, r;
  28.  
  29.    for (i=0; i<NPTS; i++)  {
  30.       x[i] = (i-(NPTS/2))/(double)(NPTS/2);
  31.       y[i] = (i-(NPTS/2))/(double)(NPTS/2);
  32.    }
  33.  
  34.    for (i=0; i<NPTS; i++)  {
  35.       xx = x[i];
  36.       for (j=0; j<NPTS; j++)  {
  37.          yy = y[j];
  38.          r = sqrt(xx*xx + yy*yy);
  39.          z[i][j] = exp(-r*r) * cos(2.0*3.141592654*r);
  40.       }
  41.    }
  42.  
  43.    plstar(1,1);
  44.  
  45.    for (k=0; k<4; k++)  {
  46.       pladv(0);
  47.       plvpor(0.0,1.0,0.0,0.9);
  48.       plwind(-1.0,1.0,-0.9,1.1);
  49.       plcol(1);
  50.       plw3d(1.0,1.0,1.0,-1.0,1.0,-1.0,1.0,-1.0,1.0,alt[k],az[k]);
  51.       plbox3("bnstu","x axis",0.0,0,"bnstu","y axis",0.0,0,
  52.              "bcdmnstuv","z axis",0.0,0);
  53.       plcol(2);
  54.       plot3d(x,y,&z[0][0],NPTS,NPTS,NPTS,opt[k],1);
  55.       plcol(3);
  56.       plmtex("t",1.0,0.5,0.5,title[k]);
  57.    }
  58.  
  59.    pltext();
  60.    plend();
  61. }
  62.