А можно как-то напрямую подставить тот вид решения в уравнение и проверить?
Очень просто: подставить и проверить. Вне начала координат получается (как и должно было бы быть) тождественный ноль, а скачок первой производной в нуле (и только в нуле) как раз и обеспечивает после повторного дифференцирования дельта-функцию.
Это если нужно ехать. А если нужно шашечки, то да, надо записать с помощью функции Хевисайда и сугубо формально дифференцировать в рамках теории обобщённых функций. Если без омеги (она всё равно непринципиальна), то выйдет так:

