ClearAll[f];
f[x_,y_]:=(1-y)(x^2)+2x*y;
Show[Plot3D[f[x,y],{x,-2,2},{y,-1,3},PlotRange->{All,All,{-25,8}},Axes->True,AxesLabel->{x,y,z},Mesh->None,Lighting->Automatic,
RegionFunction->Function[{x,y},(x^2)-1≤y≤3]],Graphics3D[{PointSize[0.02],Blue,Point[{3/2,3,f[3/2,3]}]}],
Graphics3D[{PointSize[0.02],Red,Point[{-2,3,f[-2,3]}]}]]