configurar las odas en python
¿Cómo configurar las siguientes odas con las condiciones iniciales correspondientes en python?
x'(t) =x(t) - y(t) - e^t
y'(t) =x(t) + y(t) + 2e^t
conx(0)= -1
yy(0)= -1
y0 <= t <= 4
Lo siguiente es lo que tengo hasta ahora:
def f(u, t):
x, y = u
return [x+y-e**t, x+y+2*e**t]
x0, y0 = [-1.0,-1.0]
t = numpy.linspace( 0,4,50 )