EEGLAB: How to plot the power spectrums for all channels or a subset of channels in a single graph across a frequency range?

Question

I'm an undergrad researcher using EEGLAB. I want to create a single plot that contains specific channels to show which channels correlate with the highest alpha power.

I tried finding the code used for the pop_prop function in EEGlab but I wasn't sure what to change in the code in order to produce a single plot that shows the power spectrums for a subset of 20 channels across the frequnecy range [2, 50].