ex8.3

import sympy as sp

sp.var('x'); y=sp.Function('y')
eq=y(x).diff(x)+2y(x)-2x**2-2*x
s=sp.dsolve(eq, ics={y(0):1})
s=sp.simplify(s); print(s)