start_matrix function

Setup of the start matrix