:

#
x, y = np.meshgrid(np.linspace(-2,2,20), np.linspace(-2,2,20))
z = x**2 + y**2
ashaxiom.Axiom.start(z).plot3D(x, y)