import seaborn as sns
import matplotlib.pyplot as plt
import numpy as np

# Load data
tips = sns.load_dataset('tips')
data = tips['total_bill']

# Calculate quartiles and Bowley skewness
Q1, Q2, Q3 = np.percentile(data, [25, 50, 75])
B = (Q1 + Q3 - 2*Q2) / (Q3 - Q1)

# Create boxplot
sns.boxplot(x=data)
plt.annotate(f'Bowley Skewness: {B:.2f}', xy=(0.05, 0.9), xycoords='axes fraction', fontsize=12)
plt.show()