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 f[0]=0
, f[10]=1
DSolve
, but I need a numerical solution andNDSolve
does cannot solve the above equation as the image shows. $\endgroup$DSolve
does give a solution in terms ofErf
andHeavisideTheta
. I don't know much about the numerical algorithms, sorry. $\endgroup$