Skip to main content
4 de 5
se añadieron 453 caracteres en el cuerpo
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

Claro que se puede, con Python el cielo es el límite... XD.

Bromas aparte, no especificas que suma de Riemann quieres aplicar (derecha, izquierda, máximo, mínimo). Voy a usar la suma izquierda de acuerdo a tu ejemplo:

Para la gráfica uso matplotlib y NumPy para los arrays:

import numpy as np
import matplotlib.pyplot as plt


def riemannplot(f, a, b, ra, rb, n):
    # f es la función 
    # a y b son los limites del eje x para graficar la funcion f
    # ra y rb son los limites del intervalo en el eje x del que queremos calcular la suma
    # n es el numero de rectangulos que calcularemos

    atenuacion = (b-a)/100
    x = np.arange(a, b+atenuacion, atenuacion)
    plt.plot(x, f(x), 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()



def f(x):
    return x**2

riemannplot(f, 0, 1.1, 0, 1, 10)

No se si tienes que hacerlo mediante SymPy obligatoriamente o no. No obstante, la idea básica es la misma uses el metodo que uses para graficar:

❶ Dividir el intervalo en n partes iguales, obteniendo los valores de x que separan cada rectangulo.

❷ Para cada triangulo calcular su altura, para lo cual basta con pasar cada valor de x obtenido antes a la función.

❸ Ya solo queda graficar el histograma usando para ello las alturas y el ancho de cada barra que es el incremento de x (longitud del intervalo/n)

La gráfica que nos crea es esta:

introducir la descripción de la imagen aquí

FJSevilla
  • 62.1k
  • 7
  • 43
  • 68