along with a successive over relaxation method so we will converge to a solution faster. The boundary conditions for the stream function will not change during the progression of the program but we will have to calculate for the vorticity function along the wall separately before we may calculate the next time step for the vorticity function inside. The formula for solving the vorticity function along the wall will involve a Taylor expansion at the boundaries and will be of the form,
\(\)