x, y, z, t = var('x, y, z, t')
f(x,y)=(x^2)*(y^2)+3*x*(y^2)-2*y
P=plot3d(f, (x,-3, 3), (y,-3, 3), color='red')
Q=implicit_plot3d(y==2, (x,-4,4), (y,-3,3), (z,-26,168))
R=parametric_plot3d((t, 2, 4*t^2+12*t-4), (t, -4, 4), color='black', thickness=10)
S=point3d((1,2,12),size=40,color='black')
show(P+Q+R+S,figsize=(6,6))