Comparto mi solución,en la que se hace la rotación alrededor del eje x del área delimitada
por la recta sobre el eje x delimitada por el intervalo [0, 2], la recta sobre el eje y
delimitada por el intervalo [0, 1] o y = x + 1 evaluada en 0, la recta resultado de evaluar
x + 1 en el intervalo [0, 1.3] siendo 1.3 el valor aproximado donde se intersecta la recta
con la parábola y finalmente el límite de la curva generada al evaluar -x^2 + 4 en el
intervalo [1.3, 2].
El código traté de hacerlo claro espero que ayude.
import numpy as np
import matplotlib.pyplot as plt
import mpl_toolkits.mplot3d.axes3d as axes3d
np.seterr(divide='ignore', invalid='ignore')
fig = plt.figure()
solido = fig.add_subplot(111, projection='3d')
x = np.linspace(0, 2, 100)
xRecta = np.linspace(0, 1.3, 100) ### Intervalos de los segmentos de curva que
xParabola = np.linspace(1.3, 2, 100) ### delimitan el área. x = 1.3 intersección
para = -x**2 + 4 ### Funciones para dibujar la gráfica 2D
rect = x + 1
v = np.linspace(0, 2*np.pi, 100)
Xr, V = np.meshgrid(xRecta, v)
Yr = (Xr + 1) * np.cos(V)
Zr = (Xr + 1) * np.sin(V)
Xp, V = np.meshgrid(xParabola, v)
Yp = (-Xp**2 + 4) * np.cos(V)
Zp = (-Xp**2 + 4) * np.sin(V)
X = np.concatenate((Xr, Xp), axis=0)
Y = np.concatenate((Yr, Yp), axis=0)
Z = np.concatenate((Zr, Zp), axis=0)
solido.set(xlabel='x',
ylabel='y',
zlabel='z',
xlim = [0, 2],
ylim = [-3, 3],
zlim = [-3, 3],
title='Sólido de revolución')
grafica = fig.add_axes([0.05,0.7,0.15,.2])
grafica = plt.plot(x, para, 'r')
grafica = plt.plot(x, rect, 'b')
solido.plot_surface(X, Y, Z, cmap='viridis')
plt.show()
Anexo una imagen del resultado