# Perform KDE
xMin, xMax = x.min(), x.max()
yMin, yMax = y.min(), y.max()
print(xMin, xMax, yMin, yMax)
kde = gaussian_kde([x, y], bw_method='silverman')
xGrid, yGrid = np.meshgrid(np.linspace(xMin, xMax, 100), np.linspace(yMin, yMax, 100))
density = kde(np.vstack([xGrid.ravel(), yGrid.ravel()])).reshape(xGrid.shape)
# Plotting
plt.figure(figsize=(10, 8))
plt.imshow(np.sqrt(density), origin='lower', aspect='auto', extent=[xMin, xMax, yMin, yMax])
plt.colorbar(label='Density')
plt.show()