how to remove 60hz noise from an ECG signal using notch filter in matlab

The code for notch filter is ready

 >>>>>>>>>>>

fs=300;
f0=60;

r=0.99;

b0=1;
b1=-2*cos(2*pi*f0/fs);
b2=1;
a0=1;
a1=-2*r*cos(2*pi*f0/fs);
a2=r*r;
b=[b2,b1,b0];
a=[a2,a1,a0];
[h,w]=freqz(b,a);
hf=abs(h);
figure(1);
plot(w*fs/(2*pi),hf);
>>>>>>>>
the ecg data is also available..

How to interface these notch filter and given ECG data ?