2D case, diverge
 

I have made the code by vorticity and stream function equations with G-S in 2D steady case. I only can get the value for 20 steps.after that float outflow. How could I sovle this problem? many thanks