UNIDAD 6 CAPTULO VII SOLUCIN EN SERIES DE

UNIDAD 6. CAPÍTULO VII. SOLUCIÓN EN SERIES DE POTENCIAS ALREDEDOR DE UN PUNTO SINGULAR REGULAR.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Considérese ahora la solución de la ecuación diferencial lineal de segundo orden coeficientes variables: alrededor de un punto singular regular, x 0. El punto singular regular se supone en el origen, x 0 = 0, lo que simplifica el problema sin pérdida de generalidad, ya que cualquier punto singular x 0 de una ecuación diferencial puede moverse al origen con la transformación t = x x 0. Así mismo se considera el intervalo x > 0 en las explicaciones, para evitar el signo de valor absoluto.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Esto no es una limitación seria, ya que el intervalo x < 0 puede manejarse cambiando la variable t = x y realizando el análisis para t > 0. De la definición se deduce que, si x 0 = 0 es un punto singular regular, entonces las funciones x. P(x) y x 2 Q(x) son analíticas y tienen expansiones de serie de Taylor en ese punto; por lo que pueden expresarse en la forma: y

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Las funciones p(x) y q(x) no requieren incluir un gran número de términos, de hecho, en la mayoría de los casos éstas resultan ser polinomios de grado uno o dos. Por ejemplo, en la ecuación de Bessel de orden cero: Se observa que p(x) = 1 y q(x) = x 2, que son polinomios simples. Al multiplicar la ecuación por x 2 se obtiene: Esta ecuación puede reordenarse en la forma:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. o Una revisión de esta ecuación indica que sus coeficientes aparecen como el producto de los coeficientes de la ecuación de Euler y los de las ecuaciones que tienen soluciones en series de potencias, lo que sugiere una búsqueda de la solución en la forma: es decir, como un producto de la solución de Euler y una solución en series de potencias.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. El procedimiento de buscar una solución de una ecuación diferencial en la forma de un producto de una potencia dada de x y una serie de potencias en x se conoce como método de Frobenius, por el matemático alemán Georg F. Frobenius (1849 -1917). Se verá que, alrededor de un punto singular regular, la ecuación diferencial tiene al menos una solución de la forma de la antes descrita, y otra de forma análoga o modificada, dependiendo de los valores de r. Para determinar el valor de r, se considera la ecuación en su forma general y se sustituyen las expansiones para p(x) y q(x), obteniéndose:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. luego se supone una solución de la forma: con a 0 ≠ 0 y derivándola se obtiene: Sustituyendo estas expresiones en la ecuación diferencial e igualando a cero el coeficiente de xr, que es la potencia mas baja de x, se obtiene:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. o Esta expresión se llama ecuación indicial de la ecuación diferencial y tiene dos raíces reales o complejas, r 1 y r 2, que denominadas exponentes del punto singular regular considerado. Si las raíces son complejas conjugadas, la solución consta de funciones complejas; sin embargo, aún se puede obtener funciones con valor real separando la parte real y la parte imaginaria de la solución.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Cuando las raíces r 1 y r 2 son reales y diferentes, siempre se considera a r 1 como la raíz mayor. Cuando p(x) = x. P(x) y q(x) = x 2 Q(x) son polinomios, las constantes p 0 y q 0 son sus términos independientes; de no ser así, se determinan mediante los siguientes limites: y Ejemplo: Demuestre que x = 0 es un punto singular regular de la ecuación diferencial: y determine las raíces de su ecuación indicial, r 1 y r 2.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Solución: Esta ecuación de segundo orden coeficientes variables expresada en su forma estándar es: entonces: Como se puede ver, Q(x) es analítica en todas partes, pero el denominador de P(x) se hace cero en x = 0, de manera que x = 0 es un punto singular. Sin embargo, este punto es singular regular, ya que p(x) y q(x) son polinomios, es decir, funciones analíticas.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. De aquí se puede observar que p 0 = 2 y q 0 = 0 y al sustituir estos valores en la ecuación indicial r 2 + (p 0 1)r + q 0 = 0, se obtiene r 2 + r = 0 o r(r + 1) = 0, cuyas raíces son r 1 = 0 y r 2 = 1. Luego, la raíz mayor es r 1. Ejemplo: Demuestre que x = 1 es un punto singular regular de la ecuación de Legendre: con a una constante. Mueva este punto singular al origen mediante la transformación t = x 1, y determine las raíces de la ecuación indicial en la ecuación transformada.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Solución: Esta ecuación 2° orden de coeficientes variables expresada en su forma estándar es: entonces: funciones cuyos denominadores se hacen cero en x = 1, por lo que x = 1 y x = 1 son puntos singulares de esta ecuación diferencial. Sin embargo, estos puntos son singulares regulares, ya que para x 0 = 1:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. y son analíticas. Análogamente se puede probar que x 0 = 1 es también un punto singular regular. Para mover el punto singular x 0 = 1 al origen, se aplica la transformación t = x 1 y al sustituir en la ecuación diferencial se obtiene:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. en donde: Como se puede ver, el punto t = 0 es un punto singular, es decir, el punto singular regular x 0 = 1 ha sido desplazado al origen. Previo al cálculo de las raíces de la ecuación indicial, es necesario determinar los valores de p 0 y q 0 y, como p(t) y q(t) no son polinomios, se requiere determinarlos a través de los límites:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Al sustituir en la ecuación indicial r 2 + (p 0 1)r + q 0 = 0, se obtiene r 2 = 0; es decir, las raíces son reales e iguales a cero. Nota: En este caso resulta más conveniente resolver la ecuación: que su equivalente: lo que evita expandir los términos en el denominador usando la expansión del binomio:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. que, en este caso, resulta: y al sustituir en la ecuación específica: que es mucho mas compleja que la ecuación simplificada: La teoría del método de Frobenius es complicada y rebasa el alcance de un texto introductorio sobre ecuaciones diferenciales, los resultados se resumen en el siguiente teorema.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Sea x = 0 un punto singular regular de la ecuación: y sea r el menor radio de convergencia de dos funciones p(x) = x. P(x) y q(x) = x 2 Q(x). Si r 1 y r 2 son las raíces de la ecuación indicial: donde y r 1 > r 2 cuando las raíces son reales, entonces existen dos soluciones linealmente independientes y 1(x) y y 2(x), con un radio de convergencia r. Para x > 0, son de las siguientes formas:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Caso 1: r 1 = r 2 + l (para l positiva no entera). Caso 2: r 1 = r 2 = r. Caso 3: r 1 = r 2 + N (N entero positivo).

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. donde la constante C puede ser cero. Entonces, la solución general de la ecuación diferencial para los tres casos se expresa en la forma: donde las constantes arbitrarias C 1 y C 2 se determinan a través del uso de las condiciones iniciales y de frontera. La demostración de este teorema requiere conocimientos más avanzados sobre ecuaciones diferenciales, por lo que no se presenta en estas notas.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Ejemplo: Resuelva la ecuación de Euler de segundo orden: Usando el método de series alrededor del punto x 0 = 0. Solución: La forma estándar de esta ecuación de segundo orden coeficientes variables es la siguiente: entonces el punto x 0 = 0 es singular regular debido a que p(x) = x. P(x) = 11/4 y q(x) = x 2 Q(x) =1/2 son funciones analíticas (polinomios de grado cero).

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. De lo anterior se puede ver que p 0 = 11/4, q 0 = 1/2 y al sustituirlos en la expresión r 2 + (p 0 1)r + q 0 = 0, se obtiene la ecuación indicial: cuyas raíces son r 1 = 1/4 y r 2 = 2, la raíz mayor es r 1. Como las raíces difieren en un número no entero, de acuerdo con el teorema, la ecuación diferencial tiene dos soluciones linealmente independientes, ambas de la forma:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Derivándola dos veces: y sustituyendo en la ecuación de Euler original: o que se cumple si y sólo si todos sus coeficientes son cero para toda x, esta condición da un número infinito de ecuaciones para determinar los coeficientes de expansión.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. que para k = 0, se reduce a [4 r(r 1) + 11 r 2]a 0 = 0 o 4 r(r 1) + 11 r 2 = 0 (a 0 ≠ 0) y puede reordenarse como: Por tanto se puede concluir que, para k = 0, la relación de recurrencia siempre se reduce a la ecuación indicial. Ahora, la primera solución de la ecuación propuesta, se obtiene usando r = r 1 = ¼ por lo que la relación de recurrencia se transforma en:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. o Entonces, ak = 0 para k = 1, 2, 3, , ya que ningún k ≥ 1 puede hacer cero esta relación. De esta manera, la 1ª solución es: La segunda solución linealmente independiente se obtiene repitiendo el procedimiento para r = r 2 = 2. La relación de recurrencia en este caso se reduce a:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. o Entonces, bk = 0 para k = 1, 2, 3, , ya que ningún k ≥ 1 puede hacer cero esta relación. luego, la segunda solución es: Así, la solución general es:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Ejemplo: Resuelva la ecuación diferencial: alrededor del punto x 0 = 2 para x > 0. Solución: Ecuación diferencial de segundo orden coeficientes variables con x = 2 un punto singular. Para resolverla es necesario desplazar el punto singular al origen mediante la variable t = x 2. Así, la ecuación se transforma en: en donde:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. La ecuación, en su forma estándar, se expresa: en donde t = 0 es un punto singular regular, ya que tanto p(t) = t. P(t) = 1 como q(t) = t 2 Q(t) = t 2 son polinomios, es decir, funciones analíticas. Así p 0 = 1 y q 0 = 2 y, al sustituir en la ecuación indicial r 2 + (p 0 1)r + q 0 = 0, se obtiene r 2 2 = 0; cuyas raíces, difieren en un número no entero (caso 1 del teorema).

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Entonces la ecuación diferencial tiene dos soluciones linealmente independientes de la forma: Derivándola dos veces con respecto a t: y sustituyendo en la ecuación diferencial:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Para igualar el exponente de t en todas las sumas, el índice de la 3ª se desplaza en 1 y se reemplaza k por k 1. Ahora, como la 3ª suma inicia en k = 1, se desarrollan los términos que corresponden a k = 0 en las otras tres sumas, para que todas inicien con k = 1. Ecuación que se satisface para toda x si y sólo si todos los coeficientes son cero.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Esta condición para la potencia más baja de t, junto con el requisito de que a 0 ≠ 0, recupera la ecuación indicial, que se satisface sólo con los dos valores de r obtenidos. El requisito de que se anulen los coeficientes de todas las demás potencias de t da como resultado: o bien: la relación de recurrencia para los coeficientes ak.

