Rozważania podzielmy na kroki.
1) Najpierw trzeba znaleźć (tu: raczej odgadnąć) jedno niezerowe rozwiązanie
\(\displaystyle{ \varphi_1}\) równania jednorodnego
\(\displaystyle{ (1-x)y''+xy'-y=0}\). Łatwo zauważamy, że można przyjąć
\(\displaystyle{ \varphi_1(x)=e^x}\).
2) Sprowadźmy równanie jednorodne do postaci normalnej:
\(\displaystyle{ y''=\frac{x}{x-1}y'-\frac{1}{x-1}y}\). Wyznaczymy rozwiązanie
\(\displaystyle{ \varphi_2}\) tego równania liniowo niezależne z
\(\displaystyle{ \varphi_1}\). W tym celu pomocnicze równanie liniowe (pierwszego rzędu) postaci
\(\displaystyle{ z'=\frac{x}{x-1}-2\frac{\varphi_1'(x)}{\varphi_1(x)}z=-2z+\frac{x}{x-1}}\).
Niech
\(\displaystyle{ \psi}\) będzie rozwiązaniem tego równania (metoda rozwiązywania równań liniowych pierwszego rzędu - patrz np. 100572.htm ). Wówczas
\(\displaystyle{ \varphi_2(x)=\varphi_1(x)\int\psi(x)dx}\) (gdzie przez całkę rozumiemy dowolną z funkcji pierwotnych funkcji
\(\displaystyle{ \psi}\)).
3) Teraz wystarczy znaleźć jedno szczególne rozwiązanie
\(\displaystyle{ \varphi_0}\) równania wyjściowego (niejednorodnego). Funkcję tę można wyznaczyć np. metodą wariacji stałych, tj. przyjąć, że
\(\displaystyle{ \varphi_0(x)=a(x)\varphi_1(x)+b(x)\varphi_2(x)}\) i różniczkując dwukrotnie stronami oraz biorąc pod uwagę, że
\(\displaystyle{ \varphi_0}\) ma spełniać dane równanie, wyznaczyć niewiadome funkcje
\(\displaystyle{ a, b}\).
4) Ogół rozwiązań równania stanowią funkcje postaci
\(\displaystyle{ \varphi_0(x)+C_1\varphi_1(x)+C_2\varphi_2(x)}\), gdzie
\(\displaystyle{ C_1, C_2}\) są dowolnymi stałymi.