Cómo resolver la ecuación diferencial [matemáticas] y ‘(x) = x + y (x) + xy (x) + y ^ 3 (x) [/ matemáticas]

La ecuación dada es una ecuación diferencial ordinaria no lineal de primer orden. No es fácil encontrar una solución simbólica para esta ecuación.

Una forma de resolver la ecuación diferencial dada es obtener una solución numérica con condiciones iniciales específicas. Esto se puede hacer utilizando la función incorporada de Mathematica NDSolveValue [].

Por ejemplo, escribiendo el código

ysol1 = NDSolveValue [{y ‘[x] == x + y [x] + x * y [x] + y [x] ^ 3, y [0] == 3},
y, {x, 0, 30}]

calcula la función de interpolación para la condición límite especificada.

Las soluciones se pueden trazar escribiendo el código:

Trazar [ysol1 [x], {x, 0, 25}, PlotRange -> Todos]

La trama resultante es la siguiente:

Se pueden encontrar y trazar varias soluciones con diferentes condiciones iniciales o límite. Podemos encontrar soluciones para las condiciones de contorno [matemáticas] y (0) = k [/ matemáticas] escribiendo el código:

ysoln0 = Tabla [
NDSolveValue [{y ‘[x] == x + y [x] + x * y [x] + y [x] ^ 3, y [0] == k},
y, {x, 0, 30}], {k, 1, 10}

Esto calcula las funciones de interpolación para diez soluciones con diez condiciones de contorno diferentes.

Las soluciones se pueden trazar escribiendo:

Trazar [Evaluar [Tabla [ysoln0 [[k]] [x], {k, 1, 10}]], {x, 0, 25},
ImageSize -> Large, PlotRange -> Automatic]

La trama resultante es (haga clic en la imagen a continuación para ampliarla):

Se pueden encontrar y trazar múltiples soluciones con múltiples condiciones iniciales o límite [matemáticas] y (i) = j [/ matemáticas].

Mecanografía:

ysoln = Tabla [
NDSolveValue [{y ‘[x] == x + y [x] + x * y [x] + y [x] ^ 3, y [i] == j},
y, {x, 0, 30}], {i, 1, 3}, {j, 1, 3}]

calcula las funciones de interpolación para nueve soluciones con nueve condiciones de contorno diferentes.

Luego escribiendo el código:

Trazar [Evaluar [Tabla [ysoln [[i, j]] [x], {i, 1, 3}, {j, 1, 3}]], {x, 0,
30}, Tamaño de imagen -> Grande]

proporciona la siguiente gráfica de las soluciones (haga clic en la imagen para ampliarla):

Descargo de responsabilidad: soy un estudiante universitario y solo he tenido un semestre de diferencias, por lo que no soy una gran autoridad aquí.

Esto es lo que creo que es la mejor respuesta a este problema:

Aparentemente, la mayoría de las ecuaciones diferenciales ‘no lineales’ no pueden resolverse explícitamente. Esta ecuación es ‘no lineal’ porque contiene el término y ^ 3. Si el término dependiente (en este caso y) tiene un operador, como ln (y), sqrt (y), cos (y), etc., entonces toda la función es ‘no lineal’.

De todos modos, dado que esta es una extraña función no lineal para la que no parece haber una solución obvia, mi solución sería simplificarla un poco y luego integrar ambos lados. Para muchas ecuaciones diferenciales, esta parece ser la mejor solución.

¿No significa esto que y es una función derivada principal de x? Y … que y también varía de acuerdo con el valor de x?

La función y ‘(x) es igual / igual que x + y + xy + y ^ 3. ¿Son proporcionales entre sí?