x^y = e^(y*ln(x))
log_base_x(y)=ln(y)/ln(x)

bc -l:
define p(x,y) {return e(y*l(x));}
define r(x,y) {return l(y)/l(x);}