Plot solution for y' + 2xy = 1
f = desolve(diff(y,x) + 2xy - 1, y, ics=[0,0]); f plot(f) # error message ... unable to simplify to float approximation.
Tried plotting real_part of solution: -1/2Isqrt(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.