Hi Katie,
you have to evaluate gfu at the desired points, this can be done by "gfu(mesh(x,y))".
As you suggested, it is probably the easiest to use matplotlib to plot the data.
Code:
import matplotlib.pyplot as plt
N = 100
xpnts = [i/N for i in range(N+1)]
vals = [gfu(mesh(xi,0.5)) for xi in xpnts]
plt.plot(xpnts,vals)
plt.show(block=False)
Best,
Christoph