Python 符号计算：用新表达式或符号替换旧表达式

问：

我正在学习python来计算偏导数。我有功能

y = x1**2 + 2*x1 

跟

x1 = 2 *x + 1 

并希望获得 的值。diff(y, x)

编码如下;

import sympy as sp
import numpy as np  

x = sp.symbols('x')
x1 = 2*x + 1

y = x1**2 + 2 * x1

ans = sp.diff(y, x)

我得到了

ans = 8*x + 8

我想在最终结果
中使用替换，即"x1""x""4 * x1 + 4"

我该怎么做？