# -*- coding: utf-8 -*- """ Created on Wed Jun 21 18:15:32 2017 Exercice 2-4 : calcul de J_0 par le developpement en série @author: grivet """ import scipy.special as ss tol = 1.0e-7; imax = 40 n = 0 x = float(input("valeur de x: ")) somme = 0.0; terme = 1; i = 0 xx = x*x while (abs(terme) > tol) & (i < imax): somme = somme + terme terme = -xx/(4*(i+1)*(n+i+1))*terme i = i + 1 # tt(i) = terme print("valeur calculée: %12.9f: " % somme) print("valeur donnée par Scipy: %12.9f" % ss.jn(0,x));