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

np.random.seed(0)
data = np.random.randint(10, 60, size=100)
data_series = pd.Series(data)

mean_value = data_series.mean()
median_value = data_series.median()
mode_value = data_series.mode()[0]  # Get the first mode

print(f"Mean: {mean_value}")
print(f"Median: {median_value}")
print(f"Mode: {mode_value}")

plt.figure(figsize=(10, 6))
plt.hist(data_series, bins=15, color='lightblue', alpha=0.7, edgecolor='black')

# Correcting the f-string syntax by keeping it on a single line
plt.axvline(mean_value, color='red', linestyle='dashed', linewidth=2, label=f'Mean: {mean_value:.2f}')
plt.axvline(median_value, color='green', linestyle='dashed', linewidth=2, label=f'Median: {median_value:.2f}')
plt.axvline(mode_value, color='blue', linestyle='dashed', linewidth=2, label=f'Mode: {mode_value:.2f}')

plt.title('Histogram of Ages with Mean, Median, and Mode')
plt.xlabel('Age')
plt.ylabel('Frequency')
plt.legend()
plt.grid(axis='y', alpha=0.75)
plt.show()