Romeo y Julieta

Suponga que en cualquier tiempo t podemos medir el amor u odio que siente Romeo hacia Julieta r(t) y el amor u odio que Julieta siente por Romeo j(t). Los valores positivos de estas funciones indican amor y valores negativos, odio. De manera simple, supongamos que el cambio del amor de Romeo por Julieta es una fracción del amor que tiene hacia ella mas un fracción del amor que tiene ella hacia él. De manera similar, el cambio del amor de Julieta por Romeo es una fracción del amor que siente por Romeo más una fracción del amor que siente Romeo por ella. Estas suposiciones generan las siguientes ecuaciones:
r'(t) = a r(t) + b j(t),
j'(t) = c r(t) + d j(t),
donde a, b, c, d son parámetros con valores numéricos.
Nota: El modelo aquí mostrado se obtuvo del libro "The Mathematics of Marriage" de Gottan, Murray, Swanson, Tyson y Swanson, Ed: MIT, 2002, páginas 111-114.

Vamos a suponer el caso en que a = d = 0, esto es, que el cambio del amor de Romeo por Julieta sólo depende del amor que ella le tiene a él y que el cambio que sufre el amor de Julieta por Romeo sólo depende del amor que él le tiene a ella.
Sin entrar en debates sexistas, supongamos que b = 1 y c = -1, lo que significa que el amor de Romeo por Julieta se incrementará tanto como Julieta lo ame, pero el amor de Julieta sólo se incrementará si el amor de Romeo hacia ella decrece. De esta forma tenemos que si derivamos la primera ecuación y sustituimos la segunda en ésta tenemos que
r''(t) = b j'(t) = b c r(t) = -r(t).La transformada de Laplace de la ecuación anterior es
s²R(s) - s r(0) - r'(0) = -R(s),
por lo que
R(s) = (s r(0) + r'(0)) / (s²+1).Nuevamente, sin entrar en debates, supongamos que Romeo fué el primero en sentir algo por Julieta, algo así como amor a primera vista, en este caso r(0) ≠ 0. Por otra parte, Julieta al no saber de Romeo no podia sentir ni odio ni amor, i.e. j(0) = 0. Observe que de la primera ecuación y con nuestros parámetros, tenemos que r'(t) = j(t), por lo que j(0) = r'(0) = 0. Con esto tenemos las dos condiciones iniciales que necesitamos. Reescribiendo, tenemos queR(s) = s r(0) / (s²+1).De esta forma tenemos que la solución es
r(t) = r(0) cos (t).Nuevamente, según la primera ecuación y nuestros parámetros, tenemos que r'(t) = j(t), es decir
j(t) = -r(0) sen(t).En la figura se muestra el campo vectorial del par de ecuaciones descritas al comienzo y la solución encontrada a partir de las condiciones iniciales dadas. Se puede apreciar la oscilación de los sentimientos que se tienen entre ellos, algo así como: los dos andan muy bien (cuadrante I), Julieta ya está molesta con Romeo (cuadrante IV), los dos están molestos (cuadrante III), ahora Romeo es el que está molesto con Julieta (cuadrante II) y se repite. Les suena familiar?
Observación: Hubiera tenido tanto éxito la obra de Shakespeare si los personajes se hubieran comportado como aquí se describió? Qué parámetros serían necesarios para que la situación se diera como en la obra?