home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / CONTOUR.C < prev    next >
Encoding:
C/C++ Source or Header  |  1988-12-29  |  1.9 KB  |  87 lines
  1. #include <math.h>
  2. #include <graph.h>
  3. #include "worlddr.h"
  4. #include "segraph.h"
  5.  
  6. int  err,i,j;
  7. float ContourMap[61][61];
  8.  
  9. float Calc3DFunction(float x,float y)
  10. {
  11.   float Calc3DFunctionResult;
  12.   Calc3DFunctionResult = 3.0 * sin(sqrt(x*x+y*y));
  13.   return(Calc3DFunctionResult );
  14. }
  15.  
  16. float Calc3DFunction2(float x,float y)
  17. {
  18.   float Calc3DFunctionResult;
  19.  
  20.   Calc3DFunctionResult = 7.0 * exp(-0.1 *(x*x+y*y));
  21.   return(Calc3DFunctionResult );
  22.  
  23. }
  24.  
  25. void  CreateFunction(){
  26. int i,j;
  27. float  x,y,z ;
  28.  
  29.   for ( i = 0; i <= 60; i++ )  {
  30.     x = (i-30)/4.0;
  31.     for ( j = 0; j <= 60; j++ ){
  32.       y = (j-30)/4.0;
  33.       z = Calc3DFunction(x,y);
  34.       ContourMap[i][j] = z;
  35.     }
  36.   }
  37. }
  38.  
  39.  
  40. void main(){
  41.  
  42. int ContourColors[maxContours];
  43. int ContourHatch[maxContours];
  44.   for ( i = 0; i <= maxContours - 1 ; i++ ){
  45.     ContourColors[i] = (i%15)+ 1;
  46.     ContourHatch[i] = 1;
  47.   }
  48.   CreateFunction();
  49.   InitSEGraphics(6);
  50.  
  51.   for (i = 2; i <= 2; i++){
  52.     SetPercentWindow( 0.05,0.05, 0.95,0.80,4);
  53.     SetPercentWindow( 0.05,0.83,0.95,0.95,3);
  54.  
  55.     SetCurrentWindow(4);
  56.     SetWin2PlotRatio(4, 0.15, 0.15, 0.1, 0.15 );
  57.     BorderCurrentWindow(2);
  58.     SetAxesType(0,0);
  59.     ScalePlotArea( -20.0, -20.0, 20.0, 20.0);
  60.     SelectColor(15);
  61.     SetXYIntercepts(-20.0,-20.0);
  62.     DrawYAxis(1.0,0);
  63.     LabelYAxis(10, 0);
  64.     DrawXAxis(1.0,0);
  65.     LabelXAxis(10,0);
  66.     
  67.     SetXYIntercepts(20.0,20.0);
  68.     DrawYAxis(1.0,1);
  69.     LabelYAxis(10,1);
  70.     DrawXAxis(1.0,1);
  71.     LabelXAxis(10,1);
  72.     TitleWindow("Contour Plot z = Sin(Sqrt(x*x + y*y))");
  73.     TitleXAxis("X Axis Label");
  74.     TitleYAxis("Y Axis Label");
  75.     setlinestyleXX(1,1,1);
  76.     DrawGridX(10);
  77.     DrawGridY(10);
  78.     ContourPlot(&ContourMap[0][0],61,61,1.0, ContourColors,ContourHatch);
  79.     SetCurrentWindow(3);
  80.     BorderCurrentWindow(5);
  81.     ContourPlotLegends(&ContourMap[0][0],61,61,1.0, ContourColors,ContourHatch);
  82.  
  83.   }
  84.   getch();
  85.   CloseSEGraphics();
  86. }
  87.