solve initial value problem x2y'-xy=x2+4,
x>0, y(1)=0