real procedure nthroot(x, n);
  value x, n; real x, n;
begin
  nthroot := x ** (1/n);
end;