def TochnoeReshenie(f,a,b,n,y0):
    res=[]
    h=(b-a)/n
    ilist = np.arange(0, n, 1)
    xlist = [(a+h*i) for i in ilist]
    for x in xlist:
        func = 1/(exp(x)+x+1)
        res.append(func)
    return res