### Plot solution for y' + 2xy = 1

f = desolve(diff(y,x) + 2*x*y - 1, y, ics=[0,0]); f
plot(f) # error message ... unable to simplify to float approximation.

Tried plotting real_part of solution: -1/2*I*sqrt(pi)*e^(-x^2)*erf(I*x),
but get same error message.

Tried using list_plot, but error about symbolic expression.
Haven't been able to get implicit_plot or sol.simplify_full to work.

Thanks for any hints.