Matlab: como representar uma função trignométrica em modo tridimensional?

É um processo muito simples.
x=0:pi/1000:2*pi  (definição do intervalo de variação)
[x,y]=meshgrid(-4:0.1:4) (definição da malha bidimensional)
z=sin(x) (definição da função seno a ser projectada no eixo vertical)
surf(x,y,z) (instrução para fazer o gráfico tridimensional)
xlabel('x') (legenda do eixo x)
ylabel('y') (legenda do eixo y)
zlabel('sin(x)') (legenda do eixo z)

surf(x,y,z) (instrução para fazer o gráfico tridimensional)

Eis o resultado: