conj.cpp

   1 // Complex conjugate
   2
   3 #include "stdafx.h"
   4 #include "defs.h"
   5
   6 void
   7 eval_conj(void)
   8 {
   9         push(cadr(p1));
  10         eval();
  11         p1 = pop();
  12         push(p1);
  13         if (!find(p1, imaginaryunit)) { // example: (-1)^(1/3)
  14                 polar();
  15                 conjugate();
  16                 clockform();
  17         } else
  18                 conjugate();
  19 }
  20
  21
  22 void
  23 conjugate(void)
  24 {
  25         push(imaginaryunit);
  26         push(imaginaryunit);
  27         negate();
  28         subst();
  29         eval();
  30 }