added horzcat and vertcat to design_matrix()