newtonsolveΒΆ

newtonsolve(mu, xold, sigoldsq, N, Nmax)