x_newtonsolve

x_newtonsolve(muone, xold, sig2old, cornum, z, alpha, beta, rho, gamma, sig2e)