h=1;
x=[0:0.1:h];
y=[0:0.1:h];
u=(1-(x/h))'*(1-(y/h));
surf(u)
title('Pagoda Function')