a = (E^(-((x^2)/2)))
F0 = a
F1 = a x
F2 = a ((x^2) - 1)
Plot[{F0, F1, F2}, {x, -4, 4}]   
b = (E^(-((y^2)/2)))
G0 = b
G1 = b y
G2 = 
      b ((y^2) - 1)
  
ContourPlot[F0 G0, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F0 G1, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F1 G0, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F0 G2, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F2 G0, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F1 G1, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F1 G2, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F2 G1, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 40]
ContourPlot[F2 G2, {x, -5, 5}, {y, -5, 5}, 
      PlotPoints -> 60]