#374 Hexbin plot with marginal distributions using seaborn library

import numpy as np
import matplotlib.pyplot as plt
#conda install -c anaconda seaborn
import seaborn as sns
sns.set_theme(style="ticks")

def main():
    rs = np.random.RandomState(11)
    x = rs.gamma(2, size=1000)
    y = -.5 * x + rs.normal(size=1000)
    sns.jointplot(x=x, y=y, kind="hex", color="#4CB323")
    plt.savefig("ex374.png", dpi=100)
    plt.show()

if __name__ == '__main__':
    main()