The aoctool() only allows for one grouping variable. Instead, anovan() might be an option for what you're describing.
anovan() is an n-way analysis of variance allowing for multiple grouping variables. When at least one of those grouping variables are labeled as a continuous variable, an ANCOVA is performed (Matlab calls it ANOCOVA, same thing).
