I need a python program for Taylor series solving a differential equation dy/dx=x^2+y^2 and y(0) =1

Question

from sympy import *
x = symbols(\'x\')

def taylor(f,x,x0, n):  
    t = 0  
    for i in range(0,n):
        t = t + ((diff(f,x,i).subs(x,x0))/(factorial(i))) * (x