
For complex equations, I tend to use a symbolic approach because it allows me to verify their correctness. You should be able to write your equations accordingly. Please update your code and display the surface.
% declare symbolic variables
syms x y
% write the surface equation, z = f(x, y)
xi = @(x, y) (x + 1/4./gamma(y + 1).*2.^y)/10;
z = - 1/4*(1 + erf(sqrt(3)/2*xi(x, y)) - 1/2*exp(- sqrt(3)/2*xi(x, y).^2));
% display the equation
disp(z)
% plot the surface
fsurf(z, [-20, 20, 0, 10])
view(30, 30)
xlabel('x')
ylabel('y')
zlabel('z')