Probability Density

import matplotlib.pyplot as plt
import numpy as np
from scipy import stats

# Data - Multiple distributions
x = np.linspace(-4, 8, 200)
dist1 = stats.norm.pdf(x, 0, 1)
dist2 = stats.norm.pdf(x, 3, 1.5)
dist3 = stats.gamma.pdf(x, 3, scale=1)

# Figure - LIGHT THEME
fig, ax = plt.subplots(figsize=(10, 6), facecolor='#ffffff')
ax.set_facecolor('#ffffff')

# Probability densities with NEON colors
ax.fill_between(x, 0, dist1, alpha=0.4, color='#27D3F5', label='Normal(0,1)')
ax.fill_between(x, 0, dist2, alpha=0.4, color='#F5276C', label='Normal(3,1.5)')
ax.fill_between(x, 0, dist3, alpha=0.4, color='#6CF527', label='Gamma(3,1)')
ax.plot(x, dist1, color='#27D3F5', linewidth=2)
ax.plot(x, dist2, color='#F5276C', linewidth=2)
ax.plot(x, dist3, color='#6CF527', linewidth=2)

# Styling
ax.set_xlabel('x', color='#1f2937', fontsize=11)
ax.set_ylabel('Probability Density', color='#1f2937', fontsize=11)
ax.set_title('Probability Density Functions', color='#1f2937', fontsize=14, fontweight='bold', pad=15)
ax.tick_params(colors='#374151', labelsize=9)
for spine in ax.spines.values():
    spine.set_color('#e5e7eb')
ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.yaxis.grid(True, color='#f3f4f6', linewidth=0.5)
ax.set_xlim(-4, 8)
ax.legend(facecolor='#ffffff', edgecolor='#e5e7eb', labelcolor='#1f2937')

plt.tight_layout()
plt.show()