איך פותרים משוואות דיפרנציאליות רגילות ליניאריות אי-הומוגניות מסדר שני?

הכוונה שלי היא למשוואות מהצורה:
y''+p(x)y'+q(x)y=r(x)
למשל:
y''+(x^2)y'-4xy=1/x
אין הסבר בויקיפדיה. הדבר הכי קרוב שמצאתי הוא הסבר של פתרון משוואות דיפרנציאליות רגילות ליניאריות הומוגניות מסדר שני שהמקדמים בהן הם קבועים.
רוב המחשבונים שמצאתי באינטרנט לא מצליחים לפתור משוואות מהסוג הזה. מצאתי אחד שנותן פתרון אבל לא מוכן לרשום את הדרך (בדרך כלל הוא נותן דרך). כנראה הוא מצא קירוב בעזרת שיטות ממוחשבות.