return val black = [0,0,0] red = [255,0,0] green = [0,255,0] blue = [0,0,255] w = 600 h = 600 # display list of line segments from contouring routine gr = racg.Graphics(w=600,h=600) gr.Clear() s = 0.5*10**8 k = 4 x0,y0 = 200,100 gr.Clear() G = lambda s: lambda x,y: F(x,y)-s ix = 100 iy = 100 for y0 in range(1,600,iy): for x0 in range(1,600,ix): ctr.ShowLevelSet(gr,G(s), x0, y0, blue, tmax = 50) gr.Point([x0,y0],red) ch = gr.Show("result",15) if ch == ord('e'): break gr.Show("result",-1) fn = "ImplicitFunction-ContourVelocity-%02d.jpg" % k #gr.Save(fn) gr.Close()
black = [0, 0, 0] red = [255, 0, 0] green = [0, 255, 0] blue = [0, 0, 255] w = 600 h = 600 # display list of line segments from contouring routine gr = racg.Graphics(w=600, h=600) gr.Clear() s = 0.5 * 10**8 k = 3 x0, y0 = 200, 100 gr.Clear() G = lambda s: lambda x, y: F(x, y) - s ix = 50 iy = 50 for y0 in range(1, 600, iy): for x0 in range(1, 600, ix): ctr.ShowLevelSet(gr, G(s), x0, y0, blue) gr.Point([x0, y0], red) ch = gr.Show("result", 15) if ch == ord('e'): break gr.Show("result", -1) fn = "ImplicitFunction-ContourVelocity-%02d.jpg" % k #gr.Save(fn) gr.Close()