How to implement LH, HL, HH in DWT of an image ?

I want to implement the subbands without using the built in functions in Matlab

I have calcaluted the LL only but don't know what should be changed to calculate LH, HL, HH subbands

what I have reached so far..

X=double(imread('image.jpg'));
[row,col]=size(X);
Scale_filter = [1/sqrt(2) 1/sqrt(2)]; %Haar scaling filter
Y = conv2(X,Scale_filter(:)','valid'); %filter along the columns
Y(:,end+1) = Y(:,end); % periodically extend the matrix for the next step
Y = Y(:,2:2:end); % downsample along the columns
Y(end+1,:) = Y(end,:); %periodically extend the rows for the next step
Y = conv2(Y',Scale_filter(:)','valid'); % filter along the rows
Y = Y'; 
Y = Y(1:2:end,:)  %downsample