Skip to main content
se añadieron 1110 caracteres en el cuerpo
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68

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)

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)
se añadieron 453 caracteres en el cuerpo
Origen Enlace
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. La 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í

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. La gráfica que nos crea es esta:

introducir la descripción de la imagen aquí

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í

se eliminaron 12 caracteres en el cuerpo
Origen Enlace
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)')
    area_text = 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. La gráfica que nos crea es esta:

introducir la descripción de la imagen aquí

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)')
    area_text = 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. La gráfica que nos crea es esta:

introducir la descripción de la imagen aquí

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. La gráfica que nos crea es esta:

introducir la descripción de la imagen aquí

cuerpo editado
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68
Loading
Origen Enlace
FJSevilla
  • 62.1k
  • 7
  • 43
  • 68
Loading