Al terminar este artículo podrás resolver TODAS tus ecuaciones diferenciales lineales o NO lineales de 1er oeden, con valores iniciales, mediante el Método de Euler con SAGEMATH.

El método de Euler es implementado en SAGEMATH el cual es un lenguaje de programación de alto nivel, construído sobre python y otros lenguajes de acceso libre, es decir; SAGEMATH es software libre y se utiliza básicamente para la simulación científica. La programación y simulación con dicho lenguaje es muy sencilla, incluso podrás simular tus ejercicios aquí mismo.

Se ponen ejemplos resueltos con SAGEMATH de los ejercicios vistos en artículos anteriores. Ver los enlaces específicos para cada artículo en cada ejercicio resuelto.

Para saber cómo editar (utilizar) las celdas de SAGEMATH, ve al siguiente enlace: Simulación con SAGEMATH, da click aquí

Para entender a detalle el código de SAGEMATH para resolver ecuaciones diferenciales con valores iniciales mediante el método de Euler ve la siguiente presentación: De donde sale el método de Euler.

El código para el Método de Euler escrito en SAGEMATH es el siguiente:

CÓDIGO PARA EL MÉTODO DE EULER CON SAGEMATH

### Metodo de Euler
def Euler(fun, a, b, N, y0):
    h = (b - a)/N
    x = [a]
    y = [y0]

    for k in range(N):
        x.append(x[k]+h)
        y.append(y[k]+(h)*fun(x[k], y[k]))
    return list(zip(x, y))

Para Utilizar las celdas para simular otras ecuaciones diferenciales solo es necesario editar el apartado de «datos iniciales» y «Solución numérica para h=…». Al final de este artículo les dejo una versión simplificada de esta celda para que la modifiquen y puedan simular otras ecuaciones diferenciales de primer orden con valores iniciales. 😉

Sigue leyendo →