scores <- c(55, 67, 72, 75, 78, 80, 82, 85, 88, 90, 
            92, 95, 60, 68, 74, 79, 81, 83, 86, 91)

hist(scores,
     main = "Distribution of Student Scores",
     xlab = "Scores",
     ylab = "Frequency",
     col = "skyblue",
     border = "black",
     breaks = 7)