draw a function from double integral

See more:

MatLab

Hi!
I want to draw a two-variable function from X and Y as with below m-file:
it is under a double integral and there are some errors.
Does any body help me????

syms z h tetha0 tetha x y;
tetha = pi/4;
n = [h*cos(tetha)./sin(tetha); h*abs(cos(tetha)); 0];           %3*1 matrix
ft = ((tetha - pi/2)./sin(tetha)).^4;
Rt = [cos(tetha) -sin(tetha); sin(tetha) cos(tetha)];
zt = [cos(tetha0) -sin(tetha0); sin(tetha0) cos(tetha0)];
lt = [x;y];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h = 4;
tetha0 = 0;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

lambda = ft.*(Rt*zt*lt);                                     %result a 2*1 matrix
lm = [lambda;0];                                                %result a 3*1 mtrix
integrand = (cross(lm,n));
PhiHat = @(X,Y)(dblquad(integrand,0,X,0,Y));

ezsurf(PhiHat,[0,5]);