In this paper, we give matrix formulae for non-orientable surfaces that provide the Laurent expansion for quasi-cluster variables, generalizing the orientable surface matrix formulae by Musiker-Williams. We additionally use our matrix formulas to prove the skein relations for the elements in the quasi-cluster algebra associated to curves on the non-orientable surface.