Newton-Picard method