It seems that Mathematica can solve an initial value problem with a
DiracDelta, but not a boundary value problem with a
DiracDelta. Is there a workaround?
Consider for example a differential equation with this form:
f''[x]+x f'[x]== x^2 DiracDelta[x-1]
in the range
0,<x<10 with the boundary conditions