function r = m_nthroot(n, A)
  r = A / n;
  m = n - 1;
  do
    d = (A ./ r .^ m - r) / n;
    r+= d;
  until (abs(d) < abs(r * 1e-9))
endfunction