# -*- 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));