# -*- coding: utf-8 -*- """ Created on Thu Dec 18 10:50:10 2014 @author: grivet Exercice 5.1dSP : résolution d'une équation non-linéaire par la méthode de Newton, fournie par Scipy """ from pylab import * from scipy.optimize import * def fn(x): return x - 0.2*sin(x) - 0.5 def dfndx(x): return 1 - 0.2*cos(x) x = eval(input('approximation de la racine: ')) z = newton(fn,x,fprime=dfndx,args=(),tol=1.0e-08,maxiter=50,fprime2=None) print ('\t z fn(z)') print(' %14.6f\t%14.6f' %(z,fn(z)))