p = plot(sin(x)*sin(y), sin(x)*cos(y), cos(x), (0,pi), (0,2*pi))
p[0].surface_color = lambda a, b : cos(4*a)*sin(2*b) p.show()