For the initial value problem, $y'=2x(1+y),y(0)=2$, the Picard iteration is defined by $$y_{n+1}=2+\int_0^x2x(1+y_n(x))\ dx \quad (1)$$
I claimed that $$y_n=-1+3\sum_{k=0}^n\frac{x^{2k}}{k!}$$ Now I need to use mathematical induction to prove this is true. Skip the base case and for induction hypothesis, we have $$y_n=-1+3\sum^{n}_{k=0}\frac{x^{2k}}{k!}=2+\int_0^x2x(1+y_{n-1})\ dx$$ For $n+1$, then \begin{align*} y_{n+1}&=-1+3\sum^{n}_{k=0}\frac{x^{2k}}{k!}+\frac{x^{2(n+1)}}{(n+1)!} \\ &=\int_0^x2x(1+y_{n-1})\ dx + \frac{x^{2(n+1)}}{(n+1)!}\quad \text{(induction hypothesis)} \end{align*} I stuck at this point, such as, how to 'combine' these two terms, to make it become $\int_0^x2x(1+y_n)\ dx$? Alternatively, I was simply using (1), but I also not feel it will achieve mathematical induction. Any helps will be great. Thanks in advance.