Une version un peu améliorée du script de base de Numworks “polynomial.py”
from math import * # racines(a,b,c) determine les solutions de l'equation a*x**2+b*x+c=0 def racines(a,b,c): delta = b*b-4*a*c print("delta=",delta) if delta == 0: #donne la valeur exacte de la solution print ("x_0=","-",b,"/ 2 *",a,"=",-b/(2*a)) elif delta > 0: x_1 = (-b-sqrt(delta))/(2*a) x_2 = (-b+sqrt(delta))/(2*a) #donne la valeur exacte de la solution print ("x_1=","( -",b,"- sqrt(",delta,")",")","/ 2 *",a,"=",x_1) print ("x_2=","( -",b,"+ sqrt(",delta,")",")","/ 2 *",a,"=",x_2) else: print ("pas de solutions reelles")