home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS 1992 September / Simtel20_Sept92.cdr / msdos / spredsht / prof2d.arc / EXAMPLE2.PR2 < prev    next >
Text File  |  1989-02-26  |  2KB  |  78 lines
  1. : Prof2d, The 2D Scientific Spreadsheet 
  2. : by Gertjan L. Ouwerling, Electrical Materials Lab, Delft University
  3. : -------------------------------------------------------------------
  4. : Example 2: Making contour plots from internally generated data.
  5. :
  6. : If you have a HP7550A or HP7475 plotter connected to your PC, 
  7. : you can try the contour plots by changing cview commands in cview.
  8. :
  9. : Execute this example by typing:  prof2d example2.pr2 <enter>
  10. : 
  11. : NOTE: some steps may be slow without a co-processor!
  12. :
  13. : -------------------------------------------------------------------
  14. : In this example, synthetic data is generated internally
  15. : and two contour plots of this data are drawn.
  16. : These plots are depicted following this example.
  17. : The data is also output to file in several formats.
  18. :
  19. : Declaration of program data types:
  20. :
  21. row x1 x2 $
  22. column y1 y2 $
  23. field z1 z2 $
  24. :
  25. : Data space is allocated and x- and y-axes are made;
  26. : x and y both run between 0 and 2*pi.
  27. :
  28. nx=20
  29. ny=30
  30. x1=xcount/(nx-1)*2*pi
  31. y1=ycount/(ny-1)*2*pi
  32. :
  33. : The value of the field z1 is computed using x and y
  34. :
  35. z1=sin(x1)+cos(y1)
  36. :
  37. : A contour plot of z1 is made in one color, without markers,
  38. : and on A4 size paper. The contour values are chosen by setcont.
  39. : xscl and yscl are not altered, hence autoscaling with identical
  40. : scales for x and y will be used.
  41. :
  42. : ENABLE PLOTTING BY CHANGING COMMAND cview IN cplot !
  43. :
  44. asize=4
  45. color=0
  46. marker=0
  47. setcont -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 $
  48. cview x1 y1 z1 sin(x)+cos(y) | Demonstration plot
  49. :
  50. : A 2nd contour plot is made with distorted axes.
  51. : The other condition are equal, except that now markers
  52. : are drawn.
  53. :
  54. marker=1
  55. x2=sqrt(x1)/sqrt(2*pi)
  56. y2=sqrt(y1)/sqrt(2*pi)
  57. cview x2 y2 z1 sin(x1)+cos(y1) <-> x2,y2 | x2=sqrt(x1)/sqrt(2*pi) | Demo plot 2
  58. :
  59. : Some data is saved in listed form.
  60. : Before doing so, the fields are thinned to save disk space.
  61. :
  62. redx 2
  63. redy 2
  64. put list.dat x1 x2 y1 y2 z1 $
  65. :
  66. : And some data is saved in matrix form.
  67. :
  68. mput mat.dat x1 y2 z1
  69. :
  70. : End of this prof2d program
  71. :
  72. : -----------------------------------
  73. :
  74. : Give  quit <enter> to leave Prof2d
  75. : or continue playing around a little
  76. :
  77.  
  78.