Tengo que programar una suma de Reimann. Tenía la idea de hacerla con un ciclo `for` que fuera calculando los valores de las areas de los rectangulos para guardarlo en una variable.
Pero no sé cómo podría sumar los resultados obtenidos del ciclo `for`. El código es en Python 

    import math
    import numpy as np

    print("- Calculadora de la suma de Reimann para la función seno -")
    print("")
    print("Por favor, especifique los siguientes parametros para hacer sus 
    cálculos: ")
    print("")
    a=float(input("1.- El punto a de su intervalo: "))
    b=float(input("2.- El punto b de su intervalo: "))
    n=float(input("3.- El número n de rectángulos con el cual desea hacer 
        su cálculo: "))

    dn=(b-a)/n

    def f(x):
        y=np.sin(x)
        return (y)
    
    def RiemannINF( a, b, f, dn):
        for i in range ( a, b, dn):
            ap=... 

Es lo que llevo de código.
No estoy muy seguro de cómo poder plantear la suma con ciclos `for` y me vendría bien una mano.