How do I devise bandpass chebyshev filter ?

I designed low pass chebyshev filter n=6 and ripple 0.1-dB Ripple (ε = 0.15262) -> (s2 + 0.22939s + 1.12939) (s2 + 0.62670s + 0.69637) (s2 + 0.85608s + 0.26336).

H= (1.12939/((〖(i*w)/1000)〗^2+(0.22939 *i*w)/1000+1.12939) * 0.69637/((〖(i*w)/1000)〗^2+(0.62670*i*w)/1000+0.69637)*0.26336/((〖(i*w)/1000)〗^2+(0.85608*i*w)/1000+0.26336))

H= (1129390/(〖(i*w)〗^2+229.39*i*w+1129390) * 696370/(〖(i*w)〗^2+626.70*i*w+696370)*263360/(〖(i*w)〗^2+856.08*i*w+263360))

But I also designed bandpass filter via chebyshev.

How do I add extra things to my code included in low pass chebyshev filter?

Here is my code below.

wn=0;
w=0;
wi=0;
H=0;
Ha=0;
for w=1:2000
wi(w)=w;
wn=w;
H(w)=(1129390)*(696370)*(263360)/((-wn*wn + 229.39*i*wn + 1129390)*(-wn*wn + 626.70*i*wn + 696370)*(-wn*wn + 856.08*i*wn + 263360));
Ha(w)=abs(H(w));
plot(wi,Ha);
end

[Daniel kiracofe]

Why not just use cheby1() or cheby2() functions?