He estado probando la IA de qwen3, la última que ha salido, y me ha venido bien en este caso para resumir la teoría para hallar la solución general a una ecuación diferencial de la forma \( y''(x)=c\cdot y(x)+f(x) \) para \( c \) constante. Éste es el resumen que me ha dado (esta técnica la había olvidado):
...
La ecuación es un ejemplo concreto de EDO lineal completa con coeficientes constantes, en este caso de orden dos. Como dices que tenías olvidado el método... vamos a recordarlo...

En general, para resolver una EDO lineal completa de orden \( n \) con coeficientes constantes,
\( y^{n)}(x)+a_{n-1}y^{n-1)}(x)+\cdots +a_1y'(x)+a_0y(x)=f(x) \)
la solución general viene dada como suma de una solución particular de la completa, \( y_p \), más la solución general, \( y_h \), de la EDO lineal homogénea asociada.
Para obtener \( y_h \) se resuelve la
ecuación característica asociada a la EDO y en función de sus raíces y de sus órdenes de multiplicidad se obtienen \( n \) soluciones linealmente independientes de la EDO lineal homogénea y se expresa \( y_h \) como combinación lineal de todas ellas:
\( y_h=C_1y_1+\cdots +C_ny_n \).
Para hallar una solución particular \( y_p \) de la EDO lineal completa, hay varios métodos. Uno de ellos es el llamado
variación de las constantes y consiste en buscar una solución particular de la EDO lineal completa de la forma \( y_p(x)=C_1(x)y_1(x)+\cdots +C_n(x)y_n(x) \).
Resolviendo el sistema
\( \left\{\begin{array}{lcl}C_1'y_1+\cdots +C_n'y_n&=&0\\C_1'y_1'+\cdots +C_n'y_n'&=&0\\\vdots&\vdots &\vdots\\C_1'y_1^{n-1}+\cdots +C_n'y_n^{n-1}&=&f\end{array}\right. \) se hallan las funciones \( C_j,\,j=1,\ldots, n \) y se concluye expresando la solución general de la EDO lineal completa como \( \boxed{y=y_h+y_p} \)
Hay otros métodos para hallar \( y_p \), como el llamado
método de los coeficientes indeterminados, pero éste solo es válido para ciertos tipos de funciones \( f \),
Saludos