c
a*np.exp(-((x+x_m)**2/(2*std_x**2)+(y-y_m)**2/(2*std_y**2)))