n = size(0:.01:pi*2);
n = n(2);
z = ones(n);
aux1 = 0;
aux2 = 0;
for i = 1:n
    for e = 1:n
        aux1 = aux1 + .01;
        z(i,e) = sin(aux1) * cos(aux2);
    end;
    aux2 = aux2 + .01;
end;
surf(1:i, 1:e, z);
shading interp;