No estoy seguro si usted tiene su corazón puesto en las soluciones exactas, pero, ¿qué acerca de un perturbativa de la solución? Yo sólo puedo hablar en relativamente términos generales, por ahora, pero espero que este esquema es lo suficientemente claro para orientar sobre qué hacer.
La reescritura de la primera ecuación, de la que conocemos la solución
$$x y''+2 x y' + x^2 y = 0$$
Podemos escribir la solución deseada para esta ecuación como $y_0(x)$. Escribe una forma general de la solución anterior; Mathematica, sin embargo, escribe la forma más compacta como
$$y_0(x) = e^{-x} (a \, \mathrm{Ai}{[e^{i \pi/3} (1-x)]} + b \, \mathrm{Bi}{[e^{i \pi/3} (1-x)]}) $$
$\mathrm{Ai}$ $\mathrm{Bi}$ siendo las funciones de Airy y las constantes $a$ $b$ determinado por las condiciones iniciales no se especifica.
Podemos reescribir la 2ª ecuación de la siguiente manera:
$$x y''+2 x y' + x^2 y = (4 x+1) y'$$
El primer paso en un enfoque perturbative es introducir una "pequeña" parámetro " $\epsilon$ que representa una desviación de la solución inicial. Es decir, podemos escribir la ecuación anterior como
$$x y''+2 x y' + x^2 y = \epsilon (4 x+1) y'$$
Este es un habitual de perturbación debido a la perturbación no se hace en la mayor de derivados (es decir, no hay nuevas soluciones serían introducidas por la perturbación). La solución de $y(x)$ puede entonces escribirse como
$$y(x) = \sum_{n=0}^{\infty} \epsilon^n y_n(x)$$
sin saber muy bien si la suma converge. El $y_n(x)$ satisfacer a una relación de recurrencia:
$$x y_n''+2 x y_n' + x^2 y_n = \epsilon (4 x+1) y_{n-1}'$$
donde, para $n \ge 1$, $y_n(0) = 0$, $y_n'(0) = 0$.
Por supuesto, usted no quiere tener que calcular más de un par de términos. También es necesario determinar si el anterior de la serie converge o diverge; si diverge, entonces usted necesita para averiguar cómo muchos términos que usted puede calcular antes de la aproximación ya no funciona. Finalmente, por supuesto, para obtener una aproximación a la solución,$\epsilon = 1$.