The simulated Newton Raphson method