% Create matrix indices for TIPSH2D.m

function [rowind, colind] = indices(j,k1,k2)

rowind = [(2^j)*k1+1:(2^j)*(k1+1)];
colind = [(2^j)*k2+1:(2^j)*(k2+1)];

%
%  This code has been developed by Eric Kolaczyk, 
%  The University of Chicago, 1996.
%  Copyright (c) Eric Kolaczyk, 1996.
%