U-6. CAP. VII. Para: por tanto: SOLUCIÓN EN SERIES: PUNTOS SINGULARES.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. y, generalizando: Tomando a 0 = 1, la 1ª solución resulta: La 2ª solución linealmente independiente se obtiene repitiendo el procedimiento para:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Con coeficientes bk, la relación de recurrencia en este caso se reduce a: Esta relación difiere de la anterior por un signo, por lo que: y así, si b 0 = 1, la 2ª solución linealmente independiente es:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Entonces la solución general es y = C 1 y 1 + C 2 y 2, pero considerando que t es x 2, tal solución también puede expresarse en la forma: que converge para cualquier x positiva excepto para x = 2.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Ejemplo: Resuelva la ecuación diferencial: alrededor del punto x 0 = 0 para x > 0. Solución: La ecuación diferencial es de segundo orden coeficientes variables que, en su forma estándar, es: en donde x = 0 es un punto singular regular, ya que las funciones p(x) = x. P(x) = (x 2 3) y q(x) = x 2 Q(x) = 4 son analíticas.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. En estas relaciones se observa que p 0 = 3 y q 0 = 4, con lo que la ecuación indicial r 2 + (p 0 1)r + q 0 = 0 resulta en r 2 4 r + 4 = 0 o (r 2)2 = 0 y sus raíces son r 1 = r 2 = 2. De acuerdo con el Caso 2 del teorema, la ecuación tiene dos soluciones linealmente independientes; la 2ª contendrá un término logarítmico y la primera es de la forma: Derivando dos veces con respecto a x, se tiene:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. y sustituyendo en la ecuación diferencial: El exponente de x se iguala en todas las sumas, corriendo en 2 los índices de la 2ª y reemplazando k por k 2. Ahora, como la 2ª suma inicia con k = 2, se desarrollan los términos que corresponden a k = 0 y k = 1 y el resultado se combina con las otras sumas (iniciando con k = 2) para obtener:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Esta ecuación se satisface para toda x si y sólo si cada uno de los coeficientes es cero, condición que, para la potencia más baja de x y el requisito de que a 0 ≠ 0, recupera la ecuación indicial, entonces: Para r = 2 resultan a 1 = 0 y la relación de recurrencia:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. que expresa un coeficiente en términos del 2° anterior. Por tanto, todos los coeficientes con índice par se expresan en términos de a 0, y los de índice impar, de a 1, que es cero. Entonces, a 1 = a 3 = a 5 = = 0 y los coeficientes con índice par se expresan como: y en general: Al tomar a 0 = 1, la primera solución resulta:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. La 2ª solución linealmente independiente, y 2, se determina, de acuerdo con el teorema, en la forma: Derivando dos veces con respecto a x, se tiene:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. al sustituir ambas derivadas en la ecuación diferencial y ordenarla, se obtiene: 1) El término logarítmico se anula porque su coeficiente es la ecuación diferencial y y 1 es una solución. Esto ocurre siempre para el caso 2.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. 2) Ahora, los términos con y 1 y su derivada son: 3) Por último, las sumas que contienen a bn se pueden simplificar a la forma:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. De esta manera se tiene: o el primer término (k = 0) de la primera suma es cero, por lo que:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. un reordenamiento de la segunda suma da: así, al combinar éstas dos se obtiene: luego, al combinar y ordenar las sumas con coeficiente bk resulta:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. de esta manera se obtiene la ecuación: que se cumple si y sólo si todos los coeficientes son cero. Aplicando esta condición a los dos primeros términos se obtiene b 1 = 0 y b 2 = ¼. Luego, si se considera que los términos de potencia impar proceden sólo de la segunda suma (siendo valores impares de k) se tiene que k 2 bk + kbk 2 = 0 para k impar. Pero b 1 = 0, por lo que bk = 0 para k impar. Es decir:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Entonces la 2ª suma puede modificarse para excluir todos los términos de potencia impar, reemplazando k por 2 k, lo que resulta en: o y al anular los términos entre corchetes para toda k se obtiene la siguiente relación de recurrencia:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Recordando que b 2 = ¼ , se tiene: y en general: por tanto: La solución general es de la forma y = C 1 y 1 + C 2 y 2, y convergirá para cualquier x > 0.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Ejemplo: Resuelva la ecuación diferencial: alrededor del punto x 0 = 0 para x > 0. Solución: La ecuación diferencial es de segundo orden coeficientes variables que, en su forma estándar, es: en donde x = 0 es un punto singular regular, ya que las funciones p(x) = x. P(x) = 2 y q(x) = x 2 Q(x) = x 2 son analíticas.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. de las relaciones anteriores, p 0 = 2 y q 0 = 0, con lo que la ecuación indicial r 2 + (p 0 1)r + q 0 = 0 resulta r 2 + r = 0 o r(r + 1) = 0 y sus raíces son r 1 = 0 y r 2 = 1. De acuerdo con el Caso 3 del teorema, la ecuación tiene dos soluciones linealmente independientes; la primera es de la forma: En este caso se adoptará el menor valor de r = r 2 = 1. Obtener ceros para los coeficientes indicará que esto es incorrecto. Obtener una solución con dos constantes arbitrarias indicará que se tienen las dos soluciones.

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Derivando dos veces la solución propuesta: y sustituyendo en la ecuación diferencial: Para igualar el exponente en todas las sumas, se corre en 2 el índice de la última, reemplazando k por k 2:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Así las potencias de x son iguales, pero la segunda suma inicia con k = 2. Para hacer que las demás también lo hagan, se desarrollan los términos para k = 0 y k = 1, y se combinan las sumas, Esta ecuación se satisface para toda x si y sólo si cada uno de los coeficientes es cero, condición que, para la potencia más baja de x y el requisito de que a 0 ≠ 0, recupera la ecuación indicial, entonces:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Para r = r 2 = 1 la primera relación da 0 ∙ a 1 = 0, lo cual indica que a 1 es arbitrario, al igual que a 0. La segunda relación da la relación de recurrencia: que expresa un coeficiente en términos del segundo coeficiente anterior; así, los coeficientes con índice par se expresarán en términos de a 0 y los de índice impar de a 1. Los coeficientes de valores pares de k son:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. y generalizando se tiene: Los coeficientes para valores impares de k son: generalizando:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Finalmente, al sustituir en la solución propuesta: se obtiene la solución general: o bien:

U-6. CAP. VII. SOLUCIÓN EN SERIES: PUNTOS SINGULARES. Dado que estas soluciones son linealmente independientes, se puede concluir que ésta es una solución general de la ecuación diferencial planteada. Al repetir el análisis para r = r 1 = 0, se obtiene como la 1ª solución linealmente independiente: que es una de las soluciones encontradas. Al continuar con r = r 2 = 1 y la 2ª solución con término logarítmico, resulta C = 0. Esta solución convergirá para x > 0.