Lo que veo que pasa entre espejos:

KdV equation u_t + u_xxx + 6uu_x= 0

La solución de la KdV para 2 solitones crea esto:

Y lo que a mi me interesa ver en detalle es lo que pasa cuando a uno de los solitones, les van llegando los otros :

Pero la aproximación de la solución debe estar mal calculada, ya que el programita que hice no da bien:

-----------------------------------------------------------------------------

zoomx = 20
zoomy = .8

FOR t = -2 TO 6 STEP .01
 

FOR p = 1 TO 100000
NEXT p
CLS

FOR x = -10 TO 10 STEP .1

k1 = 2
k2 = 1
x1 = -8
x2 = -2

eta1 = k1 * (x - x1) - k1 ^ 3 * t
eta2 = k2 * (x - x2) - k2 ^ 3 * t
c12 = (k1 - k2) ^ 2 / (k1 + k2) ^ 2
f = 1 + EXP(eta1) + EXP(eta2) + c12 * EXP(eta1 + eta2)

IF x > 320 THEN x = 320
IF f > 240 THEN f = 240

IF x < -320 THEN x = 0
IF f < -240 THEN f = 0

PSET ((320 + x * zoomx), (240 - f * zoomy)), 15

NEXT x
NEXT t
-----------------------------------------------------------------------------

Así que a ver si alguien me tira los códigos correctos, ya que este simple resultado es muy útil para la docencia y la investigación, y a mas de uno/a le puede venir bien... el problema creo que pasa por malinterpretar esto:

                                                      u = 2*(Log f)_{xx}    que significa exactamente???

La que si funciona es:

Para x entre -17 y 17     y para t entre -5 y 5     con b1= 0.5 y b2= 1

                                                                                                                            Cristián Antiba

Dejar Mensajes !!!