[Maxima-lang-es] Problema al resorver una EDO con comando "problema de valor inicial (1)"

[E. S. <ser...@ho...>]

Hola:

Estoy tratando de resolver la EDO

y ’= y (xy ^ 3-1)



con valor inicial y (0) = 1



ec: 'diff (y, x) = y * (x * y ^ 3-1);



Con el comando ode2 obtengo la solución general:





oda2 (ec, y, x);



y =% e ^ ​​(- x) / (((3 * x + 1) *% e ^ ​​(- 3 * x)) / 3 +% c) ^ (1/3)



Pero para encontrar la solución particular, utilizo "Problema de valor inicial":



ic1 (y =% e ^ ​​(- x) / (((3 * x + 1) *% e ^ ​​(- 3 * x)) / 3 +% c) ^ (1/3), x = 0, y  = 1);



y obtengo la solución particular:



y = 3 ^ (1/3) / (3 *% e ^ ​​(3 * x) + 3 * x + 1) ^ (1/3)



que no satisface el valor inicial y (0) = 1.

(se obtiene y (0) = 3 ^ (1/3) / 4 ^ (1/3) en lugar de y (0) = 1)



¿Es esto un error del comando ic1 o estoy haciendo algo mal?



¡Gracias!