http://www.perlmonks.org?node_id=771382

y"-4xy'+(4x^2-1)y=-3e^(x^2)sin(x);