p1:=plot::Surface([x,y,x^2-y^2],r=0..1,t=0..a,a=0..2*PI,Mesh=[5,51])
x:=r*cos(t);y:=r*sin(t);
plot(p1);