Python somme de riemann

[novicepython]

Bonjour tous le monde,
J'ai fais un exercice de maths et à la fin on me demande d'écrire un programme python qui permet de récupérer la somme de riemann ( f est connue )
J'ai écris ce programme :

from math import sqrt
def f(x):
    return sqrt(1-x**2)
n=input("choisir un nombre entier")
di = 1/n
S=0
for k in range(0,n):
   S +=f(k*di)*di
return(S)

il m'affiche qu'il y a une erreur avec return mais je ne comprends pas pourquoi et si ce problème était réglé je pense que le programme tient la route non ?
Merci d'avance
PS: dans mon exo a=0 et b=1 c'est pour j'ai écris 1/n
-----

[eudea-panjclinne]

1/n renvoie le quotient d'une division euclidienne de deux entiers. Pour palier à cela écrivez 1.0/n vous aurez alors un division entre réels.
A la fin remplacer return(S) par print(S).