SymPy проблемы при работе с функцией

import numpy as np
import matplotlib.pyplot as plt
from sympy import *
x = symbols('x')
 
dots1 = [(2, 12.258), (4, 17.24), (8, 30.151), (3*x+5)]
dots2 = [(2, 3.688), (4, 10.791), (8, 20.705), (0.25*x**2+0.75*x+1.25)]
dots3 = [(2, 4.872), (4, 29.707), (8, 246.971), (10, 485.727), (12, 840.658), (0.5*x**3-0.25*x**2+0.75*x+1.25)]
all_dots = [dots1, dots2, dots3]
 
for dots in all_dots:
    fig, ax = plt.subplots()
    print("dots", dots[:-1])
    f = dots[-1]
    print('function', f)
    summ=0
    count=0
    for dot in dots[:-1]: # сумма квадратов разницы между истинным значением и предсказанным
        summ += (dot[1]-f.subs(x, dot[0]))**2
    MSE = summ/(len(dots)-1)
    
    #соберём значения x и y заданных точек
    x0, y0 = [x for x, y in dots[:-1]], [y for x, y in dots[:-1]]
    plt.scatter(x0, y0, c='r') # отрисуем точки
    
    # зададим область определения функции (от -10 до 10 с шагом 200)
    x = np.linspace(-10, 10, 200) 
    plt.plot(x, f) # отрисуем функцию
    print('RMSE =', sqrt(MSE), '\n')
    plt.show()

f=3*x+5
x = np.linspace(-10, 10, 200)
plt.plot(x, f)
plt.scatter(x0, y0, c='r')

    xx = np.linspace(-10, 10, 200) 
    
    y=[f.evalf(subs={x:it}) for it in xx]
    for a,b in zip(xx,y):
        plt.plot(a, b) # отрисуем функцию
    print('RMSE =', sqrt(MSE), '\n')
    plt.show()

for dots in all_dots:
    #fig, ax = plt.subplots()
    x = symbols('x')
    print("dots", dots[:-1])
    f=dots[-1]
    print('function', f)
    
    summ=0
    count=0
    for dot in dots[:-1]: # сумма квадратов разницы между истинным значением и предсказанным
        summ += (dot[1]-f.subs(x, dot[0]))**2
    MSE = summ/(len(dots)-1)
    
    #соберём значения x и y заданных точек
    x0, y0 = [xx for xx, yy in dots[:-1]], [yy for xx, yy in dots[:-1]]
    plt.scatter(x0, y0, c='r') # отрисуем точки
    
    # зададим область определения функции (от -10 до 10 с количеством значений 200)
    xn = np.linspace(-15, 15, 200) 
    # Найдём значения функции на определённой области
    y=[dots[-1].evalf(subs={x:it}) for it in xn]
    plt.plot(xn, y) # отрисуем функцию
    print('RMSE =', sqrt(MSE), '\n')
    plt.show()