home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / mbug / mbug003.arc / JRTPAS-1.LBR / JRT38.PAS < prev    next >
Pascal/Delphi Source File  |  1979-12-31  |  768b  |  29 lines
  1. {Program 6.2;
  2. extend program 4.9 to print x-axis}
  3.  
  4. Program graph2;        
  5.  
  6. const    d = 0.0625;        {1/16, 16 lines for interval [x, x + 1]}
  7.     s = 32;            {32 character widths for interval [y, y + 1]}
  8.     h1 = 34;            {character position of x-axis}
  9.     h2 = 68;            {line width}
  10.     c = 6.28318;              {2 * pi}  
  11.     lim = 32;
  12. var    i, j, k, n : integer;
  13.     x, y : real;
  14.     a : array[ 1..h2 ] of char;
  15. function exp ( x : real ): real; extern;
  16. function sin ( x : real ): real; extern;
  17.  
  18. begin
  19.     for j := 1 to h2 do a[ j ] := ' ';
  20.     for i := 0 to lim do
  21.     begin
  22.         x := d * i; y := exp( -x ) * sin( c * x );
  23.         a[ h1 ] := ':'; n := round( s*y ) + h1; a[ n ] := '*';
  24.         if n < h1 then k := h1 else k := n;
  25.         for j := 1 to k do write ( a[ j ] );
  26.         writeln; a[ n ] := ' '
  27.     end
  28. end.
  29.