Bimodal Distribution

import matplotlib.pyplot as plt
import numpy as np

COLORS = {'bars': '#8B5CF6', 'peaks': '#EF4444', 'background': '#FFFFFF', 'text': '#1E293B', 'text_muted': '#64748B', 'grid': '#F1F5F9'}

np.random.seed(42)
data1 = np.random.normal(35, 8, 400)
data2 = np.random.normal(65, 10, 600)
data = np.concatenate([data1, data2])

fig, ax = plt.subplots(figsize=(10, 6), dpi=100)
ax.set_facecolor(COLORS['background'])
fig.patch.set_facecolor(COLORS['background'])

counts, bins, _ = ax.hist(data, bins=30, color=COLORS['bars'], edgecolor='white', linewidth=0.5, alpha=0.8)

# Mark peaks
ax.axvline(35, color=COLORS['peaks'], linewidth=2, linestyle='--', alpha=0.7)
ax.axvline(65, color=COLORS['peaks'], linewidth=2, linestyle='--', alpha=0.7)
ax.scatter([35], [max(counts)*0.95], marker='v', s=60, c=COLORS['peaks'], zorder=5)
ax.scatter([65], [max(counts)*0.95], marker='v', s=60, c=COLORS['peaks'], zorder=5)

ax.spines['top'].set_visible(False)
ax.spines['right'].set_visible(False)
ax.spines['left'].set_color(COLORS['grid'])
ax.spines['bottom'].set_color(COLORS['grid'])
ax.yaxis.grid(True, color=COLORS['grid'], linewidth=1, zorder=0)
ax.set_axisbelow(True)
ax.tick_params(axis='both', colors=COLORS['text_muted'], labelsize=9, length=0, pad=8)
ax.set_xlabel('Value', fontsize=10, color=COLORS['text'], labelpad=10)
ax.set_ylabel('Frequency', fontsize=10, color=COLORS['text'], labelpad=10)

plt.tight_layout()
plt.show()