Si no quieres usar una función auxiliar puedes pasar tu función como un string a riemannplot()
usando SymPy
para ello:
import numpy as np
import matplotlib.pyplot as plt
from sympy import S, symbols
from sympy.utilities.lambdify import lambdify
def riemannplot(f, a, b, ra, rb, n):
x = symbols('x')
f = lambdify(x, S(f),'numpy')
atenuacion = (b-a)/100
xs = np.arange(a, b+atenuacion, atenuacion)
plt.plot(xs, f(xs), color='green')
delta_x = (rb-ra)/n
riemannx = np.arange(ra, rb, delta_x)
riemanny = f(riemannx)
riemann_sum = sum(riemanny*delta_x)
plt.bar(riemannx,riemanny,width=delta_x,alpha=0.5,facecolor='orange')
plt.xlabel('x')
plt.ylabel('f(x)')
plt.title('Suma izquierda de Riemann para f(x)')
plt.figtext(0.1,0.02, "Suma de Riemann: " + str(riemann_sum), color='r')
plt.show()
Ahora puedes llamar a riemannplot()
de esta forma:
riemannplot('x**2', 0, 1.1, 0, 1, 10)