function [wll,wlh,whl,whh]=get_2DWavelet_coif1
[Lo_D,Hi_D,Lo_R,Hi_R]=wfilters('coif1');
L=[Lo_D 0];
H=[0 Hi_D];
Nw=size(L,2);
xl=zeros(Nw,Nw);
xh=zeros(Nw,Nw);
yl=zeros(Nw,Nw);
yh=zeros(Nw,Nw);
zl=zeros(Nw,Nw);
zh=zeros(Nw,Nw);
for i=1:1:Nw
    xl(:,i)=L;
    xh(:,i)=H;
    yl(i,:)=L;
    yh(i,:)=H;
end
wll=xl.*yl;
wlh=xl.*yh;
whl=xh.*yl;
whh=xh.*yh;