import numpy as np import matplotlib.pyplot as plt def f(x): # Esta función evalúa el polinomio. Puedes modificarla según tus necesidades. return eval(polynomial) def suma_riemann(a, b, n): delta_x = (b - a) / n suma = 0 x_values = [] for i in range(n): x = a + i * delta_x suma += f(x) * delta_x x_values.append(x) return suma, delta_x, x_values # Solicitar entradas al usuario polynomial = input("Ingresa el polinomio (por ejemplo, 'x**2 + 2*x + 1'): ") a = float(input("Ingresa el límite inferior de integración (a): ")) b = float(input("Ingresa el límite superior de integración (b): ")) n = int(input("Ingresa el número de rectángulos (n): ")) # Calcular la suma de Riemann resultado, delta_x, x_values = suma_riemann(a, b, n) # Mostrar resultados print(f"\nSuma de Riemann: {resultado}") print(f"Delta x: {delta_x}") print(f"Número de rectángulos: {n}") # Graficar el polinomio y los rectángulos x = np.linspace(a, b, 100) y = f(x) plt.plot(x, y, label='Polinomio') for i in range(n): plt.bar(x_values[i], f(x_values[i]), width=delta_x, alpha=0.5, align='edge', color='orange') plt.title('Suma de Riemann') plt.xlabel('x') plt.ylabel('f(x)') plt.legend() plt.grid() plt.show()

Please keep input under 1000 characters