real procedure nthroot(x, n);
  value x, n; real x, n;
begin
   nthroot := exp(ln(x) / n);
end;