function D = loadimage(filename)
% D = loadimage(filename)
%
% returns a column vector representation of the image read from the
% file filename

D = double(imread(filename));
D = reshape(D,prod(size(D)),1);
