add subsitute function for ginac matrices