fs = 8000;
f = [0 355 (415+355)/2 415 1200 (1200+1270)/2 1270];
a = [0 0.25 0.75 1 0];
rp = 0.25;
rs = 48;

dev = [(10^(-rs/20)) (10^(rp/20)-1)/(10^(rp/20)+1) (10^(-rs/20))];

[n, f0, a0, w] = firpmord(f, a, dev, fs);
b = firpm(n, f0, a0, w);
freqz(b, 1, 1024, fs);