s: streamfunction
```

```
def rhs(nx,ny,dx,dy,re,w,s):
aa = 1.0/(dx*dx)
bb = 1.0/(dy*dy)
gg = 1.0/(4.0*dx*dy)
hh = 1.0/3.0
f = np.zeros((nx+3,ny+3))