Domain coloring example had to change from fig = plt.figure() to fig, ax = plt.subplots() to make it work
jwkvam
