So wp is a constant and w is the frequency you want to input?
Create a vector of frequencies to sample. You probably want this logarithmically:
w = logspace(12,16,100) % 100 points between 10^12 and 10^16
Then, evaluate that frequency response for EACH frequency. Whenever you do element-by-element calculations, you want to use dot multiply/dot divide/dot exponential.
e = 1 - (wp./w).^2;
Finally, you can plot this in a log-log plot:
loglog(w,e)
- Sebastian