v=linspace(-2,2,50)
[x,y]=meshgrid(v,v)
u=y+x-x*(x**2+y**2)
v=y-x-y*(x**2+y**2)
x0=[-2 ,-2, -2, -2, -0.5, -0.5, 0.5,0.5, 2,2, 2, 2, -0.01, -0.01, 0.01, 0.01]
y0=[-2 ,-0.5, 0.5, 2, -2, 2, -2, 2, -2, -0.5, 0.5, 2, -0.1, -0.1, -0.1, 0.1]
plot(x,y)
plot(x0,y0)
plot(u,v)

Kreyszig-6.9_ex (last edited 2010-12-18 10:14:53 by eshwara)