import numpy as np from matplotlib import pyplot as plt def dF(x,y): #y' = f(x,y) f=np.exp(-3*x)-2*y return f def fun(x): #Resultado analitico fun = 2*np.exp(-2*x)-np.exp(-3*x) return fun x0 = 0 y0 = 1 xf = 5 n = 100 h = (xf-x0)/n x = np.array([x0]) y = np.array([y0]) for i in range(1,n+1): yn=y0+h*dF(x0,y0) x0 = x0+h y0 = yn x = np.append(x,x0) y = np.append(y,y0) xp = np.linspace(0,5,100) plt.figure() plt.plot(x,y) plt.plot(xp,fun(xp)) plt.legend(['Euler','Analitico']) plt.grid() plt.show()
¡No tienes idea de cuanto me ayudaste!
Excelente video, super completo y claro👌🏻
Muchísimas gracias ✨
muy buen video!
import numpy as np
from matplotlib import pyplot as plt
def dF(x,y): #y' = f(x,y)
f=np.exp(-3*x)-2*y
return f
def fun(x): #Resultado analitico
fun = 2*np.exp(-2*x)-np.exp(-3*x)
return fun
x0 = 0
y0 = 1
xf = 5
n = 100
h = (xf-x0)/n
x = np.array([x0])
y = np.array([y0])
for i in range(1,n+1):
yn=y0+h*dF(x0,y0)
x0 = x0+h
y0 = yn
x = np.append(x,x0)
y = np.append(y,y0)
xp = np.linspace(0,5,100)
plt.figure()
plt.plot(x,y)
plt.plot(xp,fun(xp))
plt.legend(['Euler','Analitico'])
plt.grid()
plt.show()