home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / simtel / cpmug / cpmug083.ark / TRGNEG%.ALG < prev    next >
Encoding:
Text File  |  1984-04-29  |  2.1 KB  |  79 lines
  1.  file TRGNEG%.ALG                            GAE - Feb. 1982
  2.  
  3. Make a new version of file TRGNEG.ALG as follows:
  4. (1)  Begin with the original file TRGNEG.ALG.
  5. (2)  Delete everything from  FOURONPI: 4 / #PI $
  6.           up to (but not including) PROPERTY SIN, *,
  7.           (more than 50 lines).
  8. (3)  Replace the deleted material with this file starting
  9.           below this line:
  10. -------------------------------------------------------------
  11.  
  12. FUNCTION SIN (EX1,
  13.   % Local: %  EX2),
  14.   WHEN NEGLT (EX1), -SIN(-EX1) EXIT,
  15.   WHEN NUMBER(EX2:EX1/#PI),
  16.     COSPI(HALF - EX2)  EXIT,
  17.   WHEN POSMULT (TRGEXPD, 7),
  18.     #I * (1/(EX2:#E^(#I*EX1)) - EX2) * HALF  EXIT,
  19.   SIMPU ('SIN, EX1),
  20. ENDFUN $
  21.  
  22. FUNCTION COS (EX1,
  23.   % Local: %  EX2),
  24.   WHEN NEGLT (EX1), COS(-EX1) EXIT,
  25.   WHEN NUMBER(EX2:EX1/#PI),
  26.     COSPI(EX2)  EXIT,
  27.   WHEN POSMULT (TRGEXPD, 7),
  28.     (1/(EX2:#E^(#I*EX1)) + EX2) * HALF  EXIT,
  29.   SIMPU ('COS, EX1),
  30. ENDFUN $
  31.  
  32.  
  33. ROOT5: 5^HALF $
  34. SIXRT5: 6*ROOT5 $
  35.  
  36. COSPION5:  ( 1+ROOT5)/4 $
  37. COS2PION5: (-1+ROOT5)/4 $
  38.  
  39. COSPION15:  (-1+ROOT5+(30+SIXRT5)^HALF)/8 $
  40. COS2PION15: ( 1+ROOT5+(30-SIXRT5)^HALF)/8 $
  41. COS4PION15: ( 1-ROOT5+(30+SIXRT5)^HALF)/8 $
  42. COS7PION15: (-1-ROOT5+(30-SIXRT5)^HALF)/8 $
  43.  
  44.  
  45. FUNCTION COSPI(EX1,
  46.   % local : % EX2),
  47.   WHEN EX1 < 0,  COSPI(-EX1)  EXIT,
  48.   BLOCK
  49.     WHEN EX1 > 2,
  50.       EX1: EX1 - 2*QUOTIENT(NUM(EX1),2*DEN(EX1))  EXIT,
  51.   ENDBLOCK,
  52.   BLOCK
  53.     WHEN EX1 > 1, EX1: 2-EX1  EXIT,
  54.   ENDBLOCK,
  55.   WHEN MULTIPLE(DEN(EX1),2),
  56.     WHEN EX1 > HALF,  -((1+COSPI(2*EX1))*HALF)^HALF  EXIT,
  57.     ((1+COSPI(2*EX1))*HALF)^HALF  EXIT,
  58.   EX1: DIVIDE(NUM(EX1), EX2: DEN(EX1)),
  59.   WHEN MOD(POP(EX1),2) = 1, -COSPI1(EX1,EX2)  EXIT,
  60.   COSPI1(EX1,EX2),
  61. ENDFUN $
  62.  
  63. FUNCTION COSPI1(EX1,EX2),
  64.   WHEN EX1 EQ 0, 1  EXIT,
  65.   WHEN LESSER (EX2, 2*EX1), -COSPI1(EX2-EX1,EX2)  EXIT,
  66.   WHEN EX2 EQ 3, HALF  EXIT,
  67.   WHEN EX2 EQ 5,
  68.     WHEN EX1 EQ 1, COSPION5  EXIT,
  69.     COS2PION5  EXIT,
  70.   WHEN EX2 EQ 15,
  71.     WHEN EX1 EQ 1, COSPION15  EXIT,
  72.     WHEN EX1 EQ 2, COS2PION15  EXIT,
  73.     WHEN EX1 EQ 4, COS4PION15  EXIT,
  74.     COS7PION15  EXIT,
  75.   LIST('COS,#PI*EX1/EX2),
  76. ENDFUN $
  77.  
  78.  
  79.