eigenmath/polar.cpp

44 lines
580 B
C++

/* polar(z) = mag(z) * (-1) ^ (arg(z) / pi)
For example polar(exp(i pi/3)) gives the result (-1)^(1/3)
*/
#include "stdafx.h"
#include "defs.h"
void
eval_polar(void)
{
push(cadr(p1));
eval();
p1 = pop();
push(p1);
mag();
push_integer(-1);
push(p1);
arg();
push(symbol(PI));
divide();
power();
multiply();
}
static char *s[] = {
"polar(exp(i pi/3))",
"(-1)^(1/3)",
"polar(exp(-i pi/3))",
"-(-1)^(2/3)",
"rect(polar(3+4*i))", // needs sin(arctan(x)) and cos(arctan(x))
"3+4*i",
};
void
test_polar(void)
{
test(__FILE__, s, sizeof s / sizeof (char *));
}