#---------------------------------------------------------------
#
#      Sample code
#                ----------- The Riemann surface for log(z)
#
#  think x as the  radius and y as the argument
#---------------------------------------------------------------

zz(x,y) = x *exp(i * y)

logzX(x,y) = real(zz(x,y))
logzY(x,y) = imag(zz(x,y))
logzZ(x,y) = 0.2*y

logz = sur{
         [logzX(x,y),logzY(x,y),logzZ(x,y)]
         [x=0:1][y=-13:13] [samp=20:100]
        };

#---------------------------------------------------------------
plot logz;
#---------------------------------------------------------------
