The code I used that ended up working was as follows:
%Fit Plot
fitobject = fit([X1,Y1],Z1,'poly11')
%Degree Polynomial ' poly11', poly'22' etc. each Number represents the degree of each variable
figure (3)
plot(fitobject,[X1,Y1],Z1)
hold on
grid on
title('title')
set(gca, 'ZLim',[0 6])
shading interp
hold off